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ABSTRACT 


The usual procedure in structural analysis involves the deter- 
mination of the strength of a certain structure under specified 
loading conditions. From the design point of view the real 
problem is to determine the lightest practical arrangement of 
material which will transmit the required loads through the speci- 
fied distances. This may be accomplished through the use of the 
structural index. The engineering significance of this index and 
its applications to design problems are discussed. Examples 
include the tubular column, stiffened flat plate under axial load, 
the unstiffened plate with edge support, shear-carrying members, 
and cylindrical shells in bending. Methods of comparing various 
materials on a weight basis are developed and applied to several 
aluminum alloys, magnesium alloy, and stainless steel. 


SYMBOLS 


= area of cross section 
width of a box beam or column 
span; width of a panel 
coefficient of fixity for columns 
maximum distance from neutral axis 
= diameter of a circular shell or tube 
modulus of elasticity (Young’s modulus) 
tangent modulus (rate of change of stress with strain) 
stress 
= equivalent stress (actual stress divided by density rela- 
tive to aluminum alloy) 
= depth of box beam 
moment of inertia 
= constant 
kilopound (= 1,000 Ibs.) 
kips per sq.in. (= 1,000 Ibs. per sq.in.) 
actual column length 
= effective (pin-ended) column length = L/+/¢c 
moment (bending or torsional) 
= force or load (Ibs.) 
pressure 
= load per unit length (lbs. per in.) 
Presented at the Annual Summer Meeting, I.A.S., Los Angeles, 
July 15-16, 1948. 
* This paper is based on a report prepared for The Rand 
Corporation (No. RA-15069, dated February 27, 1948). 
t Consultant and Lecturer, respectively. 


S = wing area 

t thickness 

t average thickness (area per unit width) 
W weight (Ibs.) 

w specific weight (Ibs. per cu.in.) 

p radius of gyration (= +»/J/A) 

T tangent modulus ratio (= E,/E) 

# = Poisson’s ratio 


(1) Optimum DESIGN PRINCIPLES 


| toe ANY structure that fails as a result of instability 
under compressive (or shear) loading, the allow- 
able ultimate stress (average P/A) depends in a complex 
manner on the properties of the material and the geo- 
metrical proportions of the structure. For a given 
structure, such as a stiffened sheet panel, it is possible 
to determine the average failing stress by a direct test 
or to estimate it by theoretical methods. This is es- 
sentially the approach used in stress analysis—i.e., the 
problem starts with a known or assumed shape and 
size of structure for which the allowable stress is to be 
determined. 

It is possible to approach this problem from the 
design point of view by starting with the load, the free 
length over which it must be transmitted, and the 
length over which the material must be distributed. 
For any given material and configuration (proportions 
of the cross section) the allowable stress can be ex- 
pressed as a function of these design variables. Further- 
more, if a number of different types of structure are in- 
vestigated on this basis, it becomes possible to observe 
which type develops the highest allowable stress and 
therefore represents the optimum design from the stand- 
point of minimum structural weight. Finally, an en- 
velope curve can be obtained which will represent the 
best possible allowable stress for any given loading 
conditions and for a given material. If such a curve is 
used for weight analysis, the result is a basic weight that 
represents the absolute minimum—i.e., the minimum 
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based on the knowledge available at the time—for the 
material in question and for all types of structure 
covered by the investigation. This portion of the total 
structural weight is governed directly by the design 
variables, and it should be possible to derive formulas 
or charts that accurately predict the effects of changes 
in these variables, over any desired range. The effects 
of secondary variables must be handled separately. 

This approach is analogous to the situation in aero- 
dynamics, where the induced drag represents the abso- 
lute minimum for a given weight, span, and air speed, 
while the remaining (parasite) drag is controlled by 
other factors and, in fact, responds to changes in air 
speed in an entirely different manner. It will be shown 
later, for example, that the weight of a shell loaded in 
bending responds to changes in diameter in a manner 
similar to the relationship between wing drag and air 
speed—i.e., a portion of the structural weight decreases 
with increasing diameter while another portion increases. 
As a result there is an optimum diameter for a given 
bending moment. 


(2) THe STRUCTURAL INDEX 


As previously noted, the allowable stress can be ex- 
pressed as a function of load and distance. References 
1, 2, 3, and 4 contain examples of quantities that are 
used as parameters in investigating different types of 
loading conditions. Cox and Smith? use the term 
“structure loading coefficient” suggested by de Bruyne.’ 
Goodier and Thompson* propose the term “strain 
number.’ Zahorski‘ uses the term ‘‘structural index,” 
originally proposed by Wagner." 

Since the terms ‘‘coefficient” and ‘‘number’ 
dimensionless quantities and since it will be found de- 
sirable from an engineering point of view to work with 
dimensional quantities, the term ‘structural index’”’ 
is used here and is suggested as standard terminology. 

The structural index is a measure of the loading in- 
tensity; it should therefore have the dimensions of a 
stress. For a known load, P, the index has the form 
P/ab, where a and b are two dimensions of the structure. 
In the case of the simple column, it can be shown 
(references 2 and 7) that both a and 6 take on the value 
of the effective column length, so that the structural 
index becomes P/Lo?. It will be shown later that for 
the ‘“‘wide column’’ the index becomes P/BL», where B 
is the width of the column. For a circular shell in 
bending, a suitable index is 1//D*. 

The great value of the structural index is that it elimi- 
nates the effects of size in dealing with allowable 
stresses. Once a structural element of certain size has 
been designed (or tested) for a given load (P;), its 
proportions may be used for any other combination of 
load (P:) and size which results in the same value of the 
structural index. It is only necessary to increase all 
dimensions by a multiplying factor, which has the value 


V P;/P;. Since all dimensional ratios remain unchanged 


imply 


in such an expansion process, the allowable stress wil] 
remain unchanged. This holds true even in the inelastic 
range, because for a given material both the tangent 
modulus and the secant modulus are functions of the 
stress, which does not change. (Note: This type of 
analysis assumes that the stress-strain diagram is not 
affected by size, which is not strictly true. The effects 
of strain gradient on the allowable stress are also ig- 
nored. Both of these effects are of secondary impor. 
tance and need not be considered in a general weight 
analysis. For extremely large structures, however, a 
supplemental investigation of size effects would appear 
to be in order.) 

By making use of dimensional similarity principles 
in this manner, it is possible to extend available test 
data to cover a wide range of sizes. Conversely, it is 
possible to condense a large mass of data, plot the data 
on one chart, and thereby obtain an envelope curve for 
the optimum allowable stress as a function of the struc- 
tural index. For each type of loading there will result 
an optimum curve for each different material investi- 
gated. From such curves the effects of using different 
materials can be evaluated. 

The curves of allowable stress, plotted against struc 
tural index, not only give an envelope for the optimum 
structure but also reveal how various types of structure 
compare on an efficiency basis. This enables studies 
to be made to determine the weight penalties involved 
in using a ‘“‘nonoptimum’’ type of structure. 

Although test data are usually used as a basis for the 
allowable stress curves, it is possible in some cases to 
express the relationships mathematically. The simple 
column can be handled in this manner and is used as 
an introductory example. 


(3) MrntmumM WEIGHT FOR SIMPLE COLUMNS 


The pin-ended column represents the simplest case 
in which the allowable stress is affected by magnitude 
of the applied load and the length of the transmission 
path. It is used as an example in presenting the 
general methods by which all types of structures may be 
analyzed on a minimum-weight basis. 

The round tube is used as a basis for the investiga- 
tion, since it is possible to express its strength in mathe- 
matical form. It should be noted, however, that the 
method can be based directly on test data for cases in 
which analytical expressions are either not available or 
are considered to be too inaccurate.’ 

The design parameters in this case are the load P and 
the length L, through which the load must be trans- 
mitted. Fora pin-ended column the allowable stress for 
primary buckling failure may be obtained from the 
tangent-modulus formula: 


f = wE,/(L/p)? (1) 


(NoTE: The analysis may be applied to columns 
having end restraint by replacing LZ by Ly.) 
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CURVES ARE 
TYPICAL STRESS-STRAIN DIAGRAMS 
AND SHOULD NOT BE USED EOR 
MINIMUM GUARANTEED VALUES 
E STEEL * 263 210 
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Fic. 1. Tangent modulus curves for typical materials. 


Since the weight will be a minimum when the allow- 
able stress, f, is a maximum, it is desirable to show how 
the allowable stress varies as a function of P, L, and 
some shape parameter that is under the control of the 
designer. In Eq. (1) the tangent modulus, £,, is directly 
related to the stress, f, for a given material, and it is 
assumed that this relationship is known. For the 
purposes of this report the tangent modulus curves 
shown in Fig. 1 are used. These represent typical ma- 
terials and are not to be used as a basis for stress analy- 
sis reports for which minimum guaranteed values are 
required. 

Eq. (1) contains one of the design variables (ZL). 
The term p? is governed by the shape of the cross sec- 
tion, but it cannot be.used for the shape parameter be- 
cause it is dimensional. It is necessary, therefore, to 
replace p® by a nondimensional factor that depends on 
the proportions of the cross section. 

Since p* has the dimension of an area, it may be 
divided by the cross-sectional area of the column to 
obtain the shape parameter, ky. 


ki = p?/A . (2) 
The area, A, may be expressed in terms of load and 
allowable stress: 

A = P/f 
Eq. (1) may now be written: 

f = wWE ki P/L*f 
from which 
f= eV ENEVP/L? (3) 

0 
P/L? = f?/wE ky (4) 


In order to show the direct effects of the tangent 
modulus, E,, it will be used in the development of all 
equations, as above. In applying the final equations, 
however, it is convenient to use the tangent-modulus 
ratio: 


r= E/E (5) © 


All final equations are therefore repeated in this 
form. Egs. (3) and (4) become: 


f= aVVEV RV P/L? (3a) 
P/L? = f?/x*rEk; (4a) 


The term (P/L?) is the structural index for the case 
of primary column failure. Together with the material 
properties and the shape parameter (f;), it controls the 
allowable stress and therefore the structural efficiency 
of the column. 

Because of the fact that E, (or r) depends on the value 
of f, Eq. (3) would have to be solved by a “cut-and-try” 
procedure. It is therefore preferable to assume values 
for f and to solve for P/L? by using Eqs. (4) and (4a). 

For thin-walled tubes, the factor k; is given by 


ky = p*/A = (1/8r)(D/t) (6) 


where D = average diameter and ¢ = wall thickness. 
Eq. (3) then becomes, for the tube, 





f = rVEN (1/8) (D/t)V P/L? (7) 
or 
P/L? = 8f?/nE{D/t) (8) 


Eq. (7) shows that, for a given value of the struc- 
tural index, the allowable stress increases as the value 
of (D/t) increases (the increase will be proportional to 
the square root of (D/t), up to the stress at which E, 
begins to decrease). For primary failure of the column 
there is obviously no upper limit for (D/t), and it is 
therefore necessary to investigate other forms of failure 
that would be involved as (D/t) is increased. 

The next type of failure that must be considered is 
local buckling (‘‘crushing,” ‘‘crippling,’”’ etc.). The 
theory of buckling of plates and shells in the inelastic 
range has not yet been developed to the same degree 
as has that for the simple column. The approximate 
methods used in this report are based on Timoshenko’s 
effective modulus (reference 5, p. 442) and are con- 
sidered satisfactory for the present purposes, which are 
primarily to develop methods and to indicate the man- 
ner in which different materials affect the structural 
weight. For actual design purposes the basic curves 
should be determined from test data whenever any 
doubt exists as to the accuracy of the equations. 

By substituting the tangent modulus for the re- 
duced modulus used in reference 5, an equation for 
local buckling of a round tube is obtained: 


For. = FalV EE,/(D/t)] (9) 


This shows that the allowable stress for crippling 
failure will decrease as (D/t) is increased. Hence, there 
will be an optimum value of (D/t) at which the allow- 
able stress is a maximum. This is found by equating 
the crippling stress [Eq. (9)] to the column stress (Eq. 
(7)] and solving for D/t, which gives the equation: 
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Fic. 2. Optimum diameter/thickness ratio for round tubes used 
as columns. 


(D/t)opt. = 2[kx*E/a(P/L*)]" (10) 


It is interesting to note that in the derivation of Eq. 
(10) the term £, canceled out, leaving only the elastic 
modulus E. This means that for all alloys having the 
same base material the optimum (diameter/thickness) 
ratio is a function of the structural index only. Assum- 
ing a value of 0.40 for k2, Eq. (10) has been evaluated 
and plotted for various classes of alloys in Fig. 2.* 

The maximum stress that can be developed in a tu- 
bular column may be found by substituting the value 
for optimum D/t in Eq. (8), giving 


P/L? = 8f3/rk. EE,” (11) 
which may be expressed as: 


P/L? = 8f3/rk.E?r'”* (11a) 


For a given material £, (or 7) is a known function of 
f, and Eq. (11) may be plotted as a curve of f against 
P/L*. Such curves are shown in Fig. 3 for several 
different materials, assuming k; to be equal to 0.40. 
(To contract the scale, the square root of the structural 
index has been used for plotting.) To facilitate weight 
comparisons, the allowable stresses have been reduced 
to a common weight basis by considering aluminum 
alloy as the standard. For this purpose the equivalent 
stress is used, as defined by the equation: 


Sea ond f/(w/wWa) (12) 


where w = specific weight of the material and w, = 
specific weight of aluminum alloy (taken as 0.10 Ib. 
per cu.in. in this report). 

If a given value of P/L? is divided by the corre- 
sponding value of f and multiplied by the specific weight 

* The theoretical value for kz, given by the classic buckling 
theory, is 1.212, but test data indicate that much lower values are 
actually developed at failure (see reference 5, p. 462). This has 
been explained by Tsien,®)!° who showed that a lower “buckling 
load” is developed when imperfections or disturbances of finite 
magnitude exist. A value of k, equal to 0.40 appears to be suit- 
able for present purposes, The extension of this theory into the 
inelastic range is to be regarded as an engineering approximation. 


1949 


(w), the resulting quantity will equal the total weight 
of the column divided by the cube of the column length 
W m wAL wA w(P/L?*) 


an 7s Ca 7 ™ 





Curves showing this relationship are plotted in Fig. 4. 
Eq. (13) obviously applies only to pin-ended columns 
To use it for columns having end restraint, the term 
W/L* must be modified as follows: 


We _ We _ w(P/L?) 
LoL Lt ff 


where c = coefficient of fixity. 

By substituting specific values for f and w in Eq. 
(13), a family of straight lines is determined on Fig. 4 
(the curves for aluminum alloys are shown). It can be 
seen that the minimum weight curves tend to become 
straight lines as the allowable stress approaches a con- 
stant value. The effect of a variable allowable stress 
is to produce a sort of “hump” at low values of the 
structural index. 

The same general methods may be used for any type 
of cross section, such as an I-section. The shape pa- 
rameter, ki; = p*/A, can be found for any given section, 
and its variation with the proportions of the section 
can be determined. The influence of the cross-section 
proportions on the local buckling or crushing stress is 
less easily determined than in the case of round tubes; 
in fact, test data would probably be necessary for reli- 
able results. 
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columns. 
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Fic. 4. Minimum weight for round tubes used as columns, 


It is important to note that this type of analysis is 
fundamentally “‘three-dimensional’’—that is, the re- 
sults can be used for members of any size provided that 
all dimensions (length, width, and thickness) remain 
proportional. 

It can also be seen that the optimum design will have 
a cross section that is symmetrical about the centerline 
of the column—i.e., the value of J (or p) must be the 
same about any axis through the center of the cross 
section. The round tube not only meets this require- 
ment but also provides the maximum possible value of 
o°/A. It will therefore be the best possible cross sec- 
tion as far as primary buckling is concerned. At ex- 
tremely low values of P/L?, however, the local buckling 
stress tends to predominate, especially for high-density 
materials, and it is therefore desirable to modify the 
design in such a way that the effective thickness is in- 
creased. (Longitudinally corrugated pipe is an example 
of this.) 


(4) MInImuM WEIGHT FOR WIDE COLUMNS 
(FREE EpGEs) 


In the shell type of structure, such as is used for air- 
plane wings and fuselages, a large portion of the ma- 
terial must resist axial loads caused by bending. If 
this material is also used to form the shell, the column 
must be spread out over a considerable width, which, 
for a given design, remains constant. Since this width 
is usually many times the thickness of the column, 
buckling can occur only in a direction normal to the 
plane of the column. If it is further assumed that the 
edges are unsupported or that the effect of such sup- 
Port is negligible, the bending stiffness of the column 
across its width may be neglected. The term “column” 


is used for members that fall in this class. The term 
‘“‘panel’”’ is restricted (as in reference 2) to members for 
which the buckling strength depends to some extent 
on the lateral bending stiffness. The wide column 
may be thought of as a series of individual columns 
placed side by side and equally loaded. 

Since the width of the column is assumed to be con- 
stant, the total end load may be divided by the width 
to obtain a unit load (load per unit width). The analy- 
sis can then be made on the basis of a column of unit 
width, using the same general method employed for the 
simple column. There is one important difference, 
however. Since the width of the column is held con- 
stant, dimensional expansion is restricted to depth and 
length only. The cross-sectional area of the column of 
unit width is equal to the total area divided by the 
width of the column, which may be thought of as an 
“average thickness,” (‘). The shape parameter is now 
obtained from the relationship 


k = p l (14) 


where 7? = A/B, A = cross-sectional area of column, 
and B = width of column. 

In reference 2 the analysis of the wide column was 
carried out for a solid section, a series of tubes placed 
side by side, and various other configurations. Refer- 
ence 1 applied this type of analysis to the flat sheet 
stiffened with stringers of Z-section. Eqs. (3) and (4) 
for the simple column are replaced by the following 
equations for the column of constant width, constder- 
ing primary buckling failure only: 


f = 2(E,)'(ks)(q/L)”" (15) 
or = 
/L = f?/nV Ecks (16) 


where g = P/B = load per unit width. 

If Eq. (15) or (16) is applied to a solid flat plate used 
as a wide column (constant width), the shape param- 
eter becomes: 


ks = p/t = 1 ‘4/12 


In the elastic range, the effective value of E is in- 
creased, for a wide plate, to the value E/(1 — yu’). 

Assuming that this effect carries over into the in- 
elastic range and that 4 = 0.30, Eqs. (15) and (16) 
become, for the flat plate: 


f = 0.97E,/(q/L)*” (15a) 
q/L = 1.05(f'”*/E,”) (16a) 
Using the alternative notation, 
f = 0.97E'“r'"(q/L)”* (15b) 
q/L = 1.05( f'*/E'r'”) (16b) 


These approximate equations are probably on the 
safe side, since it is known that the value of yu increases 
as the stress exceeds the proportional limit. (Refer- 








138 JOURNAL OF THE AERONAUTICAL SCIENCES—MARCH, 1949 


LOCKHEED DATA 


24 SRT J STIFFENER 24ST SKIN 
24STJ STIFFENER 24ST SKIN* 





75 ST Y STIFFENER 75 ST IN NACA 1389 
24 ST Y STIFFENER 24ST SKIN NACA TN 1389 


24 ST ENVELOPE 


24 SRT COR 
24 ST SKIN 


24ST CORRUGATION 24ST SKIN® 


24ST Z STIFFENER 24ST SKIN NACA TN 1439 
24ST SL STIFFENER 24ST SKIN NACA TN 1439 
24ST STIFFENER MG FS-ih SKIN NACA TN +274 


24ST Z STIFFENER 24ST SKIN NACA ARRLSFI5 


24ST PLATE AS EULER 


ie) 100 200 300 400 500 600 700 800 


4 
L> Psi 
Optimum (maximum) stress for wide columns (sheet-stiffener type). 


Fic. 5. 


ence 10 shows that the increase is not very great in the 
range of stress usually involved in structures.) 

It is interesting to note that the equations for the 
structural index of the column of constant width do not 
involve the same combination of material characteris- 
tics as do those for the simple column. The relation- 
ships are as follows: 
f(f/Ev) 

I(f/E)”” 

This would indicate that, other things being held 
constant, a change of the material properties (involving 
a change of f/E,) does not have so much effect on the 
constant-width column as on the simple column. 

The structural index also changes for the constant- 
width column, becoming g/L instead of P/L*. Since 
q = P/B, the structural index may be written as P/BL, 
where B is the width of the panel. This shows that 
when the width is held constant the length does not 
influence the allowable stress so much as in the case of 
the simple column. 

The physical significance of the index P/BL is of 
special interest. For a given load, P, the same allow- 
able stress could be developed for any combination 
of width and length which would have the same prod- 
uct. Thus, a long narrow column could have the same 
allowable stress as a short wide column. For a given 


Simple Column 
Const.-Width Column 


bending moment, a box beam of constant depth there- 
fore requires closer rib spacing as the width is increased, 
if the allowable stress is to be held constant. 


Eq. (16) could be plotted as a family of curves for a 
given material, each curve corresponding to a different 
value of k3. The highest stress would be obtained by 
using the maximum possible value of k3. As in the case 
of the tubular column, the limiting value of k3 is reached 
when some other type of failure begins to occur. This 
usually takes the form of local buckling of some por- 
tion of the structure 


As pointed out in reference 7, it is difficult to deter- 
mine a satisfactory theoretical relationship between the 
local buckling stress and the value of the shape param- 
eter, ks. Test data are needed for this purpose and 
may be obtained from ‘‘short column’”’ tests. Reference 
1 gives a complete treatment of the wide column con- 
structed from flat sheet and Z-section stiffeners, re- 
sulting in charts of allowable stress plotted against the 
structural index for various proportions of the cross 
section. 


It is possible to plot in this manner the allowable 
stress for any form of cross section and for any material, 
assuming that necessary test data are at hand. Such 
curves will indicate the optimum type of structure for 
any given value of the structural loading coefficient, 
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g/L (or P/BL). A few such curves are shown as an 
example on Fig. 5. 

Fig. 5 shows that there is a definite tendency toward 
an envelope curve, A, for “conventional’’ stiffened 
panels made from 24 ST sheet and stringers. Data 
obtained from various sources and representing differ- 
ent types of stringers were used in this study. Curve 
B, from reference 11, shows a definite increase in allow- 
able stress obtained by using a special type of stringer 
(extruded Y) developed by the N.A.C.A. Curve C 
shows how the change from 24 ST to 75 ST material 
raises the allowable stress envelope. 

On this chart there is also shown a curve, obtained 
from Eq. (16a), for the flat plate used as a wide column. 
The curve is shown for the elastic (Euler) case, which 
also covers aluminum alloys 24 ST and 75 ST for the 
loading range shown on the chart. It is obvious that 
at low values of the structural index, the flat plate, 
used as a column, is inefficient. For low values of 
loading, g, a small column length, Zo, would be required 
in order to raise the structural index into the region 
where relatively high allowable stresses can be de- 
veloped. 

Although the points on a chart such as Fig. 5 are 
obtained from tests of specimens having certain values 
of sheet thickness, stringer depth, length, etc., they 
represent only the proportions of the specimens. The 
actual size required for any specific application de- 
pends on the intensity of the loading, g. The average 
thickness, ?, of the column (sheet and stringers) is 
found by dividing g by the allowable stress f. The 
actual dimensions are then chosen so that the proper 
average thickness is obtained. This amounts to 
contracting or expanding all the dimensions of the test 
specimen (sheet thickness, stringer thickness, stringer 
spacing, rivet diameter and spacing, etc.) by the ratio 
of the new effective column length, Lo, to the effective 
column length of the test specimen. The overall width 
of the column, however, is held constant at the value 
used in determining g. (This is the important difference 
between the simple column and the “‘wide’”’ column.) 


It is apparent that low loading intensities (low 
values of g) may require actual dimensions that are too 
small to be used in actual design. These limitations 
must be provided for in using the optimum stress 
curves for weight analyses. A method of doing this 
will be discussed later. 


It is important to note that, when a material of 
higher strength is used, the equivalent thickness of a 
column will decrease, since the higher allowable stress 
permits the use of less cross-sectional area. This de- 
crease in the amount of material tends to lower the 
allowable stress that can be developed, thereby par- 
tially offsetting the increase in material properties. 
As a result, the effect of material properties on weight 
isnot linear. This is automatically provided for when 
the structural index is used. 


Up to this point the column length, Z, has been 
treated as a variable. For a given value of g the maxi- 
mum possible allowable stress will be obtained by using 
the lowest possible value of L (see Fig. 5). This would 
correspond to a ‘‘continuous”’ type of support, such as 
that provided for a heavy cap strip by webs at right 
angles. Under such conditions the upper limit for the 
allowable stress is no longer determined by buckling 
but will depend on such things as allowable permanent 
set, buckling of other portions of the structure, etc. 
The compressive yield stress (at 0.002 permanent set) 
is usually regarded as a suitable value. 

When the bulk of the flange material (resisting axial 
loads due to bending) is spread out over the width of 
the box structure, the column length, Zo, is determined 
by the rib spacing and the fixity coefficient. A de- 
crease in Lo, therefore, corresponds to an increase in 
the number of ribs. It is obvious that for a given design 
there must exist some optimum combination of column 
design and rib spacing which will result in minimum 
weight. The allowable stress alone is therefore not a 
reliable measure of the strength/weight efficiency of 
the structure. It could be made to serve as such by 
including in the assumed cross-sectional area an addi- 
tional amount representing the average area of the ribs 
or equivalent stabilizing material. 

The true conditions for minimum weight evidently 
require a further “‘optimizing’’ process, in which the 
best relation between L and B is determined, for a given 
value of P/B. If such a relation is determined, the 
optimum column length (or rib spacing) may become 
a function of B, as well as P/B, since the weight of a 
rib will depend on the distance across which it must pro- 
vide support. The depth of the structure may also 
have to be considered, since it is obvious that a rib 
will be lighter in an extremely thin wing than in a thick 
one, other things being equal. 

From a practical point of view it will probably be 
found unnecessary to treat the weight of ribs (or equiva- 
lent stabilizing structure) with such refinement. For 
most structures the weight of this portion is relatively 
small in comparison with the weight of the directly 
loaded material. It may also be found that the true 
optimum structure involves impractical conditions, 
such as a large number of light ribs at close spacing. 
In any case, the envelope curves of allowable stress 
versus structural index will provide a basis for estimat- 
ing the weight under any assumed condition. 


The assumption that the structural material is uni- 
formly distributed over the width of the column may 
not represent the best practical structure in special 
cases. (For example, a “‘two-spar’’ type of design may 
be found most efficient for folding wings or for structures 
with large cutouts.) It does, however, give a reasonable 
basis for weight estimation, -since it eliminates any 
question as to the amount of material required simply 
to maintain the external contour. Other types of con- 
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struction can later be investigated separately and com- 
pared with the distributed type. 

It should also be noted that the trend toward thinner 
wings and higher speeds points toward the type of 
structure in which the material is uniformly distributed 
around the cross section for maximum torsional stiffness 
and surface smoothness. 


(5) MINIMUM WEIGHT FOR PANELS (SIMPLY SUPPORTED 
EDGES) 


The term ‘“‘panel”’ is used here (as in reference 2) to 
define a special form of wide column in which the edges 
are restrained and in which an appreciable part of the 
buckling strength is derived from the transverse bend- 
ing stiffness (resistance of transverse strips to bending 
under normal pressure). The edge restraint, in its 
simplest form, compels the edges to remain straight; 
it may also provide some degree of edge fixity, which 
further increases the buckling stress. 

The essential difference between the panel and the 
wide column is the fact that the width of the panel 
plays an important part, if not the major part, in 
determining the buckling strength. It is no longer 
possible to treat the panel as a series of adjacent 
columns of unit width; it must be treated as a whole. 

In the elastic range the buckling stress for a flat 
plate acting as a panel is given by an equation of the 
following type: 


f = kyE(t/b)? (17) 


where ¢ = the thickness of the plate and 6 = the panel 
width. 

The factor k, depends on the proportions of the panel 
(ratio of length to width) and on the degree of edge 
restraint. For panels that are considerably longer than 
they are wide, the factor k, becomes a constant. The 
physical significance of this is that the length (distance 
between ribs) no longer has any effect on the buckling 
stress. 

To determine the proper index to use for panels, it 
will be assumed that k, is a constant. Timoshenko’s 
effective modulus (reference 5, page 442), VEE, (with 
E, replaced by £,), will be used to provide for the in- 
elastic range. The buckling formula may then be 
written: 





f = ki(t/b)?WEE, (18) 


For a given load, P (per panel), and allowable stress, 
f, the required thickness is 


t = A/b = P/fb 
Substituting this in Eq. (18) gives 
f = ki(P*/Pb) VEE, 
from which 
f = ky(EE,)'"(P/b?)” (19) 


or 
P/b* = f'"/(EE,) “ky” (20) 


This shows that the structural index is similar to 
that for the simple column (P/L*) except that length 
is replaced by width. This result agrees with the assump- 
tion that the buckling stress is governed by the width 
and that length has no effect. Even in the case of 
relatively low ratios of length to width, where &, is no 
longer independent of length, the index P/b? may be 
used. It is then necessary to consider the length-width 
ratio (L/b) as a design constant, and the constant 
k, will have a different value for each value of L/B. 

In the equations for the flat sheet there is no shape 
parameter, such as is found in the previous equations 
for columns. 
sheet, the designer has no control over the disposition of 
material. Equations (19) and (20) therefore require no 
further modification to provide for other types of failure: 
they define the conditions for initial buckling. If 
initial buckling is taken as the criterion for failure 
Eqs. (19) and (20) may be used for design by con 
sidering P as the ultimate load. 

If b is taken as the width of the box beam in a con 
ventional single-cell wing structure, the corresponding 
value of the allowable stress will be entirely too low 
for any normal structural alloy. To obtain a suitable 
value, the width B must be reduced by using a number 
of ‘‘spanwise’’ members that divide the panel into nar- 
rower strips. When the load, P, is large in comparisot 
with the total width of the structure, it may be found 
practicable to use plain sheet material in this manner. 

When the total width of the box structure is sub- 
divided into strips by means of ‘‘shear’’ webs or similar 
stiffening members, it becomes necessary to distinguish 
between the total width of the structure and the actual 
width of each panel. For this reason the term } was 
used for the panel itself. The symbol B may be con- 
sidered to refer to the total width, as in Fig. 6. 


This simply means that, for a solid 














Fic. 6. Terminology for panel construction. 


If P is now considered to be the total load acting 
over the entire width of the structure, the load per inch 
of width is denoted as before by g = P/B. 

The load on each panel would be equal to gb. Sub- 
stituting this value for P in Eqs. (19) and (20) gives the 
equations 


f = ki'(EE,)'“(q/b)” 
q/b = f"/(EE,) "ky" 


(19a) 
(20a) 


Replacing E, by rE: 
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(19b) 
(20b) 


f “ ke’? E'*7'*(q/b)” 
q/b = f'0/Er ky 


For the “long” plate (L/B of unity or greater), the 
constant ky is equal to 3.6 (see Fig. 1-6 of ANC-5), 
and the equations become 


f = 1.532E'7'"(g/b)*”* (19c) 
q/b = 0.527(f/E'“r”) (20c) 


These equations should be compared with Eqs. 
(15b) and (16b). The physical difference is that the 
former equations apply to the plate used as a column 
with free edges, supported by ribs or bulkheads at 
spacing L. Eqs. (19c) and (20c) apply to the plate 
used as a panel, with no ribs or bulkheads, the edges 
being simply supported by webs (or similar structure) 
at spacing d. 

For a given value of f, in the elastic range, Eq. (20c) 
gives a value of g/b twice that given by Eq. (16b) for 
the column. This means that, for a given value of 
g aud to develop the same buckling stress, a plate may 
have “‘spanwise’”’ stiffeners spaced twice as far apart as 
the spacing required if “‘chordwise’’ ribs are used. Con- 
versely, for the same spacing of stiffening elements, a 
higher elastic buckling stress will be developed by the 
“spanwise’”’ stiffened plate, as shown by comparing 
Eqs. (19c) and (15b). (The ratio is equal to 2”* or 
1.58.) 

It should also be noted that the reduction of the 
tangent modulus with increasing stress has a less harm- 
ful effect in the case of the panel than in that of the 
column, since r appears to the one-half power in Eq. 
(16b) but only to the one-quarter power in Eq. (20c). 

It is obvious from these observations that flat plates 
should be stiffened in the direction of compressive 
loading. When so used, they represent an efficient 
type of structure in the regions of high loading intensity 
(high values of the structural index). 

To obtain a comparison of the flat plate with the 
stiffened panel, Fig. 7 has been prepared. Since the 
buckling stress in the elastic range varies as the two- 
thirds power of the structural index, the curves are 
plotted accordingly. 

It can be seen that the sheet-panel type of construc- 
tion begins to compete with the stiffened-column type 
when g/L reaches values in the neighborhood of 2,500 
lbs. per sq.in., provided that 75 ST or equivalent ma- 
terial is used. 

To show how various materials compare on this basis, 
the panel curves for some aluminum, magnesium, and 
steel alloys have been prepared on Fig. 8. The allow- 
able stress has been reduced to an equivalent weight 
basis, using aluminum alloy as a base. 

In the elastic range magnesium alloy is superior, 
but because of the falling off of the tangent modulus, 
the present commercially available alloys cannot com- 
pete with aluminum alloy in the range of high loading 
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Fic. 7. Comparison of wide-column and panel construction. 
intensity. It can also be seen that steel plate could be 


used efficiently only at extremely high values of the 
structural loading index and that high compressive 
vield stresses would be required. 

The equations for the flat plate, if properly modified, 
may be considered as generally applicable to other types 
of panel construction in which the resistance to buc- 
kling depends on both longitudinal and transverse stiff- 
ness. It can be seen that the ‘sandwich’ type of 
construction may be thought of as a flat plate in which 
the material is so arranged as to give an increased value 
of moment of inertia for the same cross-sectional area 
(or equivalent weight). This is comparable to the use 
of stringers on the wide column, except that it is neces- 
sary to maintain a high value of moment of inertia in 
all cross sections (cut in any direction). The shape 
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Fic. 8. Comparison of various materials used as panels (simply 
supported edges). 
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parameter p/?, or its equivalent, then becomes a design 
variable and the objective is to obtain high values of 
the shape parameter without introducing other types of 
failure, such as local buckling of the surface sheets or 
yielding of the core material. 

As in the case of the wide column, the optimum over- 
all design for panels must include the weight of the sta- 
bilizing members. Thus, if it were necessary to use 
closely spaced spanwise stiffening webs in order to de- 
velop a high allowable stress, the overall efficiency 
might be low. The depth of the box or shell structure 
will affect the weight of the stabilizing members, up to 
a certain point at least. 

The flat sheet stiffened in both directions by a grid 
arrangement of members may also be regarded as a 
special form of panel. The analysis of such a structure 
is taken up in detail in reference 2. 

Since qg is the load per unit of width, g/f gives the 
cross-sectional area per inch of width, or the average 
thickness 7. Dividing the structural index g/L» by the 
corresponding values for f gives values of 7/Lo. 


t/Lo = (q/Lo)/f (21) 


The values of 7/Ly have been plotted on Fig. 9 for 
the materials and types of structure previously in- 
vestigated. For weight-comparison purposes the 
values for steel and magnesium alloy have been reduced 
to the thickness of aluminum-alloy sheet, which would 
have the same weight. Some of the uses of these curves 
are shown later. 


(6) OpTIMUM OVERALL DESIGN 


For weight-analysis purposes, as well as for design 
work, it is desirable to know the optimum type of 
structure for any given case. Assuming that the total 


width and depth of a box beam have been determined 
and that the axial load intensity due to bending is 
known, it is possible to arrange the structural material 
in various ways, such as the following: 

(a) Heavy cap strips at corners, with light sheet 
and stiffening members forming remainder of box. 
(The cap strips act as simple ‘‘columns,”’ at zero column 
length. The auxiliary material must be included in 
the weight analysis.) 

(b) Wide column construction, with sheet, stringers, 
and full-depth ribs. (Most widely used at present for 
large airplanes.) 

(c) Panel construction with plate or sandwich 
covering, supported by spanwise full-depth webs or 
equivalent. 

(d) Wide column construction, with partial-depth 
ribs or frames forming a gridwork. 

Even though optimum arrangement of material for 
any one of the above types of construction may be 
determined in some manner, these types are not all 
comparable on the same basis. The final step in the 
“optimizing”’ process is therefore to determine how the 
minimum weights for each type of construction com- 
pare, as a function of certain design parameters. 

It can be seen that in every case the structural index 
contains the load intensity factor, g (=P/B). The 
dimensions that determine the overall efficiency of a 
structure, for a given load intensity, are the width of 
the box beam, B, and its depth, 4. The depth may be 
considered as a secondary variable, since it influences 
only the stabilizing structure. The proper overall 
parameter to be used for general design and weight 
analysis, therefore, appears to be g/B. 

To account for the effects of the depth of the struc- 
ture, the nondimensional ratio h/B may be used. For 
a given value of this parameter a series of structural- 
efficiency curves might be plotted against the index, 
q/B. To account for the variation in specific weight 
of different materials, the structural-efficiency curves 
may be plotted in terms of the weight per unit of load 
per unit of length. (For a given material this is equal 
to the specific weight of the material divided by the 
average allowable stress.) For certain types of con- 
struction, the effects of combined loading conditions 
may have to be considered. This can usually be accom- 
plished by means of interaction curves, which indicate 
how much the allowable compressive stress must be 
reduced to account for other types of simultaneous load- 
ing (such as normal pressure, shear due to torsion, etc.). 
It is also recognized that the requirements for high 
torsional stiffness may in some cases override those for 
adequate strength. This is considered to be an entirely 
separate problem, since the optimum type of structure 
is not necessarily the same in both cases. 

The collection of test data, weight data, and other 
information from which to construct the types of 
curves mentioned above represents the ultimate solu- 
tion to the optimum weight problem. Pending the 
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completion of this long-range project, however, some 
useful results may be derived from the knowledge al- 
ready at hand. 


(7) THE STRAIGHT-LINE AREA APPROXIMATION 


On Figs. 4 and 9 it can be seen that some of the 
curves can be approximated quite well by straight lines 
that do not pass through the origin. Fig. 10 shows such 
an approximation for the two cases of the Y-type stiff- 
eners of 24 ST and 75 ST aluminum alloy. (The basic 
curves are obtained from Fig. 5 by dividing abscissae 
by ordinates.) The equation for the straight lines has 
the form 


B/Ly = Ca + [(q/L0)/Fa] (22) 


where C, represents the intercept value and F, is a 
fictitious allowable stress. 
If Eq. (22) is multiplied by Lo, it becomes: 


P= C,Lo + (q/Fa) (23) 


Since the cross-sectional area is equal to ¢B and 
since g = P/B, Eq. (23) may also be written 


A = C,Bly + (P/F,) (24) 


The last term of this equation is the usual expression 
for the cross-sectional area required to transmit the load, 
P, at an allowable stress, F,. The other term repre- 
sents an additional cross-sectional area that depends 
only on the column width, B, and on the effective 
column length, 2». The values of C, and F, are deter- 
mined by the optimum design envelope for a given 
material and type of construction. For the materials 
shown in Fig. 10 the values are given in the following 
table: 


Material Co F,, Lbs. per Sq.In. 
24 ST 0.0012 49,000 
75 ST 0.0015 74,000 


Fig. 9 shows that it would be difficult to fit the entire 
curve for conventional 24 ST construction by means of a 
single straight line. Up toa value of g/Lo of about 700 
lbs. per sq.in., a line having the constants 0.0015 and 
51,000 Ibs. per sq.in. fits quite well. Beyond that 
point the weight would be overestimated. A line hav- 
ing the constants 0.0050 and 66,700 Ibs. per sq.in. would 
fit the upper part of the curve. 

The physical significance of this approximate method 
is interesting. In Eq. (24) the term, F,, may be 
thought of as a constant allowable compressive stress, 
provided that the other term is also used. This ‘‘allow- 
able stress” is rationally determined by the optimum 
design envelope for the material and type of construc- 
tion involved and not by an arbitrarily defined compres- 
sive yield stress. 

The constant C, also depends on the material and 
type of construction but has approximately the same 
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Fic. 10. 


value for the three cases shown. It represents a mini- 
mum value of //Lo, which is very much like establish- 
ing a minimum practical value for the slenderness ratio 
of columns. The portion of the cross-sectional area 
represented by C,BLo might be thought of as the 
area required for ‘‘self-support’’ against buckling. 
This portion of the area is needed regardless of 
load. 

The situation is analogous to that in aerodynamics, 
in which drag is composed of two main parts, parasite 
and induced. The fact that the true optimum curve is 
based on using certain proportions of the cross section 
is comparable to the fact that the induced drag equa- 
tion is based on an optimum shape of the wing (ellip- 
tical). 

Turning back to the simple column, the curves of 
Fig. 4 can likewise be expressed by equations of the 
type 

A/Ly? = C,. + [(P/L0?/F:] (25) 
A = C,Ly*? + (P/F,) (26) 


Here the term, C,L”, represents the portion of cross- 
sectional area needed for “‘self-support’’ against buc- 
kling. For the thin-walled round tube the constant C, 
will be found to correspond to some maximum values 
of slenderness ratio (L/p))o and diameter/thickness 
ratio (D/t)o. The following equation defines the rela- 
tionship between these quantities, when the load, P 
is zero: 


be = S2/(L/p)o?(D, t)o (27) 


A value of C, = 0.000025 is obtained, for example, 
by letting (L/p)o = 120 and (D/t)) = 70. 

The values for the constants to be used in Eq. (26), 
as obtained from Fig. 4 for three aluminum alloys, are 
as follows: 
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Fic. 11. Efficiency of shear-carrying members. 
Material C. F, Lbs. per Sq.In. —_ be overestimated by thisapproximation. For example, 
24 ST 0.000025 42,000 the constant C, = 0.003 means that a shear web 
14 ST 0.000025 57,500 of 10 in. depth would have a minimum average thick- 
75 ST 0.000040 76,000 


It should be noted that the curves used for columns 
were determined theoretically and that the above con- 
stants are therefore subject to modification as a result 
of tests. 


(S) SHEAR-CARRYING MEMBERS 


In some of his earliest papers on shear beams Her- 
bert Wagner introduced the parameter W/S/h, where S 
is the shear load and / the depth of the member over 
which the shear is transmitted. Fig. 11 shows how the 
average allowable shear stress (ultimate shear load 
divided by average cross-sectional area) varies with the 
structural index. Most of the points are obtained from 
N.A.C.A. tests. A solid line has been drawn to repre- 
sent an envelope for webs with beaded holes. Such a 
curve may be converted to a curve of A/h? plotted 
against S/h? by dividing abscissas by ordinates. The 
resulting curve is shown in Fig. 12 and can be approxi- 
mated by equations of the type: 


A/h? = [(g,/h)/F;] (28) 
A = C,h? + (S/F,) (29) 


For the straight line shown in Fig. 12, C, 0.003 
and F, = 21,500 lbs. per sq.in. It can be seen that the 
minimum weight of lightly loaded shear members might 


ness of 0.030 in. for webs with beaded holes. 

The data on optimum design of shear members are 
obviously far from complete, and it would be unwise to 
draw general conclusions from Fig. 11, except to note 
the marked absence of test data for certain types of 
construction, especially in the higher loading ranges. 
If all reliable test data on shear-carrying members 
could be collected and plotted as in Fig. 11, it might 
be possible to establish more definite envelopes. It is 
hoped that this can be done in the future. 
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Fic, 12. Unit weight of shear webs with beaded holes. 
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Fic. 13. Modulus of rupture for round tubes and shells in pure 
bending (24 ST material). 


(9) MInIMuM WEIGHT FOR CYLINDRICAL SHELLS IN 
BENDING 


The unstiffened shell of circular cross section and 
uniform wall thickness is used as a basis for deriving 
the structural index and obtaining typical curves of 
optimum bending stress. The allowable stress for pure 
bending is usually expressed in the form of a modulus of 
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to express the modulus of rupture by means of a formula. 
It is assumed only that the modulus of rupture depends 
on the value of D/t for a given material. An arbitrary 
curve is drawn to represent this relationship, and the 
analysis is carried out on the basis that such curves are 
available for various materials. 

The modulus of rupture is assumed to depend on 
D/t in some manner such as that indicated in Fig. 13. 

This curve has a shape roughly approximating an 
equation of the type f, = k VEE, /(D/t), but because 
of the effects previously noted, the value of k will prob- 
ably vary with D/t. 

For a given value of applied bending moment, the 
thickness required is given by the equation: 





t = 4M/nD*f (30) 
Dividing both sides by D, 
t/D = 4M/nrD*f (31) 


This equation may be used to determine the required 
value of D/t if the allowable stress f is known. The 
structural index is found to be the quantity M/D*. 


To obtain a curve of allowable stress as a function of 
the structural index 1//D*, Eq. (31) may be written 


M/D* = (x/4)[f/(D/d] (32) 


For any assumed value of f, the corresponding value 
of D/t may be read from the curve of Fig. 13, and the 


0 rupture, which is defined as 
fy = M/(/c) 
where JJ = ultimate bending moment and J/c = sec- 


tion modulus. 
In the region of high stresses the actual stress dis- 
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Fic. 14. Optimum modulus of rupture for round tubes and shells in pure bending (24 ST material). 
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bending (24 ST material). 


index can be calculated and the modulus of rupture ob- 
tained from a curve such as that shown in Fig. 14a. 
The required thickness and cross-sectional area are 
obtained from the formulas: 


t = 4M/nrD*f (33) 
A = rDt = 4M/Df (34) 


From Fig. l4a it might appear that the minimum 
weight for a given bending moment would be obtained 
by using the smallest possible diameter, since this would 
give the highest stress. Eq. (84), however, shows that 
the required cross-sectional area varies not only in- 
versely as the stress but also inversely as the diameter. 
Hence, for a given bending moment, the required cross- 
sectional area will be a minimum when the product Df 
is a maximum. 

To obtain the minimum value for the cross-sectional 
area, a Straight line may be drawn on Fig. 14b, tangent 
to the curve. This line determines the minimum value 
of the quantity 1/'”/(Df), which for the curve in ques- 
tion is equal to 1.52 K 10~*. 

Eq. (34) is divided by M™, giving 


A/M™” = 4(M'“/Df) (35) 
If the minimum value for the term (1/'*/Df) is sub- 
stituted in the equation, we obtain 
(A/M*")min, = 4(M'"/Df) min. = Rs (36) 
For the case illustrated, 
(A/M") nin. = ks = 6.08 & 10-4 


Thus, for the optimum design of an unstiffened 
cylinder in bending, the cross-sectional area will vary 
according to the following equation: 


Aopt. = ks” (37) 


where &; is a constant that depends only on the materia] 
involved. 

Fig. 14b shows also that there is an optimum stress 
for minimum cross-sectional area (50,000 Ibs. per sq. in, 
in this case), from which it follows that there is a single 
value.for the optimum D/t (=90). 

The optimum diameter for a given bending moment 
may be found from the optimum value of M/‘”*/D, 
Let 

ke = (M“/D)ont 
then 
Don, = M/k (38) 

For the case illustrated in Fig. 14b, kj = 7.6. A 
different value would be found for other materials or 
for a different assumed curve of modulus of rupture 
versus D/t. 

Fig. 14b may be converted to a curve of A/M“ 
versus the structural index, ./'”,D, by using Eq. (35). 
Such a curve is shown in Fig. 15. This curve shows the 
optimum value for the structural index and the mini- 
mum value for the term A/M*™*. It also shows how the 
cross-sectional area varies with the structural index for 
a given value of the bending moment. In this case the 
curve is relatively flat in the region of minimum cross- 
sectional area, indicating that the area is not sensitive 
to changes in diameter. For example, the limits of 
M'“/D, within which the area will not exceed the 
optimum by more than 5 per cent, are found from Fig. 
15 to be approximately between 5.2 and 10.8. From 
Fig. 14b it can be seen that this corresponds to a range 
of stress (modulus of rupture) approximately between 
33,000 and 67,000 Ibs. per sq.in. It is useful to know 
that a variation in diameter of approximately 30 per 
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Fic. 16. Design type of chart for unstiffened circular cylinders 
in bending (24 ST material). 
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cent on either side of the optimum will not increase the 
cross-sectional area (weight) by more than 5 per cent. 

For design purposes it is convenient to show the re- 
sults in terms of actual bending moment, diameter, 
thickness, and cross-sectional area. Fig. 16 shows one 
way of doing this. From such a chart the designer can 
determine the diameter at which the weight will be a 
minimum. If this cannot be used, he can select any 
desired diameter and note the corresponding increase in 
required cross-sectional area. 

A chart such as that shown in Fig. 16 may be pre- 
pared for any material and for any cross-sectional shape, 
if the variation of bending modulus of rupture with D/t 
(or equivalent ratio) is known. 

It can be seen from Fig. 16 that the diameter for 
minimum weight is relatively small (for a bending 
moment of 10,000,000 in.Ibs. it is in the region of 30 
in.). This, of course, shows why the unstiffened cir- 
cular shell is not efficient for such structures as fu- 
selages, for which ther equired diameter exceeds the 
optimum by a large amount. It should be more suit- 
able for long slender members, such as booms. 

Since the bending moment for a concentrated load 
varies directly as the length of transmission, Eq. (38) 
gives the variation in optimum diameter with length. 
The general shape involved is shown in Fig. 17a. The 
shell thickness is shown exaggerated to indicate that it 
should vary at the same rate as the diameter (for 
constant D/t). 

If the load being transmitted is uniformly distrib- 
uted along the path of transmission, the bending mo- 
ment will vary as L’, and the optimum diameter will 
therefore vary as L’”*. The general shape involved is 
shown in Fig. 17b. 

(10) MinimuM WEIGHT FOR STIFFENED SHELLS IN 
BENDING 


As previously noted, for a typical large fuselage the 
diameter is far beyond the optimum for an unstiffened 
shell. Since the design bending moment and the diam- 
eter are usually specified, the only thing that can be 
done to improve the efficiency is to increase the effect- 
ive thickness of the shell. This is done by the same 
general methods used for wide columns atid panels— 
that is, by utilizing some of the shell material in the 
form of stringers, frames, or both. 

In the usual type of construction the shell is stiffened 
by longitudinal stringers that are supported by frames 
or bulkheads. For shells of large diameter the strength- 
ening effect of the curvature of the sheet can be neg- 
lected without introducing appreciable errors. The 
strength of the shell is then governed by one of two 
types of failures: (a) wide column failure of sheet and 
stringers between frames; and (b) general instability 
of the shell, including frames. 

The most efficient construction will represent a com- 
promise between these two types of failure. 





THICKNESS GREATLY EXAGGERATED 





(a) CONCENTRATED LOAD AT ONE END 
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(b) UNIFORMLY DISTRIBUTED LOAD 


Fic. 17. Optimum proportions for cylindrical shell in bending. 











Conversion of shell structure into equivalent axially 
loaded structure. 


Fic. 18. 


By considering the frame spacing to be a design con- 
stant, it is possible to apply to the shell the results 
already obtained for the wide column. The axial 
compressive load per inch of circumference will be a 
maximum at only one point, but it will remain sub- 
stantially constant over a considerable fraction of the 
circumference on either side of this point. (For ex- 
ample, the value at 20° on either side of the centerline 
will have dropped to approximately 94 per cent of the 
maximum. ) 

For a circular shell of uniform cross section, the 
beam theory (f = Mc/I) gives the following equation 
for the maximum value of q: 

gq = (4/m)(M/D?*) (39) 

The optimum stress for wide columns has been deter- 
mined (Fig. 5) as a function of the structural index g/L. 
Eq. (39) is therefore divided by L, giving 

q/L = (4/x)(M/D°*L) (40) 

This is the equation by which the uniform circular 
shell in bending can be converted into the equivalent 
wide column under axial compression. To give a more 
direct physical picture, the following relationships can 
be used: 

Equivalent axial load: 


Peg, = M/D (41) 


Equivalent width: 
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Beg, = (r/4)D (42) 


Chis is illustrated in Fig. 18. 

From these relationships the designer may determine, 
on Fig. 5, the optimum stress for any combination of 
bending moment and diameter, considering L to be the 
spacing of the frames or bulkheads. 

This does not take into account any increase in allow- 
able stress caused by the curvature of the shell; nor 
does it allow for the effects of shear lag, which occurs 
when a large part of the bending moment is caused by 
a high shear load. For typical fuselage construction, 
however, these effects are usually negligible; in any 
case, the use of the wide column data in the above 
manner is likely to be on the safe side. 

In Section 9 it was shown that for the unstiffened 
circular cylinder there is an optimum diameter corre- 
sponding to any given value of bending moment. At this 
diameter the optimum stress has a constant value. 
The same methods will now be applied to the stiffened 
evlinder. 

It is assumed that a given bending moment, J/, is to 
be carried by a stiffened shell of circular cross section 
and that the diameter D may be varied so as to obtain 
the lightest possible construction. It is further as 
sumed that the same construction is to be used over the 
entire circumference of the cross section. (This agrees 
fairly well with actual construction, especially where 
the design bending moment must be applied about more 
than one axis.) 

The average thickness of the shell plating and stiff- 
eners (excluding frames) is given by 


i = q/f (48) 
where f is the allowable compressive stress. The total 
cross-sectional area is then 


A = rDi = xDq/f 


Substituting for g from Eq. (39), 
A = 4(M/Df) 


This is, of course, identical with Eq. (34) derived for 
the unstiffened cylinder. 

Eq. (35), derived for the unstiffened cylinder, can 
also be used here: 


A/M” = 4(M'“/Df) (35) 
Let 
L = kD (44) 


Then Eq. (40) may be written as 
g/L = (4/2k:)(M/D*) (45) 
(q/L)'* = (4/mk:)'“(M'“/D) (46) 


The term A/M™ in Eq. (35) has a minimum value 
when the quantity (1/'"/D)/f is a minimum. The 
conditions under which this occurs for a certain type of 
construction can be found by replotting the appropriate 
envelope curve from Fig. 5 against (q¢ L)‘“, as shown in 
Fig. 19. The method used for the unstiffened shell 
can now be used to find the optimum stress by drawing 
a line tangent to the curve, as shown. For the partic- 
ular envelope curve used in Fig. 19 (24 ST conven- 
tional construction), the optimum stress turns out to 
be approximately 25,000 Ibs. per sq.in. 

The value for (q ) Bid is found to be about 3.9. 
Substituting this term in Eq. (46), 


(AL'"/D)opt. = ke(wk:/4)' (47) 


where kg = (q Loe. obtained from the tangent to the 

optimum column stress curve, plotted against (q L)* 

and k; = the frame or bulkhead spacing ratio = L/D. 
Eq. (47) may be written 


Dopt. = (4/m)'*(M'"/ kk”) (47a) 


It can be seen that the optimum diameter for a given 
bending moment depends not enly on the material used 
(governed by ks) but also on the frame spacing ratio, 
kz. 

As the frame-spacing ratio is increased, the optimum 
diameter tends to decrease for a given value of the 
bending moment. If the frame-spacing ratio is held 
constant, the optimum diameter will vary as M™, 
which agrees with the result already found for the un- 
stiffened shell. 

By using Eqs. (46) and (35), together with an en- 
velope curve for a particular material, curves of A/M"" 
may be plotted against the index M/'”“/D for various 
values of the frame-spacing ratio k;. Such curves are 
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shown in Fig. 20 and may be compared with Fig. 15 for 
the unstiffened shell. 

For design purposes Fig. 20 can be converted into 
terms of cross-sectional area, diameter, and bending 
moment for different values of the frame-spacing ratio. 
It is also possible to plot the curves for constant values 
of the frame spacing (effective column length). Both 
types of curves are shown on Fig. 21, which has been 
prepared for only one value of the bending moment 
‘1,000,000 in.Ibs.). It is interesting to note that, for a 
constant frame spacing and for the assumed value of 
bending moment, increases in diameter beyond 100 in. 
approx.) have practically no effect on the cross-sec- 
tional area required. 


(11) CONCLUSIONS AND RECOMMENDATIONS 


[he structural index provides a sound and con- 
venient basis for: 

(a) Presentation of test data for various types of 
construction and different materials (this will assist 
greatly in planning research programs). 

(b) Comparison of various designs for a given 
material, resulting in determination of optimum propor- 
tions for minimum weight. 

(c) Comparison of various materials on an optimum 
basis, thereby allowing for the fact that the optimum 
type of design varies with the material and the loading 
conditions. 

(d) Determination of a rational ‘allowable stress’’ 
for structures subject to buckling (by means of the 
straight-line approximation for area curves). 

The straight-line approximation provides a con- 
venient basis for weight estimation, since it permits 
the required cross-sectional area to be divided into two 





OASH LINES INDICATE AVERAGE ALLOWABLE 
STRESS BELOW I0 ks 


\ L 
\ OPTimuM FoR { 


25,000 ps 








} | | 
C/ | } | | | | 
| 
SS SS SS SS SSS SS SS SS SSS SSS SSS SSS SSS SSS SESS SSS SSE SES EEES 
2 3 4 5 


u's % 
bw: (psid 


on 


Cross-sectional area of stiffened circular cylinders in 
bending (24 ST material). 


Fic. 20. 


FOR MINIMUM WEIGHT 149 


BENDING MOMENT mst0°in Ib 


24 ST UNSTIFFENED SHELL 
(FROM FIG. 13) 





£ 
Cc 
a” 
z 
5 DASH LINES INDICATE AVERAGE ALLOWABLE 
8 STRESS BELOW 10 ksi 
yy 
” 
° 
og 
oO 
cs 
ra) 
< 
w 
« 
a 
= 25000 i 
fopt ~ 
) 20 40 60 80 100 «= 120—S ss 140—S—s«*1600 
DIAMETER OF SHELL, in. 
Fic. 21. Design type of chart for stiffened circular cylinders in 


bending (24 ST material). 


imaginary parts, one of which depends only on the 
loading and the material properties, while the other 
depends on the dimensions of the structure. 

It is recommended that more emphasis be placed on 
the ‘‘design’’ approach, through the use of the struc- 
tural index, instead of on the “‘stress analysis’’ approach. 
The latter is best adapted to checking the adequacy of 
a structure already designed, while the structural index 
method permits the selection of the best structure for a 
given force transmission requirement. 
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The Oscillating Wedge in a Supersonic 
Stream 


G. F. CARRIER* 


Brown 


ABSTRACT 


When an obstacle in a fluid stream is allowed to execute oscil- 
lations, the pressure distribution on the obstacle will be modified. 
Such a modification frequently implies a dynamic stress resultant 
that leads to the flutter phenomenon. Theoretical investigations 
of this situation for those cases where the stream is supersonic 
and the obstacle is slender have been pursued by several authors. 
A rather complete bibliography and a detailed summary of these 
It is to be noted that 
each of these investigations relies on a linearized theory (i.e., a 


investigations can be found in reference 1. 


perturbation of a given steady flow) and neglects the presence of 
the shock wave. If one adopts the point of view that the oscil- 
lating obstacle produces acoustic waves that must reflect from 
the shock, it is evident that the dynamic pressure distribution 
may be considerably modified by these wave reflections. It 
therefore seems of importance to determine the modifications that 
arise in the theory of an oscillating airfoil in a supersonic stream 
when one takes the presence of the shock into account. 

The present paper is a presentation of such a modified theory 
for an oscillating wedge. Specifically, we shall consider a wedge 
of finite chord and infinite span which executes small oscillations 
about its leading edge. The perturbation wave functiors are de- 
termined in a conventional product-series expansion, and re- 
currence formulas for the associated coefficients are readily ob- 
tained. The results are compared to those for the no-shock (i.e., 
indefinitely thin airfoil) situation. 


THE DIFFERENTIAL EQUATIONS AND BOUNDARY 
CONDITIONS OF THE PROBLEM 


iy Fic. 1, we have indicated the wedge and the field 
of flow which are to be considered here. In the re- 
gion R (above the wedge and behind the shock), the 
steady flow is uniform and supersonic. When the 
wedge is allowed to oscillate, the flow field is perturbed, 
and our problem is to find the “‘wave functions’ asso- 
ciated with this motion. For waves of small ampli- 
tude, it has been shown? that these wave functions obey 
the equations 


Ay — (d*o/dt?) = 0 (1) 
dE/dt = 0 (2) 
dS/dt = 0 (3) 
q = —dy/dt (4) 


where the particle velocity is 
V = U, + a (grad ¢ + curl E) (5) 
the pressure and density are given by 
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pb = p2(1 + q), p=p(l+q-—S) (5a) 


where a2 is the downstream steady-state sonic velocity, 
M = |U,| ‘do, d/dt = (0/dt) + (.\10/dx), and tis the time 
multiplied by a2. 

Eqs. (5) and (5a) are, in reality, definitions of the 
quantities y, E, and S, which are introduced only for 
convenience in manipulation. However, one may 
associate a physical meaning with each: yg may be con- 
sidered as representing an acoustic wave; E defines a 
Finally, S is the 
entropy increase associated with these waves. Eqs 
(1) to (4) are, of course, the result of placing Eqs. (5 
and (5a) in the conventional equations defining the con- 
servation of mass, momentum, and energy. 

The boundary conditions at the shock have also been 
developed in reference 2 and may be writtent 


[grad ¢ + curl E]* = [1 — (p:/p2)] X 
[ian (ys + nf.) + Jj awy,| (6 


[— dy/dt|* = [1 — (pi/p2)]as(vi + ny) (7) 
[S]* = [1 — (1/pe2)las(ys + ny,) (8 


wave with nonvanishing vorticity. 


ti, k,j, are respectively directed normal to the shock, parallel 
to the wedge leading edge, along R:-i. 











Fic. 1. The flow field of an oscillating wedge problem. The 
region R is that region in which the solutions defined ,here are 
valid. 
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Here, the asterisk indicates that the bracketed quan- 
tities are to be evaluated at the shock locus; y and y, are 
the shock displacement and slope as defined in Fig. 1. 
The a, are defined by 


a, = (2 + 4m*)/3(1 — m?*) (9) 
a2 = (5 + m*)/6m (10) 

az = —m(5 + m*)/3(1 — m?) (11) 
ay = —(1 — m*)/3m (12) 


and m is the component of U2/a2, which is normal to 
the shock, and m is that component tangent to the 
shock. These shock boundary conditions are, of 
course, the conventional Rankine-Hugoniot conditions 
written to apply to a curved moving shock and linear- 
ized with regard to the small disturbances. Unlike the 
analysis in reference 2, we include here only two- 
dimensional disturbances, since the wedge motion does 
not imply nonvanishing velocity components in the z 
direction. This allows us to use the single component 
vector 
E=kE 
Che boundary condition at the wedge surface requires 
that the normal component of velocity ¢, — E, should 
be consistent with the prescribed motion of the wedge. 
If the motion of the wedge is an oscillation about its 
leading edge, then the normal component of the par- 
ticle velocity at its surface can be written in linearized 
form as (see reference 1, p. 29): 

(gy — Ex)yao = [x — ((M/c)]e™ (13) 
where c is the circular frequency of the motion and 
where the amplitude of the angular velocity of the wedge 
isunity in some unit system. 


THE WAVE FUNCTIONS 
It was shown in reference 2 that if we define 7? = x? — 
8’y?, tanh 6 = By/x, where 8? = M? — 1, then solutions 
of Eq. (1) may be written in the form 


y,(x, y, t) = [a, cosh (v8) + 
b, sinh (v0)] J,(kr)e!~O%/0)) (14) 


Solutions of Eqs. (2) and (3), which will be of use here, 
are (when £2? = 1 + \? — M?): 


E.(x, y, t) = c,J, (Rev) etl! — @/) — Oy/ MB) (15) 


and 
S.(x, y, t) = h, J, (Réy)ett!— G29 —-Os/Mey (16) 


Finally, it will be convenient to develop y in a series 
whose terms are: 
v, = [1 — (p1/p2)]~! d, J, (REy et —OM9/691 (17) 


With regard to the functions defined in Eqs. (14) to 
(17), one should note that the only terms that give a 
nonvanishing contribution to ¢,(x, 0, t) — E,(x, 0, #) 
are those with coefficient b, provided v > 0. In fact, 
if we write 


Cbs neers (18) 


$e 0208 08 8 4 


then the boundary condition (13) becomes (if k = ¢/8?): 


> 5,8,J,(kx)e~*™* = x? — (iMx/c) — (19) 
v=1 
However, the fundamental expansion theorem for 
Bessel functions states that 


e/lt-/Oln — Jo(n) + le” + (—t)~"|J,(n) (20) 
v=l1 


Thus, if we define ¢ = 7(1/ + 8), n = kx, the left side of 
Eq. (20) becomes e™”, and appropriate differentiations 
with regard to M yield the expansions of ne and 
ye". The useful result is: 

b, = k-*p-* [? + (-)~"] (21) 
The lowest order nonvanishing coefficient is, as one 
should expect, 1. 

It should be noted that, if we consider the laterally 
oscillating wedge, the only change in the analysis re- 
quired would be to change the 3,. 

If we now use the functions of Eqs. (14) to (18) to 
satisfy the boundary condition (6) to (8), we obtain an 
infinite nonhomogeneous set of equalities in the co- 
efficients a,, ...d,. After the proper collecting of 
coefficients, the recurrence formulas become as shown 
on top of page 152. 

With reasonable caution, a limiting process where 
6) —+ © produces the results applicable to the case 
where the shock degenerates to a Mach line. For this 
case, one obtains 


a,/b, = —1, c =d,=0 23) 


We can readily investigate the importance of taking 
the shock presence into account by computing the a, 
(for say, vy < 4) associated with given m, m (i.e., given 
wedge angle and upstream Mach Number). If these 
values (the —a,) differ appreciably from the b,, we 
should conclude that the presence of the shock cannot 
be disregarded in computing the pressure, —dg/dt. 
For the casest m? = 0.9, = 1; m? = 0.3, n? = 1; we 
obtain, respectively, 


a/b, = —1.009,  d2/b» = —1.00005 + 0.0687 iby /be 
and 
a,/b, = —1.182, d2/bo = 01.094 + 1.0057%b,/b- 


It is evident that the weak shock solution is not appre- 
ciably altered by the shock boundary condition omis- 
sion but that the solution associated with the stronger 
shock undergoes an appreciable change. 


+ Extremely weak and extremely strong shocks, respectively. 














—M sinh v6, M sinh vA 


21M cosh vy 
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V1 — m? 
cosh vo é pda = — ae es. (may, + a2) a, +1 
se: V1 = m' 
B sinh vO 0 — E tien — na) Cott 
—M cosh v4 0 nV 1 — m? — m?* as —d,41 
till aw 
—cosh vA —&£ a. —— (ma, + a) 
wt naa 
+ + sinh v6 0 ee (maz — n*ay) ios 
+ M cosh vA 0 - nV 1 — mas —d,_1 
— in % (22) 
—2iM cosh v6 U “s (1 — m*)ma, + nae a, 
2in 
° 2187/M — — Maz — (1 — m*)a, & 
M 
27 cosh v6 0 —2i(1 — m?*)az3 —d, 
sinh v6) —sinh v6) 22A\/ cosh vA 
+ B cosh v% 8 cosh vM% 0 


where tanh # = B/X. 








It is hoped that more detailed numerical results as 
to dynamic lift, etc., will be forthcoming soon. 


APPLICABILITY OF RESULTS 


It seems evident that if the motion of the leading 
edge of the obstacle were sufficiently violent, the shock 
might detach from the wedge, and, in fact, the flow 
might simulate the flow past a blunt obstacle. We can 
obtain an estimate of when this is apt to occur by the 
following consideration. Let us suppose that the wedge 
has a uniform lateral velocity normal to the free stream 
direction. Then, effectively, we have an oblique 
steady flow relative to the wedge, and the question of 
whether the shock will be detached is readily answered 
once the wedge angle, the free-stream state variables, 
and the wedge velocity are given. In general (i.e., if 
we are not operating near the critical condition where 
the shock is detached), the wedge velocity? required 
to produce this detachment is of order unity (say 0.4, 
for example). Now, if we translate these results into 
the oscillating wedge problem, we associate this hypo- 


t I.e., wedge velocity /sonic velocity. 





thetical steady wedge velocity with the maximum os 
cillating wedge velocity, which is bw where b is the wedge 
displacement amplitude and w the circular frequency of 
the motion. Thus, dw/a ~ 1, and for model tests ir 
a moderate sized tunnel, where 0 is limited to say !/ in 

the frequency must be of order 10° rad. per sec. From 
the point of view of the /inear stability problem, 3 is 
taken indefinitely small so that w < © isa suitable cri- 
terion for a nondetaching shock. For large displace 
ments such as occur in wing flutter, any correct theory 
becomes intrinsically nonlinear, and our results would 
Thus, for the early stages of 
For small 


not apply in any event. 
wing oscillation the theory appears valid. 
finite amplitudes the validity range (of b and w) can be 
estimated in the manner described above. 
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Base Pressures at Supersonic Velocities’ 


FREEMAN K. HILLT anp RALPH A. ALPHER? 
Applied Physics Laboratory, The Johns Hopkins University 


ABSTRACT 


Results are given of flight tests in which base pressures were 
directly measured in the range of Mach Numbers 0.8< M<2.5. 
The base pressure coefficient, p,/g, was found to be nearly con- 
stant for 1.2< 4/<1.6 and then to decrease in absolute magnitude 
with increasing Mach Number. The values of )/q appear to be 
significantly larger in absolute value than were obtained from 
wind-tunnel tests at JJ = 1.73 on a model geometrically identical 
to one of the free-flight test vehicles and from wind-tunnel 
measurements on other models. However, all the wind-tunnel 
measurements were at Reynolds Numbers considerably lower 
than those in the flight tests, the latter being of the order of 
2x 107 tol X 108 at J = 1.6. This suggests that at supersonic 
velocities there may be a significant scale effect on base pressure. 
The distribution of pressure over the base was found to be 
essentially constant. No dependence of base pressure on body 
length was found in free flight and only a slight dependence on 
angle of attack. The functional dependence of base-pressure 
coefficient on the parameters Mach Number, Reynolds Number 
body length, position on the base, and body angle of attack is 


given in graphical form. 


(I) INTRODUCTION 


| BASE DRAG of a body in supersonic flight is 
known to contribute a large percentage of the 
total drag, and this base drag is the integrated under- 
pressure acting on the base. Determination of the 
pressure on the base of missiles, projectiles, or super- 
sonic aircraft is therefore of practical importance. 
So far no satisfactory theory has been developed, 
nor have many data been published of an experimental 
nature from wind-tunnel studies. The fact that there 
have been few data from wind-tunnel work is due 
primarily to the difficulty of unambiguously measuring 
base pressures in a tunnel. It was therefore believed 
that base pressures measured in free flight would be of 
value in (a) establishing the behavior of the actual 
base pressure on a missile of simple geometry and (b) 
providing a basis for evaluating the dependence of base 
pressure, as measured in wind tunnels and in free 
flight, on Reynolds Number. 

This report gives results of a series of experiments 
measuring the pressures acting on the base of a missile 
in free flight at supersonic velocities.' These results 
are believed to be unique in that they represent values 
actually measured in flight under full-scale conditions 
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as a continuous function of Mach Number. The 
experiments were designed so that the smallest full-scale 
missile in this series could also be tested in a large 
supersonic wind tunnel (test section, 27.5 by 19 in.) 
for direct correlation of the flight results with tunnel 
measurements. 

Free-flight measurements are not subject to certain 
limitations associated with wind tunnels and with 
indirect ballistic determinations. The base pressure 
data from wind-tunnel models are subject to the general 
criticism that the mounting sting or support, as well 
as shock waves reflected from tunnel walls, may produce 
some unknown and unwanted effect. It has been 
demonstrated that a base sting has at least a strong 
local effect on base pressure. In ballistics, the deter- 
mination of base pressure, theoretical head drag, and 
skin friction drag are subtracted from measured total 
drag; the remainder is attributed to base drag. Thus, 
the indirectly obtained base drag includes the errors 
from two theoretical computations and from a total 
drag measurement. Moreover, the high rate of spin 
required for stability in ballistic tests may affect the 
boundary layer and, hence, also the base pressure. 

Free-flight measurements are direct and, in the 
experiments reported here, are obtained continuously 
through a wide range of Mach Numbers in a single 
flight, in contrast to wind tunnel and ballistic tests. 
It should also be noted that wind-tunnel operation is 
impracticable near Jf = 1, whereas the data from free 
flight are continuous through the transonic region. 
Furthermore, free-flight data are obtained under: 
actual atmospheric conditions. 


(II) FRee-FLIGHT EXPERIMENTS 


(A) Description 

Each test vehicle consisted of a 14-cal. ogival nose 
adjoining tangentially a cylindrical body of 4-in. 
diameter, the base of the cylinder being flat and closed. 
Four tail fins of rectangular plan form, biconvex section, 
as shown in Fig. 1, were added so that the center of 
pressure of the missile was located well aft of the 
center of gravity, thus providing high flight stability. 
Several missiles of each of the following overall lengths 
were flown: 24.5, 30, 40, 60, 80, 100, and 120 in. 

Four-channel telemetering was installed in each 
missile. Each channel was used to record a pressure 
through a 0.040-in. opening on the missile surface. 
In most of the tests the pressure data recorded were: 
(a) stagnation pressure at the missile nose tip (point 
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Fic. 1. Free-flight missile. Pressure data telemetered were 
(1) stagnation pressure at Pitot tube; (2) differential pressure 1 
in. back of nose tip; (3) on most flights, base pressure at point 
0.5 in. from base center; (4) on base pressure survey missiles, 
base pressure at points 0, 0.5, 1, 1.5 in. from base center. All 
pressure holes 0.040 in. diameter. 


“1” on Fig. 1), to provide Mach Number for com- 
parison with optical tracking data; (b) the difference 
in pressures at two points located 180° apart and 1 in. 
back of the nose tip (point ‘‘2” on Fig. 1), as a measure 
of angle of attack; (c) pressure on the cylindrical body 
for use in another study not included in this report; 
and (d) pressure at one point 0.5 in. from the center 
of the base (point ‘‘3’’ on Fig. 1). In some of the tests 
all four channels were devoted to the measurement of 
base pressures at points located at the center and at 
0.5, 1, and 1.5 in. from the center, all points being in the 
plane of a tail fin (points “4” on Fig. 1). Signals 
were continuously transmitted by a frequency modu- 
lated subcarrier system to a receiving station on the 
ground. This, in turn, provided a continuous record of 
each of the four telemetering channels, together with 
an accurate time record, throughout flight. The 
reliability and accuracy of the telemetering system 
have been well established, and the data from the 
tests reported herein have been found to be consistent 
and reproducible. In addition to the usual laboratory 
calibration of each channel, a spot calibration was 
made on the firing ramp prior to flight, and an auto- 
matic check for zero shift was made in flight. 

The missiles were launched by rocket boosters. 
Two techniques were used; one gave maximum veloc- 
ities corresponding to a Mach Number of about 1.8, 
while the other gave a peak Mach Number of about 2.5 
At the end of rocket thrust the missiles detached from 
the boosters and coasted on in free flight, continuously 
decelerating with a low rate of spin, less than 2 r.p.s. 
The time of separation, clearly shown on tracking and 
telemetering records, marked the start of usable data. 
The terminal Mach Number was about 0.7. 

An optical tracking system was used for accu- 
rately locating the trajectories in space. It consisted 
of four phototheodolite stations, which simultaneously 
recorded the elevation and azimuth of the missile at 


1/,-sec. intervals throughout flight. The position in 
space versus time was obtained by triangulation. By 
differentiating these position-versus-time data, velocity 
as a function of time was found, and with meteoro- 
logical data that provided the local atmospheric pres- 
sure and temperature versus altitude, the velocity 
was converted to Mach Number. (The maximum 
altitudes reached were low, so that data corrections for 
variations of local pressure and temperature were 
small.) The optical tracking system proved to be 
reliable, so that in most of the flights the optical 
tracking data were used for the basic determination 
of Mach Number and dynamic pressure. This system 
has been developed and used successfully in flights 
where the maximum velocity encountered was about 
2,000 ft. per sec. In some of the tests reported here, 
where the maximum velocities were about 3,000 ft. 
per sec., the missiles were not tracked; hence, the basic 
Mach Number data were derived from the telemetered 
stagnation pressure records. It was the practice to 
make all launchings at low angles, so that the missiles 
encountered relatively small variations in atmospheric 
pressure and density during flight. The maximum 
altitude in low-velocity cases was less than 2,500 ft. 


In most of the tests, determination of Mach Number 
from both optical tracking and telemetered pressure 
data provided a check of the accuracy and consistency 
of the two measuring systems. Other checks of internal 
consistency were (1) agreement to within 1 per cent of 
position-time data computed by triangulation from 
each of three pairs of phototheodolite stations; (2) 
a direct check in a few cases of optically determined 
velocities by a continuous wave “‘doppler’’ radar system 
to within 1 per cent agreement, differences being 
random; (3) a sharp decrease in the base pressure as 
shown by the telemetering records when the missiles 
passed through sonic velocity; and (4) a dynamic 
check during flight for any zero shift of the telemetering 
channels. ‘ 

Notwithstanding the care and effort that were given 
to the preparation and testing of each round, a con- 
siderable number of missiles were flown for which the 
data were not valid. Sixteen flights? out of 38 tests 
gave completely usable data, and these form the basis 
of the present report. The remainder could not be 
used because of (a) launching failures, (b) poor or in- 
sufficient optical tracking, or (c) failure of the channel 
or channels telemetering base pressure. In most of 
the tests considered as failures, the data that were 
obtained indicate at least qualitative agreement with 
the reported results. 


(B) Data 

The data collected for each missile were so extensive 
that a recapitulation of all the experimental results is 
not attempted here; however, a sampling of the data 
is given in graphical form for a typical missile (see 
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Figs. 2, 3, and 4). The method of obtaining velocities 
and Mach Numbers has been indicated in Section A 
above, and a plot of the velocity as a function of time 
is given for a specific test (AH-302) in Fig. 2. In this 
figure, optical tracking data and stagnation pressure 
telemetering are compared. 

The pressures measured in flight at the base of the 
missile were the differences between the absolute 
pressure on the base and a reference pressure. Since 
local atmospheric pressure varied slightly as the missile 
gained altitude, trajectory and aerological data were 
used to make small adjustments in the telemetered 
base pressures. These values of base pressure, repre- 
sented by p», are therefore the differences between 
absolute pressure on the base and local free-stream 
static pressure. Typical curves of the measured base 
pressure versus time of flight are shown in Fig. 3. 

In order to reduce the data to dimensionless form, 
the values of p), were divided by dynamic pressure, 
qg(= ‘'/2 pV). The values of the base pressure co- 
efficient, p»/g, represent the basic data derived in this 
experimental work, and a sample is given in Fig. 4. 


(C) Results 


(1) Base Pressure Versus Missile Length.—Within 
the scatter of the data, the tests show no dependence of 
base pressure on missile length in the range of test 
Mach Numbers 1.0<1/<1.7. In Fig. 5 base pressure 
coefficients, p,/g, are plotted against missile length for 
several Mach Numbers in the range of flight data. 
The straight line drawn through the points at each 
Mach Number represents the average of all lengths 
for that Mach Number. One reason for plotting the 
base pressure coefficient as a function of length rather 
than of Reynolds Number is that, although all the 
tests were made under essentially the same atmospheric 
conditions with no significant variation in kinematic 
viscosity, change in length altered the geometry of the 
missiles; hence, the base pressures are not directly 
comparable from missile to missile on a Reynolds 
Number plot. The variation of Reynolds Numbers, 
computed on the basis of total length, was, for example, 
approximately 1.6 X 107 to 7.9 X 107’ at M = 1.2, and 
2.2 X 10’ to 1.1 X 10% at M = 1.6, these values corre- 
sponding to the missile length range of 24.5 to 120 in. 

(2) Dependence of Base Pressure on Mach Number.— 
In view of the fact that no dependence of base pressure 
on missile length was found, the base pressure data 
from individual tests, as measured at points 0.5 in. 
from the center of the base, were averaged. In Fig. 6 
the average base pressure coefficients, p,/q, were plotted 
against Mach Number. In the range 0.8<M<1.8, 
data from all 16 flights are included in the averages. 
Two of the 16 missiles were launched to a Mach Number 
of 2.5; consequently, data shown above M = 1.8, 
represented by a dotted line, are averages from these 
two tests only. The width of the band defined by 
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Fic. 2. Velocity versus time. 24.5-in. test missile, No. 


AH-302. Data obtained by phototheodolites and stagnation 
pressure telemetering are compared; absolute accuracy of both 
methods decreases with decreasing velocity. 
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Fic. 3. Base pressure versus time. Telemetered base pres- 
sures for two free-flight tests (AH-302 and 296). The plotted 
base pressures have not been corrected for variation of local 
atmospheric pressure. The pressure measurement point is 0.5 
in. from the center of the base on both missiles. The Mach 
Numbers are indicated at points along the records. Note be- 
havior at M = 1.0and M = 1.7. 
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Fic. 4. Base pressure coefficient versus Mach Number. 
Typical data on one missile (24.5-in. missile, No. AH-302). 
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Fic. 5. Base pressure coefficient versus missile length. The 
plotted points are taken from free-flight data, and the heavy 
lines represent average values. Point of measurement in all cases, 
25 per cent of radius from center of 4-in. diameter base. 
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Fic. 6. Base pressure coefficient versus Mach Number. The 
solid line is the average of 16 missiles with lengths from 24.5 to 
120 in. where ./ < 1.8; the dotted line is the average of the two 
missiles that attained a Mach Number of 2.5. The values shown 
for o are the standard error of the mean. 


dashed lines in the range 0.8<.1/<1.8 represents twice 
the standard error of the mean at each Mach Number. 
The reason for the increase in o toward lower Mach 
Numbers lies in the nature of the measuring equipment. 
Missile velocities, used in computing dynamic pressures, 
qg, were obtained from phototheodolite tracking data. 
The missile trajectories were such that peak velocity 
occurred in the region of maximum tracking accuracy, 
and, as the missiles flew on, the distances from the 
phototheodolites increased and accuracy correspond- 
ingly decreased. Since the velocity values are squared 
in computing g, the errors in g are double the errors in 
velocity. 

A sharp change in the base-pressure coefficient curve 
(Fig. 6) is found at M = 1, while a more gradual change 
in the curve occurs in the region where M = 1.7. 
The appearance of the change in the actual data in 
these Mach Number regions is illustrated in Fig. 3, 
where the unreduced data from the telemetering system 
are plotted against flight time. The averaged p/q 
plot of Fig. 6 is based on pressures measured 0.5 in. 


from the center of the base. In view of the fact that 
no significant dependence of base pressure on position 
on the base has been found (see below), the plot may 
be regarded as representative of any measurement 
point on the base, except perhaps near the edge. 

(3) Base Pressure Distribution.—Base pressure sur- 
veys were made on test missiles 40, 80, and 120 in. in 
length. On the 40-in. missiles one of the four tele- 
metering channels was devoted to stagnation pressure 
measurement, while the other three pressures were 
measured at points located 0.5, 1, and 1.5 in. along a 
radius from the center of the base. On the S0- and 
120-in. missiles no stagnation pressures were obtained, 
the fourth channel being used to transmit pressure 
at the center of the base. 

Fig. 7 presents pressure data on a typical flight in 
which base pressures were measured at four points. 
There appears to be at most a small variation in base 
pressure with position. The base pressure shows a 
slight tendency toward atmospheric pressure near the 
rim of the base.* In general, the total change from 
center of base to the outermost measurement point is 
less than the magnitude of the estimated probable 
error; for such practical purposes as base drag com 
putation, the base pressure may be taken as inde 
pendent of the position of measurement. The in 
dicated probable errors in the data in Fig. 7 are illus- 
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Fic. 7. Base pressure coefficient versus position. Pressure 
distribution on the base of 80-in. missile, No. AH-290. Indi 
cated probable errors, typical of free-flight missiles, are due to 
errors of velocity and telemetering measurements. Data at 
M = 1.73, marked with an asterisk, are wind-tunnel measure- 
ments on 24.5-in. free-flight model at Daingerfield; all other data 
are free flight. 
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Fic. 8. Comparison of base pressure coefficient with large 


angles of attack in the transonic region. Note the differential 
pressure is measured 1 in. back of nose tip, between two points 
180° apart, on 40-in. missile, No. AH-280, and is proportional 
to angle of attack——the factor of proportionality is a function of 
Mach Number. 


trative of all flights, combining the errors in velocity, 
in telemetering, and in data reduction. 

(4) Relation Between Angle of Attack and Base 
Pressure.—Free-flight data indicate that base pressure 
is at most only slightly affected by small angles of 
attack of missiles at supersonic speeds. While flights 
in which angles of attack were less than 1° were used in 
obtaining the data in Fig. 6, in one or two flights yaws 
Base pressure data 
There were not 


of several degrees were observed. 
on these showed only small variation. 
sufficient flight tests where yaws of considerable magni- 
tude occurred to provide a correlation of base pressure 
with angle of attack. In one test, however, where 
there was large flight instability, apparently due to 
damage of the missile in launching, it was interesting to 
note that the sensitivity of base pressure to angle of 
attack increased sharply in the transonic region. In 
Fig. 8 are plotted base pressure coefficients and nose 
differential pressures versus flight time for this case. 
Its short-duration supersonic flight was evidently at a 
considerable, although steady, angle of attack, and the 
base pressure coefficients, »,/g, were much larger than 
normal. When the missile passed through sonic veloc- 
ity, it started to oscillate irregularly, and base pressure 
data were markedly affected, in phase with the yaw. 


(III) Wunp-TuUNNEL TESTS OF A FREE-FLIGHT MODEL 
(LoNE STAR LABORATORY) 


(A) Description of Experiments 


A model of the same dimensions as the 24.5-in. 
missile used in free-flight tests was studied at M 
1.73 in the Lone Star Laboratory wind tunnel at Dain- 
gerfield, Tex. The model was supported by a 0.752- 
in. diameter sting at the center of the 4-in. diameter 
base. There were eight points of base pressure meas- 
urement, four along each of two radial lines 45° apart, 
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one row being in the plane of a tail fin. The pressure 
taps were located at points 0.25, 0.5, 1, and 1.5 in. 
from the rim, the innermost tap being only 0.125 in. 
from the sting surface. Data were obtained up to 
effective angles of attack of about +5° (see Fig. 9). 
The pressure distribution over the ines was obtained 
by combining measurements taken with the model in 
various positions of rotation about the sting axis. 


(B) Results 


(1) Dependence of Base Pressure on Mach Number.— 
The wind-tunnel data give an evaluation of base 
pressure coefficient at a single Mach Number—in 
this instance, 1.73. This is to be compared, of course, 
with the flight-test results. In Fig. 10 the average 
curve of base pressure versus Mach Number from free- 
flight data (Fig. 6) has been reproduced. The average 
base pressure coefficient at zero angle of attack as 
measured on the free-flight model in the Daingerfield 
tunnel, P,/g = —0.21, is also plotted on Fig. 10 at 
M = 1.73. The latter appears significantly smaller 
than the free-flight values at this Mach Number. 
Several possible explanations may be offered for this 
difference : 

(a) There may have been a pronounced sting effect, 
as already suggested. 

(b) Schlieren observation of the flow about the wind 
tunnel model at JJ = 1.73 showed that the bow shock 
reflected from the top and bottom walls of the tunnel 
just barely cleared the base of the model and inter 
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ANGLE OF ATTACK, (DEGREES) 


Fic. 9. Base pressure coefficient versus angle of attack. 
Variation of base pressure with angle of attack, measured on 
24.5-in. missile at \f = 1.73, in Daingerfield wind tunnel (Lone 
Star Laboratory). Notice behavior of point near sting. 
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MACH NUMBER, M 
Fic. 10. Base pressure coefficient versus Mach Number. 
Comparison of free-flight, wind-tunnel, and ballistic data, all 
at zero angle of attack. von Karman-Moore semiempirical 
theory is also shown. 


sected the sting a short distance back of the base. 
Since the flow in the wake is almost certainly subsonic 
for some distance downstream of the base, the effect 
of this shock may be communicated to the base. Fur- 
thermore, the reflected bow waves from the sides of 
the tunnel wall would be reflected to the sting even 
further upstream, because the width of the tunnel was 
only 19 in. as compared to the 27.5-in. height. Sub- 
sequent tests with the same model at M = 1.5, where 
the shocks clearly intersected the model slightly ahead 
of the pressure to be 
greatly affected. The tests at WM = 1.5 were dis- 
continued because of the large effect of the reflected 


base, showed the _ base 


shocks. 

(c) Physical conditions in the tunnel test section 
were different from those encountered in free flight. 
In particular, Reynolds Number in the test section 
was 1.42 X 107 as compared to about 2.3 X 107 for the 
same model in free flight at JJ = 1.73. It seems un- 
likely, however, that the 1S per cent difference in base 
pressure coefficient (—0.255 —0.21) 
could be due entirely to a Reynolds Number effect. 
Another difference between the physical conditions 
in wind tunnel and free flight was humidity, which may 
The dewpoints during the tests 


compared to 


contribute an effect. 
were between 0° and — 10°F. 

(2) Base Pressure Distribution. —The base pressure 
is only slightly dependent on position on the base; 
the coefficient decreases in absolute value toward the 
rim to about the same extent as observed in free flight 
(compare the a = 0° intercepts of the three regular 
curves in Fig. 9). 

(3) Relation Between Angle of Attack and Base 
Pressure.—Variation with attack up to 
about +5° is also small—of the same order of magni- 
tude as the change with position on the base. (Ab 


angles of 


solute value of f»/g increases with angle of attack.) 
This is in agreement with qualitative conclusions from 
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free-flight tests. Moreover, the effect appears to be 
symmetric with respect to angle of attack and of the 
same order of magnitude wherever the pressure is 
measured, except near the sting (Fig. 9). Nevertheless, 
it seems reasonable that the effect may become large 
and asymmetric with respect to the center of the base 
at large angles of attack. a 

(4) Effect of Tatl Fins on Base Pressure.—The 
proximity of tail fins has no apparent effect on base 
pressure; pressures measured along two lines on the 


base—namely, in the plane of a fin and 45° out of the 
plane—showed no difference in any of the measure- 
ments. 


(5) Effect of Base Sting on Base Pressure.—The 
mounting sting has a strong effect, at least locally, on 
base pressure. The two pressure taps located near the 
sting (0.125 in. away) showed marked variation with 
angle of attack and roll, whereas the taps further 
away showed only small variation. In Fig. 9 are 
plotted base pressure coefficients, p,/q, versus angle of 
attack for the four radial positions of measurement 
which were in the plane of the angle of attack. The 
erratic behavior of pressure near the sting is clearly 
shown; furthermore, it is possible that the sting may 
have had a consistent effect on all the measurements, 


but this is not evident from the data. 


(IV) COMPARISON WITH OTHER PRESSURE DaTA 


Two other types of base pressure measurements are 
available for comparison with the data reported here. 
They are (1) wind-tunnel measurements on relatively 
small models at Reynolds Numbers that are extremely 
low compared to free flight, so that there may be a 
scale effect; and (2) ballistic determinations, where 
the base pressure is not measured directly but is 
derived from total drag measurements together with 
theoretical and/or experimental estimates of drag 
components other than base drag. The work de- 
scribed below is not all that is available; unfortunately, 
some of the data are classified and therefore cannot be 
discussed here. It is believed that no other direct 
measurements of base pressure in supersonic flight have 
been published. 


(A) Wind-Tunnel Results 


The results of measurements made in several wind 
tunnels will be discussed and may serve as examples 
of the trend of results obtained. 

(1) Measurements Made in the Kochel Wind Tunnel. 
—These tests, as reported recently,* were made on 
small cone-cylinder models where the cone angle varied 
from 20° to 60°. Measurements were made at Mach 
Numbers of 1.2, 1.56, 1.84, 2.5, 3.1, and 4.38, but there 
is considerable doubt as to the validity of measure- 
ments at the high Mach Numbers (M>3). Therefore 
the data shown for comparison in Fig. 10 are for the 
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It should be noted that 


lower Mach Number values. 
In view 


the agreement at 1 = 2.5 and less is good. 
of the fact that the Kochel type of tunnel induces from 
the atmosphere, the Reynolds Numbers may be com- 
puted with reasonable assumptions for atmospheric 
pressure and temperature. For the range of reported 
Mach Numbers, the Reynolds Numbers should have 
been of the order of | X 10®to5 X 10°. Because of the 
meagerness and large scatter of the data, and especially 
since the wind-tunnel Reynolds Numbers were lower 
than free-flight values by an order of magnitude, the 
differences may not represent significant disagreement. 
It should be noted (Fig. 10) that the Kochel wind- 
tunnel studies at widely scattered Mach Numbers have 
failed to reveal either the fairly constant value of 
p,/qin the range 1<./< 1.6 or the marked changes in the 
regions of JJ = 1 and of M = 1.7. Kochel data, in 
fact, have been interpreted as indicating a continuous 
decrease in absolute value of »,/q with increasing Mach 
Number. Since all the wind-tunnel data were taken 
at Reynolds Numbers that were smaller than those 
obtaining in free flight, a comparison of p)/q versus 
logarithm of Reynolds Number is given in Fig. 11. 
The data presented indicate the trend toward an in- 
crease in absolute magnitude of base pressure co- 
efficient with increasing Reynolds Number. It is true, 
of course, that these data are not strictly comparable 
on a Reynolds Number basis because the models com- 
pared are not homologous and, hence, do not meet the 
requirements for scale-effect comparison. Neverthe- 
less, at supersonic velocities where the fineness ratio, 
length/diameter, is 5 or greater, it is not expected that 
total length and small differences in nose shape should 
affect the base pressure to a large extent. Thus, while 
the data in Fig. 11 correspond to a wide range of fine- 
ness ratios, they seem to be the best representation of 
the results that are now available in regard to possible 
Reynolds Number effects. 

(2) Guidonia Wind-Tunnel 
sults of base pressure measurements on three projectile 
shapes in the Italian Guidonia wind tunnel at a Mach 
Number of 2.06 are shown in Fig. 10. In this case, 
again, the models were small (1.2 in. diameter) cone- 
with a_ bourrelet. 


Measurements.*—Re- 


cylinders and _ ogive-cylinders, 
Data on the test conditions in the tunnel are rather 
meager, although the Reynolds Numbers must surely 
have been less than free-flight values. The average 
values of base pressure coefficients as read from the 
figures given by Ferri® are plotted in Fig. 10. This 
value is considerably different from flight-test data, and 
the difference again may be due in part to scale effect. 
Plots of the Guidonia data® gave a change of about 
—0.01 in p,/q with 4° angle of attack; this, together 
with the fact that measurements at each point changed 
in the same way, is in essential agreement with the 
Daingerfield tests. In addition, the independence of 
base pressure on position is verified by the measure- 


ments at a large number of points. 
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Fic. 11. Base pressure coefficient versus logarithm of Rey- 
nolds Number. The Kochel wind-tunnel results are shown for 
WU = 1.56 and M = 1.84 because no measurements were made 
at M = 1.73, the value for the Daingerfield wind-tunnel meas- 


urements and free-flight tests. The fineness ratio varied from 
model to model for the Kochel tests over a range 1 < length + 
diameter < 6; for the Daingerfield tests, length/diameter = 6; 
for flight tests, 6 < length/diameter < 30. The average value 
of the flight tests is transposed from Fig. 6 for comparison. 


(B) Ballistics Data 


Ballisticians’ and others have computed base drag 
and base pressure coefficients by considering the total 
drag of a simple projectile, such as a cone-cylinder com- 
bination, to be made up of the wave drag of the cone, 
the skin friction drag of the conical and cylindrical 
surfaces, and the base drag. Values of the first two 
components, computed theoretically or obtained from 
wind tunnels, are subtracted from the measured total 
drag, the remainder being attributed to base drag. 
In this way, Charters’ has obtained base pressures, and 
his values, P,/): (where /; is free stream static pres- 
sure and P, is absolute base pressure), have been 
converted for purposes of comparison here to the 
coefficient form p,/q (by taking dynamic pressure 
q = '/eyM*p,). The converted data are plotted in 
Fig. 10. It seems likely that disagreement of this 
order of magnitude might arise from the method of 
computation—i.e., errors in the value of wave drag 
and skin friction, computed or measured in a wind 
tunnel, might combine in the subtraction from total 
drag to give a base drag in proportionately greater 
error. It should be noted that the ballistic data tend 
toward good agreement with the flight-test results in 
the neighborhood of M = 2.5. 


(C) Theory 
There has not been published any satisfactory theory 
to predict the values of base pressure. However, in 
1932, von Karman and Moore*® developed a semiem- 
pirical formula that gave the base pressure coefficients 
as 
b/bo = (1 — [(y—1)/2)AM2Y"/O- 


where p is the absolute value of the pressure at the 
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base, fo is free-stream static pressure, y is the ratio of 
specific heats, and A is an arbitrary constant to which, 
on the basis of low subsonic measurements, von Karman 
and Moore assigned the value 0.2. This function, 
converted to p,/q is plotted in Fig. 10 up to M = 2.5 
(the limit of applicability given by von Karman and 
Moore) for comparison with the free-flight, wind-tunnel 
and ballistic data already described. It is not sur- 
prising that the agreement is not particularly good, 
in view of the rather oversimplified nature of the 
theory. There is obviously a need for an adequate 
theory of the pressure acting on the base of a missile 
of simple shape, and it is hoped that the experimental 
data from this study will stimulate interest in a theoret- 
ical development and afford a means of checking 
theoretical predictions. 
(V) CONCLUSIONS 

(A) The magnitude of the base pressure coefficient 
under conditions prevailing in free flight has been 
determined and recorded in the range 0.9<_M<2.5. 
There exists an unresolved disparity between these 
results and those of wind-tunnel and ballistic measure- 
ments at low Mach Numbers, whereas in the neighbor- 
hood of 4 = 2.5 the agreement is good. Further 
studies are needed on the effects of scale, which may in 
large part explain the difference and establish an explicit 
dependence of base pressure on Reynolds Number. 
The combined evidence of wind-tunnel and _free- 
flight tests thus far indicates that the base pressure co- 


efficient increases in absolute magnitude with in- 
creasing Reynolds Number. 
(B) Base pressure is independent of body length 


for the missiles considered in this study (6 < fineness 
ratio < 30). Thin streamlined airfoils mounted on the 
body have little or no effect on base pressures. 

(C) 
to a good approximation, over the base of a missile 
There is little effect 


In supersonic flight, base pressure is constant, 


(except possibly near the rim). 
of angle of attack up to about 5°, and for practical con- 
siderations this effect may be ignored. 

(D) Asa result of this study the need for certain 
specific free-flight tests becomes evident; these are: 
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(1) Tests at high altitudes where the Reynolds 
Numbers obtained can be lower than at sea level and 
thus provide data nearer the range of wind-tunnel 


tests. 

(2) Tests at higher velocities at sea level, thereby 
extending the present data to higher values of both 
Mach and Reynolds Numbers. 

(3) Tests with homologous missiles, to study scale 
effect at high Reynolds Numbers and to obtain further 


correlation with low Reynolds Number tests in wind 


tunnels. 

(4) Tests with various degrees of smoothness of the 
missile surface, to examine the possible effect of the 
condition of the boundary layer at high Reynolds 
Numbers. 

(E) Finally, there is the need for a theory to inter- 
pret the existing results and the data to be obtained in 
(D) above and to present a unified picture of the 
effects of supersonic flow past the base of a missile. 
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The Analytical Design of an Axially 
Svmmetric Laval Nozzle for a 


~ Parallel and Uniform Jet 


KUNO FOELSCH 


Consulting Physicist 


ABSTRACT 


The equations for the nozzle’s contours are derived by integra- 
tion of the characteristic equations of the axially symmetric flow. 
Since it is not possible to integrate these equations mathematic- 
illy in an exact form, it was necessary to find a way to approxi- 
mate the calculations. The approximation offers itself by con- 
sidering and comparing the conditions of the flow in a cone with 
those in a nozzle, as a linearization of the characteristic equations. 

he first part of the report deals with equations for the transi- 
tion curve by which the conical source flow is converted into a 
parallel stream of uniform velocity. The equations are derived 
by integration along a Mach line of the flow in the region where 
the conversion takes place. A factor f is introduced expressing a 
relation between the direction and the velocity of the flow along 
, certain Mach line. f remains undetermined and is not involved 
in the final equations. 

In the second part of the report, the spherical sonic flow sec- 
tion is converted into a plane circular section of the throat. 
The nozzle’s contour adjacent to the throat is formed by the arc 
of a circle connected with the transition curve by a straight line. 
The gas dynamic properties of the boundary Mach line are cal- 
‘ulated in Table 1, the use of which shortens the calculations con- 


siderably. 


SYMBOLS 


Vf = Mach Number of the flow 
x = Mach angle of the flow 
= velocity 
tp = mass flux of the flow 
9 = angle between direction of flow and axis of symmetry 


= semicone angle 


yr = radius vector from source point 

D = diameter of nozzle at exit 

s; = radius of the throat area 

R = radius of the nozzle’s wall curvature at the throat 
\ = length of a straight part of nozzle’s contour 

¥ = expansion angle of a Mach line 

u = asmall increment of Mp 

€ = asmall increment of 6p 

r =r/ro 

L = length of the supersonic part of the nozzle 


y = ratio of specific heats; f and k constants along a Mach line 
Subscripts 
P refers all quantities to a variable point along the Mach curve 
1 = fixed values at point A 
E = final values of all quantities at point Z 
) = values at the critical arc in a cone flow 


INTRODUCTION 


I AVAL NOZZLES PRODUCING parallel and uniform flow 
are a necessity for wind-tunnel testing and are 
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often required for jet propulsion units. The design of a 
two-dimensional nozzle (plane flow) has been treated 
in the past by several authors. In 1929 Prandtl and 
Busemann! worked out their famous graphical method, 
which has been much in use since. Later, Puckett* 
modified their graphical method to a semigraphical 
one. The author’ of the present report proposed at the 
J.P.L. Cal-Tech symposium in 1946 an analytical de- 
sign method which is based on the fundamental ideas 
used by Prandtl and Busemann in solving the problem 
graphically. 

However, the: case of a three-dimensional nozzle 
with axial symmetry represents a more difficult problem. 
Up to date there are known in the literature only graph- 
ical design methods as proposed by Sauer,‘ Guderley,® 
and de Haller. All these proposed methods for de- 
signing an axially symmetric nozzle are, for practical 
use, tedious and complicated. They require a great 
deal of graphical work, since two charts are made and 
used simultaneously and involve errors of uncertain 
magnitude, particularly when the graphical construc- 
tion approaches the axis of symmetry. 

It is the purpose of this report to show that a good 
approximation for an analytical solution of the prob- 
lem may be obtained. In other words, the coordinates 
of the contour of an axially symmetric nozzle, as well 
as of any streamline in the nozzle’s flow, may be im- 
mediately determined from two simple equations. 

The first part of the report deals with the derivation 
of the equations for the transition curve by which the 
conical source flow emanating from the sonic section 
is converted into a parallel and uniform flow. In the 
second part of the report the spherical sonic flow section 
is adjusted to a plane flow section at the throat of a 
nozzle. 


CONVERSION OF A RADIAL FLOW INTO A 
PARALLEL AND UNIFORM FLOW 


Part I. 


Fig. 1 shows a section through the upper supersonic 
part of a three-dimensional nozzle with axial sym- 
metry. The flow in the nozzle may be divided into three 
regions as to the character of the flow. 


(I) Region TREA: Radial flow. 
(II) Region AEB; Conversion of flow into parallel 
stream by transition curve AB. 
(III) Region BEX: Parallel and uniform flow. 
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Fic. 1. Conversion of radial flow into parallel flow. 


Analysis of Flow in Region I 


The flow enters the nozzle at the arc 7‘R with sonic 
velocity and proceeds in radial streamlines until the 
Mach line AF is reached. The expansion of the flow 
along the Mach lines in that region is given by the 
equations [see Eq. (24)] 


d@ = +(dw/2w) cot a (1) 


At point E the expansion is completed, and the flow is 
at the state the nozzle is designed for. 


Analysis of Flow in Region II 


Region II is bounded by the curved Mach line AL, 
the straight one HB, and the transition curve AB. 
A left Mach line PA, starting at any point P along ABZ, 
is described by the characteristics equations of an 
axially symmetric potential flow: 


dy/dx = tan (0+ a) (2a) 
dy/y = [(dw/w) cot a — dé] (cot a+ cot @) (2b) 


The above equations represent two relations between 
four quantities w, 6, y, and x, and the functions w(x), 
6(x), and y(x) along a Mach line may be found only if 
there is another equation stating, for example, a cer- 
tain relation between w and @ for a certain Mach line as 


F(w, 0, 0p) = O (3) 


Substituting Eq. (3) into Eqs. (2a) and (2b), we obtain 
the three above-mentioned functions by direct or step- 
by-step quadratures. 

If the nozzle problem is solved by the methods of 
characteristics, then there is no need of Eq. (3), since 
the geometric and gas-dynamic properties of a Mach 
line are established by using all four characteristic 
equations simultaneously in a slow graphical integra- 
tion process. 

Before we analyze Eq. (3), let us first consider the 
flow in the nozzle after passing the reflected Mach line 
AE and compare it with the flow in a cone. At point 
A the straight cone wall has been bent into the nozzle 
wall AB (transition curve) by which the radial stream- 
lines at AE are being turned parallel. The magnitude 
of the velocity that would have been acquired by the 
flow in a cone along the circular arc LE will be attained 
by the flow in the nozzle only after passing the line EB. 
It is obvious that the process of expansion in the nozzle 
has been considerably slowed down and that the ve- 


locity vector along PK will not differ much from the 
velocity vector at point P as to its direction and mag- 
nitude. This statement justifies the below performed 
integration of Eqs. (2a) and (2b) in a linearized férm. 

For the analysis of Eq. (3) we make the assumption 
that, along a continuous Mach line PK, the velocity x 
is uniquely determined by the direction 6 of the flow. 
Since we are dealing with an expansion and potential 
flow and since the Mach line PK is continuously inter- 
sected by Mach lines of the second family, we conclude 
that, along PK, dw must be positive and dé either posi- 
tive or negative but not zero. However, along a stream- 
line down the flow in Region II, both quantities ap- 
proach zero, a value that is attained at the straight 
Mach line BE. 

We consider for Eq. (3) the possibility of two ex- 
treme cases with the actual case somewhat between 
them. 

The assumption that d@ is positive along PK suggests 
to consider Eq. (1) with the upper sign as one extreme 
case for the differential form of Eq. (3). Substituting it 
into Eq. (2b), we obtain 


dy/y = (1/2)(dw/w) cot a (cot a + cot @) (4 


Eq. (4) indicates a fast expansion along the curve PA, 
since it is the equation for a Mach line in a conical 
source flow, as can be verified immediately from the 
properties of the flow in a cone (Fig. 2). 


y =rsin 8, dr/rd@ = cot @ 


The assumption that dé is negative along PKA leads t 
the other extreme case for Eq. (3): 


d@ = —(1/2)(dw/w) cot a (la 


Eq. (la) shows a much slower expansion along PK than 
in the first case, as can be seen by substituting it into 
Eq. (2b): 


dy/y = (3/2)(dw/w) cot a (cota + cot @) (5 


It is of interest to notice that Eq. (la) represents the 
characteristic equation for the curve AE, and, by apply- 
ing it to the curve PK as the differential form of Eq. 
(3), each point of PA is given an image point along PF 
as to the magnitude and direction of the vector w. 

If AK, is the corresponding point to A (Fig. 1), the 
entire curve PK is gas-dynamically reflected in Curve 











Fic. 2. Mach line PK in a conical source flow. 
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PK;, and, further, each point of the transition 
curve AB has gas-dynamically an image point along 
AE. 

We assume now that Eq. (3) has the differential 


form 


dé = k(dw/w) cot a (6) 


with +(1/2) > k > —(1/2), and, by substituting it 
into Eq. (2b), 


dy/y = f(dw/w) cot a (cot a + cot 6) (7) 


with (3/2) > f > (1/2). f is assumed to be constant 
ilong a Mach line but may vary from Mach line to 
Mach line as a function of the boundary conditions along 
AE. 


Calculation of the Coordinates of the Transition Curve 


Expressing the velocity w and the Mach angle by 
the Mach number J/ of the flow, Eq. (7) is 

dy _[(M?—1) + VM? — 1 cot 6]dM a 

= f-—— ever eae a 

y M{1 + ((y — )/21a} oe 

Since it was shown above that the values of M and 6 

at any point of curve PK do not differ much from those 
of Mp and 6p, respectively, at point P, we set 


M= Mp+4, 6=Op+e (8) 


und consider » and ¢ as small quantities so that their 
higher powers and products yp’, e*, we ... may be neg- 
lected in our calculations. Observing that rule for 
the integration of Eq. (7a) from point P to any arbi- 
trary point N on the curve PK, we obtain 


y = a peaaiiaatad (9) 


with 
V Mp? — 1 (WV Mp? — 1 + cot 6p) 


Mp{1 + [(y — 1)/2] Mp3} 





A line element dy of the curve PK is 
dn = dy/sin (a + @) 
or with Eq. (9) 
dn = fAy(MdM/cos 06 + VM? — 1sin@) (10) 


To compute the y coordinate for the intersection point 
K of the Mach curve PK and the transition curve AB, 
we apply to the flow in Region II the law of conserva- 
tion of mass, which states that the mass flow that passes 
the surface of the spherical zone at critical section be- 
tween the angles w and 6p and with the radius 7) must 
also pass the surface formed by rotation of PK about 
the X axis. 


@ K 
2a(wn)art ff sin 0 d@ = on f (wp) sin aydn (11) 
Oe P 


The above equation determines also the y coordinate 


for any streamline if w and K are changed correspond- 
ingly. 

The mass flux (wp) at point N may be expressed by 
the Mach Number as (reference 4, Eq. 28): 


9 ¥— 1 (y+1)/2(7—-1) 
wp = (wp) u/(— — $. ——_— ur) 12) 

Ny tL yt 
Substituting Eqs. (9), (10), and (12) into Eq. (11) and 
observing the rule we have set up for Eq. (8), the inte- 
gration of Eq. (11) gives 


To (cos Ap -— Cos w) = 
2/A(M— Mp 
eC” ( a l 


1 
ae — (13) 
2°” sin Op(V Mp? — 1 + cot Op)rp’ \ 

with 

> (y+1)/2(y7—1) 
2 y¥-1 
Tp? = | - oe M *) Mp (13a) 
° oer ytio' aes 
rp has a geometrical meaning. It represents the ratio 
of the radius at point P to the radius 7» of the critical 

are (reference 4, Eq. 33). 

Tp = Yp/To (14) 

Furthermore, 


Jr = le sin Op (15) 


With Eqs. (14) and (15), Eq. (13) is simplified to 





ef M—™P) sin Op = F(Op) (16) 

with 
F(@p) Pe = - a eal ade E 
V sin? Op + 2(cos Op — cx S w)(W Mp? — lsin Op + cos Op) 
(16a) 


and by Eqs. (9), (14), and (15) the y coordinate of point 
K is 
y = rotpF (6p) (17) 
ro, the radius of the critical are 7'R, is determined by 
applying the law of conservation of mass to the exit of 
the nozzle with the diameter D (Fig. 1). The flow at 
the exit is gas-dynamically of the same state as in point 
E. 
Aro? sin? (w/2)(wp)o = (xD?/4)(wp)e (18) 
and by Eqs. (12) and (13a) 
(wp)e = (wp)o(1/rx*) 
Substituting it into Eq. (18), we obtain 
to = D/Arg sin (w/2) (19) 
which transforms Eq. (17) into 
y = [D/4 sin (w/2)](te/t2) FOr) (20) 


y, the ordinate of the transition curve, has been ob- 
tained as a function of 6p, since rp and Mp are deter- 
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TABLE 1 i 
PROPERTIES OF BOUNDARY MACH LINE AE (Fie. 1) 
vy = 1.40 
V/M?—1 y° T VM? — 1 y° T 
0.1 0.00788 1.0000 2.65 23.2085 1.8999 
0.2 0.06188 1.0002 2.70 23.6861 1.9428 
0.3 0.2021 1.0008 2.75 24.1563 1. 9867 
0.4 0.4581 1.0024 | 2.80 24.6189 2.0317 
0.5 0.84729 1.0055 | 2.85 25.0742 2.0778 
0.55 1.09387 1.0077 2.90 25.5221 2.1251 
0.6 1.3749 1.0106 | 2.95 25.9626 2.1734 
0.65 1.68968 1.0141 3.00 26 . 3959 2.2229 
0.7 2.0369 1.0182 | 3.05 26 8220 2.2734 
0.75 2.41491 1.0231 3.10 27.2411 2.3251 F 
0.8 2.82200 1.0288 3.15 27.6532 2.3780 
0.85 3.25612 1.0352 j 3.20 28.0584 2.4320 
0.9 3.7151 1.0425 | 3.25 28.4568 2.4871 
0.95 4.19662 1.0506 3.30 28.8480 2.5434 
1.00 4.6987 1.05965 3.35 29.2336 2.6009 P 
1.05 5.21910 1.0696 3.40 29.6122 2.6596 
1.10 5.75559 1.0804 3.45 29.9814 2.7194 
1.15 6.30609 1.0921 3.50 30.3503 2.7804 
1.20 6.86885 1.1048 3.55 30.7100 2.8427 
2 1.25 7.44175 1.1184 3.60 31.0638 2.9061 p 
1.30 8.02319 1.1330 3.65 31.4115 2.9707 fl 
1.35 8.61138 1.1485 3.70 31.7534 3.0366 
1.40 9.20503 1.1650 3.75 32.0894 3. 1037 a 
1.45 9.80251 1.1824 3.80 32.4198 3.1720 Dp 
1.50 10.4026 1.2008 | 3.85 32.7449 3.2462 : 
1.55 11.0041 1.2202 | 3.90 33.0644 3.3124 tl 
1.60 11.6059 1.2406 3.95 33.3786 3.3844 Oo 
1.65 12.2072 1.2619 | 4.00 33.6876 3.4577 : 
1.70 12. 8068 1.2842 | 4.05 33.9914 3.5224 I 
1.75 13.4039 1.3075 4.10 34.2903 3.6082 fi 
1.80 13.9980 1.3318 | 4.15 34. 5841 3.6854 
1.85 14.5883 1.3571 4.20 34.8732 3.7639 tl 
1.90 15.1743 1.3834 4.25 35.1576 3.8436 E 
1.95 15.7553 1.4106 4.30 35.4373 3.9247 
2.00 16.3311 1.4389 4.35 35.7125 4.0071 
2.05 16.9011 1.4682 4.40 35.9832 4.0907 
2.10 17.4651 1.4985 4.45 36.2497 4.1758 
2.15 18.0228 1.5298 4.50 36.5117 4.2621 
2.20 18.5737 1. 5621 4.55 36. 7696 4.3498 fr 
2.25 19.1178 1.5955 4.60 37.0234 +. 4389 
2.30 19.6549 1.6299 4.65 37.2733 4.5293 
2.35 20.1848 1.6653 4.70 37.5192 4.6211 
2.40 20.7073 1.7017 4.75 37.7611 4.7142 
2.45 21.2225 1.7393 4.80 37.9994 4.8106 E 
2.50 21.7303 1.7778 4.85 38.2339 4.9047 
2.55 22.2306 1.8175 4.90 38.4648 5.0020 
2.60 22.7233 1.8582 
mined by the conical source flow (see Table 1). wand Eqs. (20) and (21) determine the transition curve by 
Tg are constants; they will be analyzed below. which a supersonic conical flow emanating from the 
The x coordinate of point K is computed from Eq. spherical surface TR is converted into a parallel and ul 
2a) in the following way: uniform jet. We notice that the unknown factor f ; 
in Eq. (7) is not involved in the final equations for x 
dy dM r 
dx = —————~ = fAy ——— - and y. ¥ 
tan (0 + a) tan (9+ a) ne , 
’ The parameter @p has values from zero to w; 7, and 
™ ; fA(M—Mp) eae ; i : 
il. al V M,*? — 1 are determined from these values as will be 
dx = fAyp aponenmeneie: aE : ; . 
zp up tan (0 + @) shown later. For #p = 0 we obtain from Eqs. (20) and 
——- (21) with 
° ss V Mp? = 1 a tan Op fu MP) It 
x — Xp = yp ————_ ————— ~1] . ; " A 
1 + tan OpV Mp? mt [F(@p)]o,<0 = 2 sin (w/2) and rp = rz 
. ” . [Ylep-o = D/2 
and by using Eqs. (9), (14), (15), (16), (19), and . D D 
- =v [x]op=0 ak a eee + —= O08 Oe V1 
Xp = Yp cot Op ‘ 4sin (w/2) 2 th 
we obtain the x equation as tia , ‘ ‘ M 
q which is in agreement with the geometrical picture al 
x = [D/4 sin (w/2)](tp/Te) X of the nozzle, since, according to Eqs. (14) and th 
1 + (cos pV Mp? — 1 — sin 6p) F(6p) - (19), ti 
sin 0pV Mp? — 1 + cos Op OE = rg = ter = D/4 sin (w/2) 
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Wy =+*NOA + & AOP 


Fic. 3. Definition of an expansion angle in a conical source 
flow, 


Properties of a Conical Source Flow. Relation Between 


M, 7, 0, w, and y 

In a conical source flow a Mach line starting at any 
point of the critical are (point N, Fig. 3) runs across the 
flow and is reflected from the walls. A point moving 
along the Mach line from the start of it to any arbitrary 
position will have covered a certain angle, which we call 
the expansion angle of the Mach line (at the position 
of the point) and denote it by y as it is shown in Fig. 3 
for point P (see also reference 3). The expansion angle 
for a certain point of a Mach line can be expressed by 
the Mach Number. The logarithmic differentiation of 


Eq. (13a) gives 


dt (ll? — 1) dM 





= —-~ — (22) 
tr 2Mi1l + [(y — 1)/2]M*} 
from Fig. 3 and Eq. (14) it is 
dr/r dp = dr/r dp = VIP = 1 (23) 
By substitution of Eq. (23) into Eq. (22), 
dy Vie = 1 (24) 


dM 2Mf1 + [fy — D2} 


and by integration of Eq. (24) along a Mach line from 
the start (YW = 0, M = 1) toan arbitrary point with the 
Mach Number JJ, 


-1 y-i 
- ES We (re -1) - 
21 ¥7r-1 rs 


tan-' VM? — 1 (25) 


It is, by definition of y along the boundary Mach line 
AE of the flow in Region I (Fig. 3), 


Op = ve — vp (26) 


ve is the final (maximum) expansion angle reached by 
the flow at point Z and is computed from Eq. (25) for 
Mr, the design Mach Number of the nozzle. yp varies 
along AE from yp, to Ye. 4 determines the length of 
the nozzle and is chosen by the designer from the rela- 
tion 


Vode S Va < Ve 


The possible shortest nozzle is obtained from 
va = "Jobe = w = 14 (27) 


which, gas-dynamically interpreted, means that the 
boundary Mach line starts at the intersection point of 
the critical arc and the axis of symmetry and is re- 
flected once at the wall, at point A, as it was assumed 
in Fig. 1. The values of Y and 7 computed fron 
Eqs. (25) and (13a), respectively (for y = 1.40), are 
civen in Table 1 as function of VM? — 1, which is a 
more convenient term than JJ since it is involved in 
Eqs. (20) and (21). 


Part II. PLANE FLOW SECTION AT THROAT 


It is usually assumed that the flow section at the 
throat (M/ = 1) of a nozzle is plane and not spherically 
curved as at the critical section of a conical source flow. 

The conversion of the spherical flow section into a 
plane circular flow section may be accomplished by 
bending the portion of nozzle wall adjacent to the 
throat into any smooth convex curve Q7D tangent to 
the rest of the wall at D and having at the throat a 
tangent parallel to the axis of symmetry (Fig. 4) 
For continuity reason the cross-section area at the 
throat must be equal to the area of the spherical section 


ws? = 4mrro? sin? (w/2) 
or by Eq. (19) 
s = D/2rtz (28) 


The change in the geometry of the sonic flow section 
will certainly affect the source character of the flow in 
Region I of the nozzle. However, by inserting a straight 
line into the nozzle contour ahead of the transition 
curve, the source character of the flow may be restored 
and secured in the portion of the nozzle adjacent to 
the Mach curve AE. This is a necessity for our nozzle 
design, since the calculation of the transition curve was 
based on the existence of a radial flow along the Mach 
curve AE. We now refer all coordinates to the new 
origin 0,, at the center of the throat and choose for the 














Fic. 4. Conversion of a spherical flow section into a plane flow 
section at the throat. 
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curve QTD the arc of a circle having radius R and the 
center at 02 (Fig. 4). The equation of the circle is 


y =s+R—RV1 — (x/R)? (29) 
and the coordinates of the point D are 
Xp = Rsin ow; Yo = s+ RI — cosw) (30) 


Through the transformation of the coordinate origin 
from 0 to 0,, the x values for the transition curve must 
be reduced by x, which, derived from Eqs. (29) and 
(28), is 


xo = (D/2rg) cot w — R tan (w/2) (31) 


The y values of Eq. (20) are not affected by the above 
transformation. 

It can be shown that between R and the length A 
of the straight part DA of the nozzles contour exists 
the relation 


.+ Rtan- = 
A+ a5 = 


D ] 
——___|,- | (a2) 
dry sin (w/2) cos (a/2) 


R must be so chosen by the designer that will stay 


positive. If we set R = X, \ and R are then calculated 

from 
D tT, cos (w/2) — 1 

R =\ = ————._ —_————__ (33) 
4rp sin (w/2) cos (w/2) + sin (w/2) 

Eq. (33) substituted into Eq. (31) gives 

xo = D/2rze X 
T4 COS (w/2) — 1 
rac0s(w/2)-1 gy 


scot w - — ta : — -. en { 
( 2 cos (w/2) [sin (w/2) + cos (w/2)] 
The length of the supersonic part of the nozzle can be 
found from 
D D *,————— ' 
= ——_ +- VM —1-—% (35) 
4sin (w/2) 2 
All equations derived above for the nozzle design may 
be divided by the diameter D of the nozzle, since D 
is involved in them only as a first power multiplication 
factor. As dimensionless quantities these functions may 
be used then for a nozzle of a certain Mach Number /, 
and of any diameter. 

The change of J/ and @ along a Mach line from point 
P to point K may be calculated from Eqs. (16), (6), and 
(7) as 
log F(6p) — log sin Op 


= M — Mp = ’ 
ld F fA 


resp. 


il (- — ‘) log F(@p) — log sin 6p 
«= _— _ = —_—_—— — ~ amen - 
: f 1VM,p* — 1+ cot 6 


SUMMARY OF STEPS IN ANALYTICAL DESIGN OF A THRER- 
DIMENSIONAL LAVAL NOZZLE WITH AXIAL SYMMETRY 


Given: Mach Number M, and diameter D. 

(1) Determine the fundamental constants Wz, tr, 
from Table 1 or by use of Eqs. (25) and (13a). Ac- 
curacy of one ten-thousandth of unity is recommended, 

(2) Determine the third constant w from 


o> 1/We = 4 = Wa 


Determine 7,4 and VM. — 1 from Table 1 for ¥, = 
‘We. 

(3) Calculate the coordinates of the transition 
curve from Eqs. (20) and (21) by varying V Mp? —] 





(4) Calculate xo by Eq. (34) and subtract it from all 
x values obtained in 3. 

(5) Calculate \ and R from Eq. (33) and determine 
the coordinates of point D (Fig. 4) by using Eq. (80). 


APPENDIX 


Eq. (20) may be obtained in a shorter and mathe- 
matically more general way than in the above analysis. 
The only assumption we make here is that the ordinate 
y of the Mach curve PK is uniquely given by 


y = ved( MV, Op) (36) 


M changes along the curve, but the way it changes will 
remain undetermined. We are concerned here only 
with the values of the function in Eq. (36) at the 
boundary points P and A, It is obvious that for point 
P 

o(Mp, Op) = 1 


For point K the value of the function ¢(.\/x, 6p) is 

given by the continuity equation as applied in Eq. (11). 
Substituting into Eq. (11) the line element 

é’(M, Op)d M* 

a °° % 

sin (@ + a) 


and performing the integration in the same manner as 
above, we have 
¢?(Mx, Op) — 1 = 
2(cos @p — cosw) [sin 6pV M,? — 1+ cos 4] 
sin? Op 
or 
o(M,, Op) = F(Op)/sin Op (37) 
* #’ is the first derivative of ¢@ with respect to the argument 
M. 


(Continued on page 188) 
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Recent N.A.C.A. Research on High Length- 
Beam Ratio Hulls 


ARTHUR W. CARTER" 


Nalional Advisory Committee for Aeronautics 


SUMMARY 

A series of related seaplane hulls having a wide variation in 
length-beam ratio has been investigated by the Langley Labora- 
tory of the N.A.C.A. An increase in hull length-beam ratio from 
6 to 15, when the length?-beam product is held constant, reduced 
the aerodynamic drag without appreciably affecting the hydro- 
dynamic qualities in smooth water, reduced the vertical accelera- 
tions and motions during landings in rough water, and reduced 
the structural weight required for a given load factor 

Extreme warping of the forebody planing bottom of length- 
beam ratio 15 resulted in large reductions in the maximum verti- 
cal and angular accelerations during landings in rough water and 
considerably improved the spray at low speeds. Warping the 
forebody and extending the afterbody resulted in further reduc- 
tions in the maximum accelerations in rough-water operation. 

The hull of high length-beam ratio, which has been shown to 
have acceptable hydrodynamic characteristics, is therefore of 
particular interest in the design of high-performance water-based 
aircraft, inasmuch as appreciable reductions in aerodynamic drag 
and hull weight may be obtained. Such high length-beam ratios 
combined with long afterbodies and warped forebodies are mark- 
edly superior to conventional length-beam ratios for rough-water 


operation. 


INTRODUCTION 


” VIEW OF THE REQUIREMENTS for increased range 
and increased speed in future flying-boat designs, 
the Langley Aeronautical Laboratory of the National 
Advisory Committee for Aeronautics has made an in- 
vestigation of the aerodynamic and hydrodynamic 
characteristics of flying-boat hulls as affected by hull di- 
mensions, and hull shape. The principal emphasis was 
placed on the reduction of aerodynamic drag without 
detriment to hydrodynamic performance. As part of 
the general investigation, a related series of hulls of 
varying length-beam ratio has been investigated. The 
series extends from the present conventional length- 
beam ratio of 6 to the high length-beam ratio of 15. 
(See Fig. 1.) The length Z used for determining length- 
beam ratio is the distance from the forward perpendic- 
ular (F.P.) to the sternpost (S.P.). The beam 3 is 
the maximum beam of the hull at the chine. 

The hulls were designed to have similar resistance 
and spray characteristics for the same gross weight and 
to be physically interchangeable on a hypothetical sea- 
plane design. Analysis of available information indi- 
cated that similar characteristics could be obtained if 
the hulls had the same value of L*b. The hulls, there- 

Presented at the Annual Summer Meeting, I.A.S., Los Angeles, 
July 14-16, 1948. 

* Hydrodynamics Division, Langley Aeronautical Laboratory. 


fore, became longer and narrower as the length-beam 
ratio was increased. The hull volume was reduced:28 
per cent and the frontal area was decreased 42 per cent 
when the length-beam ratio was increased from 6 to 15. 

The aerodynamic characteristics of the family of hulls 
were determined in the Langley 300 m.p.h. 7- by 10-ft. 
tunnel. An analysis of the change in structural weight 
with variation in length-beam ratio was made by the 
Langley structures research laboratory. The hydro- 
dynamic characteristics were determined in Langley 
tank No. 1. The present paper presents the results of 
these investigations with principal emphasis on the 
hydrodynamic qualities in both smooth and rough wa- 
ter. Effects of increasing forebody dead rise and after- 
body length on the hydrodynamic qualities of the hull 
with high length-beam ratio are included. 

The hypothetical seaplane design is a twin-engined 
propeller-driven flying boat having a design gross weight 
of 75,000 Ibs., a wing loading of 41.1 Ibs. per sq.ft., and 
a power loading of 11.5 Ibs. per b.hp. for take-off. The 
hydrodynamic qualities determined in the investiga- 
tion were the longitudinal stability during take-off and 
landing, the spray characteristics, the take-off and 
landing behavior in rough water, and the take-off per- 
formance, 


MODELS AND APPARATUS 


The parent model of the series was designed at Lang- 
ley tank No. | in accordance with well-established hy- 
drodynamic principles. Emphasis was placed on the 
necessity of accommodating the large changes in pro- 
portions required for the higher length-beam ratios, 
The overall dimensions of the hull of conventional 
length-beam ratio of 6 were based on the hull of a Navy 
twin-engined flying boat having good hydrodynamic 
qualities. 

The lines of the models of the series with high length- 
beam ratio were derived from those of the parent model 
with a minimum of change in other parameters. The 
procedure for deriving the portion of the hull below the 
chines consisted of three steps: 

(1) The spacing of the stations (transverse sections) 
was increased in order to arrive at the desired length- 
beam ratio. The forebody stations were moved parallel 
to the forebody keel at the step, and the afterbody sta- 
tions were moved parallel to the afterbody keel so that 
the angle between the keels was unchanged. The ratio 
of forebody to afterbody length was unchanged, and the 
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Fic. 1. Hulls having length-beam ratios of 6 and 15. 


Hat (constant dead rise) on the forebody bottom from 
the step forward increased in length with increase 
in forebody length. 

(2) The resulting lines were scaled down to arrive at 
the same value for L*) as that for the parent model. 
[he stations, below the chines, thus remained geomet- 
rically similar to those of the parent. 

3) The afterbody was raised vertically in order to 
irrive at the absolute depth of step of the parent (1.16 
in. for one-tenth-size models). The sternpost angle (the 
angle between the forebody keel and a line joining the 
forebody keel at the main step with the sternpost) thus 
became slightly lower as compared with that of the 
parent. 

The position of the step relative to the mean aerody- 
namic chord was held constant for all length-beam ra- 
tios. The maximum depth of hull and the length of 
hull aft of the step (afterbody plus length of tail exten- 
sion) were held constant in all the models in order to 
accommodate the same aerodynamic configuration. 
[he deck radii in all cases were held proportional to 
the change in maximum beam, and the fairing above the 
chines at the bow and tail was systematically altered to 
suit. As may be seen in Fig. 1, the high length-beam 
ratio of 15 remains closely comparable to the parent ex- 
cept for the difference legitimately introduced by the 
primary change in proportions. 

In the derivation of the series, the usual criterion for 
holding the bottom area constant (constant length-beam 
product) also was examined. As shown by extensive 
tank research, this criterion results in definitely im- 
proved hydrodynamic spray characteristics and resist- 
ance with increase in length-beam ratio. No significant 
improvement in aerodynamic or structural character- 
istics, however, could be anticipated for such a series. 
Furthermore, the resulting forms for the extreme 
length-beam ratios were obviously not interchangeable 
on the hypothetical seaplane designed around the par- 
ent hull. This criterion was, therefore, not pursued 
further. 

The wind-tunnel models! were one-tenth size and had 
length-beam ratios of 6, 9, 12, and 15. The models were 
mounted on a rectangular wing that spanned the tunnel 


test section. On all models, the wing had a 20-in. chord 
and was set at an angle of incidence of 4° to the keel of 
the hull. The hulls and wing were of laminated-wood 
construction and were finished aerodynamically smooth 
with pigmented varnish. 

The form, size, and relative locations of the aerody- 
namic surfaces of the one-tenth-size powered dynamic 
models used in the hydrodynamic investigation corre- 
sponded to those of the Navy flying boat. For conven- 
ience in making changes to the afterbodies, the fairing 
aft of the sternpost of the wind-tunnel models was 
omitted from the tank models and a slight modification 
was made to the sides of the afterbodies above the chine. 
These changes have a negligible effect on the hydrody- 
namic characteristics. 

The models were powered with three-bladed metal 
propellers driven by two variable-frequency motors. 
Slats were attached to the leading edge of the wing in 
order to delay the stall to an angle of attack more nearly 
equal to that of the full-size airplane. 

The hydrodynamic investigation was made in Lang- 
ley tank No. 1, which has been described by Truscott.? 
The apparatus used for the towing of powered dynamic 
models has been described by Olson and Land.* The 
setup of the model of high length-beam ratio on the 
towing carriage is shown as Fig. 2. The models were 
free to trim about the pivot, which was located at the 
center of gravity, and were free to move vertically but 
were restrained in roll and yaw. For self-propelled 
tests in waves they were also free to move fore and aft 
in order to absorb the horizontal accelerations intro- 


duced by the impacts. 


AERODYNAMIC CHARACTERISTICS 


The results of the wind-tunnel investigation in the 
Langley 300 m.p.h. 7- by 10-ft. tunnel to determine the 
effect of length-beam ratio on the aerodynamic charac- 
teristics of the family of hulls have been reported by 


Yates and Riebe.' 

The results of the tests are presented as standard 
N.A.C.A. coefficients of forces and moments. Yawing- 
moment and pitching-moment coefficients are given 
about the 30 per cent chord point of the wing. 





Fic. 2. Setup of model on towing carriage. 
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The coefficients and symbols are defined as follows: 


Cp = drag coefficient (drag/gS) 
Co = pitching-moment coefficient (/qSz) 
Cr = yawing-moment coefficient (.V/gSb) 
M = pitching moment, ft.lbs. 
N = yawing moment, ft.lbs. 
q = free-stream dynamic pressure, Ibs. per sq.ft. (1/29 V?) 
S = wing area, sq.ft. 
C = mean aerodynamic chord of wing, ft. 
b = wing span, ft. 
V = air velocity, ft. per sec 
p = mass density of air, slugs per cu.ft. 
a = angle of attack of hull base line, deg. 
y = angle of yaw, deg. 
Comin = minimum drag coefficient 
ma = longitudinal-stability parameter (0C,/0a) 
Cny = directional-stability parameter (0C,/dy~) 


The effects of length-beam ratio on drag and on the 
stability parameters C,,, and C,, are summarized in 
Fig. 3. The minimum drag coefficient for most condi- 
tions occurred in the angle-of-attack range between 2° 
and 3°. With transition fixed, the minimum drag co- 
efficient for length-beam ratio of 9 was less by a value of 
0.0009 (12 per cent) than the minimum drag coefficient 
for length-beam ratio of 6. Smalier reductions in mini- 
mum drag coefficient, 0.0007 and 0.0006, occurred when 
length-beam ratio was extended from 9 to 12 and from 
12 to 15, respectively. The overall reduction for an 
extension of length-beam ratio from 6 to 15 was 0.0022, 
a reduction of 29 per cent. 

Increasing length-beam ratio decreased the value of 
the longitudinal-stability parameter C,,,. This favor- 
able effect corresponded to a rearward aerodynamic- 
center shift of the flying boat of 2.5 per cent mean aero- 


dynamic chord when the length-beam ratio was 
changed from 6 to 15. The effect on the value of the di- 
rectional-stability parameter C,,, was adverse because of 
the increase of hull length and side area ahead of the 
center of moments at high length-beam ratios. 


STRUCTURAL CONSIDERATIONS 


An estimate of the relative weights of the hulls of the 
series has been made. This estimate was prepared by 
considering loading conditions in landing and assuming 
that components such as bottom, side, and deck plating 
have the same thickness at the same relative stations 
The variation in percentage reduction in hull weight 
with length-beam ratio is shown in Fig. 4. Increasing 
length-beam ratio from 6 to 15 resulted in a reduction of 
8.4 per cent in structural weight of the hull. 

The weight reduction was computed with the assump- 
tion that all hulls would be designed for the same load 
factor. The most recent Navy specifications do not 
make allowance for the fact that there will be some re- 
duction in maximum impact accelerations with high 
length-beam ratios. In current designs of seaplane 
hulls, the maximum acceleration of the impact will 
generally be reached before the chines of the hull sub- 
merge. If narrow hulls of high length-beam ratios are 
used, the chines will submerge before the maximum 
acceleration has been reached. Seaplane impact theory 
indicates that with chine immersion the maximum im- 
pact acceleration will be reduced. The estimated re- 
duction in hull structural weight is, therefore, conserva- 
tive. 


HYDRODYNAMIC QUALITIES 


In order to realize the established reductions in aero- 
dynamic drag and structural weight of the hulls of high 
length-beam ratio in the design of high-performance 
flying boats, it is necessary that the hulls have ac 
ceptable hydrodynamic qualities. The effect of length. 
beam ratio on the principal hydrodynamic qualities 
has been determined by an extensive tank investiga 
tion. 
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The hydrodynamic qualities of interest in the normal 
operation of seaplanes have been described by Parkin- 
son.‘ The qualities investigated are as follows: 

(1) Longitudinal stability. 

(a) Trim limits of stability. 
(b) Center-of-gravity limits of stability. 
(c) Landing stability. 
(2) Seaworthiness. 
(a) Spray. 
(b) Accelerations and motions in rough water. 
(3) Take-off performance. 


These qualities are discussed briefly in the order 
named, 


Trim Limits of Stability 

The trim limits of stability define the ranges of trim 
and speed in which porpoising occurs. Trim is defined 
as the angle between the forebody keel at the step and 
the horizontal. These trim limits consist of a lower 
limit; an upper limit, increasing trim; and an upper 
limit, decreasing trim. If the trim is decreased below the 
lower limit, low-angle porpoising, associated with the 
forebody alone, occurs. If the trim is increased above 
the upper limit, increasing trim, high-angle porpoising, 
associated with both the forebody and the afterbody, 
occurs. Once high-angle porpoising has been estab- 
lished, stability can be regained by decreasing the trim 
to some angle less than that at which the upper-limit 
porpoising first appeared. 

If the stable range of trim between the lower and up- 
per trim limits is known, stability during take-off be- 
comes a problem of balancing the aerodynamic and hy- 
drodynamic moments so that the trim track of the sea- 
plane (curve of trim against speed) lies within this stable 
range. 

The trim limits were determined at constant speeds 
by use of the methods described by Olson and Land.* 
In order to obtain sufficient control moment to trim the 
model to the trim limits, the lower limit was determined 
at forward positions of the center of gravity and the 
upper trim limits were determined at after positions of 
the center of gravity. 


Center-of-Gravity Limits of Stability 


The longitudinal position of the center of gravity, as 
in the case of aerodynamic stability, is a convenient 
parameter to define the stable range of pitching mo- 
ments. The range of positions of the center of gravity 
at which take-off stability is acceptable is of interest in 
determining the position of the step. By proper loca- 
tion of the step, both the hydrodynamic and the aerody- 
namic requirements for stability and control can be 
satisfied for the same range of position of the center of 
gravity. 

For a given elevator and flap deflection, the practical 
center-of-gravity limit is usually defined as that position 
of the center of gravity at which the amplitude of por- 


poising becomes 2°. Lower-limit porpoising, defining 
the forward center-of-gravity limit, usually occurs at 
intermediate planning speeds where the trim track in- 
tersects the lower limit and then re-enters the stable 
range, after which the instability damps out. Upper- 
limit porpoising, defining the after center-of-gravity 
limit, usually occurs at speeds near take-off where the 
trim track crosses the upper limit, increasing trim, and 
continues above it until the seaplane is air-borne. 

The center-of-gravity limits of stability were deter- 
mined by making accelerated runs to take-off with 
fixed elevators, full thrust, and a constant rate of ac- 
celeration of 1 ft. per sec. per sec. Trim, rise, and am- 
plitude of porpoising were continuously recorded during 
the accelerated run. A sufficient number of center-of- 
gravity positions and elevator deflections were investi- 
gated to cover the normal operating range and to define 
the center-of-gravity limits of stability. 


Landing Stability 


The hydrodynamic longitudinal stability of a sea- 
plane applies to both take-offs and landings, but the 
maneuvers differ in detail and the landing stability is 
best treated as a separate quality. The landing stabil- 
ity becomes of particular importance when the hull 
tends to leave the water in a succession of skips below 
flying speed when the seaplane is not under complete 
control. Skipping is primarily a function of landing 
speed and trim but is also influenced by the approach 
technique and the vertical speed. 

The landing stability was investigated by trimming 
the model in the air to the desired landing trim, at a 
speed slightly above flying speed. The towing carriage 
was then decelerated at a uniform rate of 2 ft. per sec. 
per sec., which allowed the model to glide onto the water 
and simulate an actual landing. The contact trims and 
behavior on landing were observed visually, and trim 
and rise were continuously recorded throughout the 
landing run. The landings were made with one-half of 
full thrust used during the take-off runs and with the 
center of gravity located at 32 per cent mean aerody- 
namic chord. 


Spray 


Spray is of importance in the operation of seaplanes 
when it obscures vision, inflicts physical damage to 
structural components, causes instability, or delays 
take-off by reducing the power of the engines. Of par- 
ticular importance in multiengined configurations is 
spray entering the propellers or striking the flaps. 

The range of speed over which spray entered the pro- 
pellers or struck the flaps was determined at constant 
speeds for gross loads from 55 to 125 per cent of the de- 
sign gross load. The minimum loads for spray wetting 
the propellers and the flaps were determined. Ob- 
servations were made of the spray striking the horizon- 
tal tail. 
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Accelerations and Motions in Rough Water 


The most severe service conditions for a seaplane are 
the rough-water take-off and landing. From an overall 
standpoint, the qualities of most interest are the verti- 
cal and angular impact accelerations and the motions 
in trim and rise. The accelerations are measures of the 
load factors for structures supporting concentrated 
masses, and the maximum trims are indicative of the 
extent of dangerous operation above the stall angle, usu- 
ally below flying speed and without lateral control. 

The landings were made with power on and with the 
thrust adjusted so that the model upon initial contact 
with a wave was approximately a free body. An accel- 
erometer mounted on the towing staff of the model 
measured the vertical accelerations. Two accelerom- 
eters were used to measure the angular accelerations. 
These accelerometers, which were mounted 1 ft. apart 
vertically, were located within the model in such a man- 
ner that their centers of gravity were in line with the 
center of gravity of the model. 

Results of several tests in rough water have shown 
that, except at dangerously low trims, there was no 
appreciable effect of landing trim on either the variation 
of trim during the landing run or the maximum acceler- 
ations. Consequently, all landings were made at ap- 
proximately 8°. An electrically actuated trim brake 
that was attached to the towing staff fixed the trim of 
the model in the air and controlled the initial approach. 
The trim brake was automatically released when the hull 
contacted the water. In order to determine the part of 
the hull contacting the water, contacts were located at 
the sternpost, at the step, and at a point approximately 
40 per cent of the forebody length aft of the forward per- 
pendicular. Wave struts forward and aft of the model 
were used to record the wave profiles and to determine 
the length between crests. 

Waves were generated by a wave maker that consists 
of a swinging plate hinged at the bottom and driven by 
a connecting rod at the top of the plate. The to-and-fro 
motions generate waves that travel from the north end 
of the tank through the test section and into an area 
where they are dissipated by a beach. The desired 
height and length of waves are obtained by a suitable 
combination of stroke and frequency of the plate. 


Take-Off Performance 


The take-off performance is conveniently defined in 
terms of take-off time and distance. The take-off time 
and distance are functions of the longitudinal acceler- 
ating force and, hence, of the longitudinal acceleration 
a. This acceleration was derived from the excess 
thrust, which was measured directly. 

The excess thrust, with the model propellers produc- 
ing scale take-off thrust, was measured by means of a 
spring dynamometer. The maximum excess thrust was 
determined at each speed by varying the trim of the 
models (elevator deflection), except at speeds where 
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porpoising was encountered. Over this speed range the 
trim was increased to remain just above the lower limit. 

The take-off time was determined from the area under 
a curve of 1/a plotted against speed V; the take-off 
distance was determined from the area under a curve of 
V/a plotted against speed. 


Basic SERIES 


In order to establish broad trends for the several 
hydrodynamic qualities described, the results of the 
tank investigation of the basic series are confined to 
the high length-beam ratio of 15 and the conventional 
length-beam ratio of 6, which is used as a basis for com- 
parison. All data are presented as full-size values to 
aid in their interpretation in terms of actual operation 
of the prototype flying boat with either of the two 
hulls. The tests were made at the design gross weight 
corresponding to 75,000 Ibs., except for the spray tests 
which covered a range of gross weights. The flaps 
were deflected 20° for all the hydrodynamic tests. 


Trim Limits of Stability 


The trim limits of stability are presented in Fig. 5. 
The upper limit, increasing trim, was approximately 
the same for both length-beam ratios. At high speeds, 
the effect of length-beam ratio on the upper limit, de- 
creasing trim, was small. The lower trim limit for the 
high length-beam ratio was shifted bodily to higher 
speeds. This shift, approximately 15 m.p.h., appreci- 
ably reduced the range of stable trim between the 
lower limit and the upper limit, increasing trim. 

Between 52 and 70 m.p.h., the upper and lower trim 
limits for the length-beam ratio of 15 were close to- 
gether. When constant-speed tests were made in this 
speed range, porpoising of the model could build up to 
such a large amplitude that the model porpoised across 
both the upper and lower limits. Under such conditions, 
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recovery from this violent porpoising by use of the ele- 
vators was not possible. Although not present with the 
length-beam ratio of 6 of the series, similar behavior has 
been noted for other models of flying boats with hulls of 
conventional length-beam ratios and, therefore, cannot 
be attributed solely to the high length-beam ratio, 
During accelerated take-offs, this violent porpoising 
was encountered only at center-of-gravity positions that 
were definitely ahead of the forward center-of-gravity 
limit and therefore is not considered of too great sig- 
nificance in the evaluation of the longitudinal stability 
during take-off. 


Center-of-Gravity Limits of Stability 


The maximum amplitudes of porpoising encoun- 
tered during take-off with fixed elevators are plotted 
against position of the center of gravity in Fig. 6. The 
maximum amplitude is defined as the difference between 
the maximum and minimum trims of the greatest por- 
poising cycle that occurred during the take-off. 

As the center of gravity moved forward, the maxi- 
mum amplitude of lower-limit porpoising increased. 
At the high length-beam ratio, a forward movement of 
the center of gravity of 3 per cent mean aerodynamic 
chord caused the amplityde to increase approximately 
8.5°, and both the upper and lower trim limits were ex- 
ceeded. The comparable increase in amplitude of 
lower-limit porpoising for the length-beam ratio of 6 
was approximately 2.5°, and porpoising between the 
limits did not occur. 

At after positions of the center of gravity, the slopes 
of the curves for constant elevator deflection approach 
zero at the most after positions of the center of gravity 
for the length-beam ratio of 15; the oscillation in trim 
with this high length-beam ratio did not exceeed 2.5° 
at the most after position of the center of gravity. The 
longer afterbody apparently was effective in damping 
the oscillations. In comparison, the amplitude of up- 
per-limit porpoising for the length-beam ratio of 6 con- 
tinued to increase with further movement of the center 
of gravity aft. 

Fig. 7 is a plot of elevator deflection against center- 
of-gravity position at which the maximum amplitude 
of porpoising was 2°. When this arbitrary limit was 
used, the high length-beam ratio had a range of stable 
center-of-gravity position which was approximately 
2.5 per cent mean aerodynamic chord less than that for 
the low length-beam ratio. There was, however, a wide 
practicable range of position for satisfactory take-off (10 
per cent M.A.C. with elevators deflected —10°) of the 
hull having the length-beam ratio of 15. 

Inasmuch as the upper-limit porpoising of the high 
length-beam ratio was not violent and was of small am- 
plitude, the after limit of position of the center of grav- 
ity does not have the same significance as that of the 
low length-beam ratio. For instance, if 2.5° amplitude 
of porpoising were selected as the maximum allowable 
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amplitude, there would be no after limit of position of 
the center of gravity for the high length-beam ratio, 
whereas the after limit for the low length-beam ratio 
would be moved aft approximately 2 per cent mean aero- 
dynamic chord. 


Landing Stability 


The maximum and minimum values of the trim and 
rise of the flying boat at the greatest cycle of oscilla- 
tion during the landing run were obtained from time 
histories of the landings in smooth water and are plotted 
against trim at contact in Fig. 8. Thehull having the low 
length-beam ratio did not skip at any contact trim inves- 
tigated, which indicated that the depth of step (9 per cent 
beam) provided sufficient ventilation. At contact trims 
up to 7.4° (the sternpost angle) the amplitude of oscilla- 
tion, both in trim and rise was small and was of approxi- 
mately a constant amplitude. At trims between 7.4° 
and 10°, the sternpost entered the water first and 
caused a slight increase in the amplitude of the trim os- 
cillation which damped out after one or two cycles. 
At contact trims above 10°, porpoising was encountered 
during the high-speed portion of the runout. This por- 
poising was the result of the trim remaining above the 
upper limit under the fixed elevator condition of the test 
and several cycles of upper-limit porpoising occurred 
before a trim within the stable range was attained. It 
would not necessarily occur in practice if the pilot let 
the trim decrease after contact. 

Similar trends in the landing behavior were noted for 
the hull of high length-beam ratio. This hull with the 
same absolute depth of step (16.5 per cent beam) also 
had no tendency to skip during landing. Porpoising at 
contact trims above 10°, however, was more severe than 
for the hull with the low length-beam ratio. This in- 
crease in violence of porpoising was attributed to the 
narrow range of stable trims between the trim limits, a 
relatively small oscillation in trim being sufficient to 
cause the trim to cross both the upper and lower trim 
limits. In normal operations, however, landings at 
trims greater than 10° would be avoided because of the 
danger of reaching a stalled attitude. 


Spray 

lhe range of speed over which spray entered the pro- 
pellers and struck the flaps is plotted against per cent of 
design gross load in Fig. 9. With the length-beam ratio 
of 6, the heavy, blister spray entered the propellers or 
struck the flaps at a lower gross load than with the 
length-beam ratio of 15. As the gross load was in- 
creased, the speed ranges over which the heavy spray 
entered the propellers or struck the flaps became 
slightly greater for the high length-beam ratio. At the 
design gross load, the spray entering the propellers or 
striking the flaps was approximately the same for both 
length-beam ratios. Based on the observation of the 
spray characteristics of a large number of models of suc- 
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Fic. 9 (top). Spray in smooth water. Fic. 10 (boltom). Spray 
in rough water. 


cessful conventional flying boats, the amount of spray 
entering the propellers or striking the flaps at the design 
gross load was considered acceptable. 


The horizontal tail of the low length-beam ratio was 
relatively clear of spray. The spray from the forebody 
of the high length-beam ratio hull struck the horizontal 
tail surfaces at high speeds. This spray might necessi- 
tate raising the horizontal tail if the spray loads became 
severe. They did not, however, damage the paper cov- 
ering of the tail surfaces of the model. 

The range of speed over which spray entered the pro- 
pellers in waves, 2 ft. high and 110 ft. long, is shown in 
Fig. 10. It was not possible to distinguish between 
light spray and heavy blister spray and, therefore, the 
comparison was made with the light-spray range of 
speeds in smooth water. The hull with the low length- 
beam ratio tended to ride over the tops of the waves, 
and the range of speed and load over which any spray 
entered the propellers was reduced for this particular 
wave. The hull with the high length-beam ratio tended 
to cut through the tops of the waves, however, and the 
range of speed and load was increased when compared 
with the range for smooth water. 


Accelerations in Rough Water 


A photograph of a typical record of a landing in 4-ft. 
waves is shown in Fig. 11. Of particular significance is 
the record of vertical accelerations where it is shown 
that the maximum acceleration occurred at an impact 
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Fic. 11. Typical record of landing in 4-ft. waves. 











VERTICAL 
ACCELERATION AT 
INITIAL CONTACT, g L/b 
15 
a § onan 
4 a 
iim 
0-4 T T T T T T — 
80 160 240 320 400 


WAVE LENGTH, FT 


MAXIMUM VERTICAL 
ACCELERATION, g 











L/b 
16- 15 
12-4 
8-4 
a 
Oo T T T T T T T ” 
80 i160 240 320 400 


WAVE LENGTH, FT 


Fic. 12 (top). Vertical accelerations at initial contact in 4-ft. 
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Fic. 14 (top). Effect of wave height on maximum vertical 
accelerations. Length-beam ratio, 15. Fic. 15 (bottom). An- 
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subsequent to the initial impact. The sinking speeds 
for the initial landing approach ranged from 125 to 325 
ft. per min., whereas, sinking speeds at the maximum 
vertical accelerations ranged from 300 to 1,400 ft. per 
min. In general, the sinking speeds at maximum verti- 
cal acceleration for the low length-beam ratio were 
greater than those of the high length-beam ratio. 

The variation of vertical acceleration at initial con- 
tact with wave length is shown in Fig. 12. The vertical 
accelerations at initial contact were approximately 45 
per cent less at the long wave lengths than at the short 
wave lengths. 

Che variation of maximum vertical acceleration with 
wave length is shown in Fig. 13. A peak in the curve of 
maximum vertical acceleration occurred at intermediate 
wave lengths. At the longer wave lengths, the accelera- 
tions were approximately 50 per cent less than the 
wcelerations at this peak. 

The increase in length-beam ratio from 6 to 15 caused 
negligible changes in the vertical accelerations at initial 
contact. The peak maximum accelerations, however, 
were reduced approximately 25 per cent. The peak 
vertical accelerations for both the low and high length- 
beam ratios occurred at approximately the same wave 
length. 

The effect of wave height on maximum vertical 
iccelerations of the hull with high length-beam ratio is 
shown in Fig. 14. An increase in wave height from 2 
to 4 ft. increased the peak accelerations approximately 
45 per cent. When wave height was increased from 4 
to 6 ft., the peak accelerations remained approximately 
the same. As wave height was increased, the peak 
accelerations occurred at longer wave lengths. 

The variation of angular acceleration at initial con- 
tact with wave length is shown in Fig. 15. The angular 
accelerations at initial contact were approximately 60 
per cent less at the Jong wave lengths than at the short 
wave lengths. 

The variation of maximum angular acceleration with 
wave length in Fig. 16, indicated that a peak was 
reached in the positive angular accelerations (bow ro- 
tating upward) at intermediate wave lengths. At the 
longer wave lengths, the accelerations were reduced 
approximately 65 per cent. The negative angular ac- 
celerations (bow rotating downward) did not vary 
greatly with wave length. The negative accelerations 
occurred when a bow-down rotation was induced dur- 
ing landing on the sternpost. The accelerations be- 
came slightly less negative with increase in wave 
length. The variation of angular accelerations with 
wave height was similar to that noted for the ver- 
tical accelerations. 

Increase in length-beam ratio from 6 to 15 increased 
the angular accelerations at initial contact approxi- 
mately 35 per cent. The peak maximum angular ac- 
celerations, however, were increased only 20 per cent, 


and the negative angular accelerations were reduced 35 
per cent. 


Motions in Rough Water 


The motions of the seaplane in trim and rise have 
been noted as qualities of interest in rough-water land- 
ings. The variation of trim, at the greatest cycle of 
oscillation, with wave length is shown in Fig. 17; both 
maximum and minimum trims obtained during the 
cycle are plotted. The greatest cycles occurred in waves 
approximately 200 ft. in length. The variation of trim 
with wave length was small; a slight reduction in the 
trim cycle was obtained, however, at wave lengths both 
shorter and longer than 200 ft. 

The maximum trims for both length-beam ratios 
exceeded the stall angle. The maximum trim of the low 
length-beam ratio was 4.5° greater than that of the high 
length-beam ratio. The maximum change in trim for 
the high length-beam ratio was approximately 25 per 
cent less than that of the low length-beam ratio. These 
reductions in the trim motion and in the maximum trim 
would make landings in waves less hazardous for the 
hull with high length-beam ratio than with the low 
length-beam ratio. 

The variation of rise, at the greatest cycle of oscilla- 
tion, with wave length is shown in Fig. 18; both maxi- 
mum and minimum rise obtained during the cycle are 
plotted. The greatest cycles occurred in waves approxi- 
mately 240 ft. in length. The maximum rise is reduced 
somewhat at wave lengths both shorter and longer than 
240 ft. 

The maximum rise was reduced when the length- 
beam ratio was extended from 6 to 15. The maximum 
rise for the low length-beam ratio was not determined 
for waves between 160 and 250 ft. in length, inasmuch 
as the rise would be in excess of the stop used to limit 
the rise of the model. The minimum rise for both 
length-beam ratios was approximately the same. 

The change in trim and rise at the greatest cycle for 
wave heights from 2 to 6 ft. followed the same trends as 
noted in the maximum vertical accelerations. In 6-ft. 
waves, however, the low length-beam ratio hit the rise 
stop at all wave lengths below 400 ft., whereas the high 
length-beam ratio did not hit the stop at any wave 
length. The maximum rise of the high length-beam 
ratio was approximately the same in both 4- and 6-ft. 


waves. 


Take-Off Performance 

The variation of the longitudinal acceleration a with 
speed V as determined from tests of powered dynamic 
models of the two length-beam ratios is shown in Fig. 
19. Comparison of the accelerations indicates that the 
water resistance of the low and high length-beam ratios 
was approximately the same. The take-off time and 
distance for the low length-beam ratio were 22 sec. and 
1,600 ft., respectively, and for the high length-beam 








176 JOURNAL OF THE AERONAUTICAL SCIENCES—MARCH, 1949 





SPEED, MPH 
. SPEED TRIM, DEG 

20 20 

10 16 

. TRIM 12 

60 
8 
4 





RISE ,FT 
5 





T errrrrttT. Tra ™ 
e€eees&+s#mp#mseRerseewnt 2s & 
TIME, SEC 


Fic. 20. Typical record of take-off in 4-ft. waves. Length 
beam ratio, 6. 














SPEED, MPH 
0 TRIM, DEC 
J SPEED js 
204 
404 
60 
TRIM 
RISE, FT fe) 
0- -4 
-54 RISE 
-10- 
"Se Odd 
© 2 4 6 8 10 12 4 16 18 20 22 24 26 
TIME, SEC 
Fic. 21. Typical record of take-off in 4-ft. waves. Length-beam 
ratio, 15. 
SPEED, MPH 
0 


RISE, 


ie) 














tes 
ss A 2a SM sk 
> 46 8 «WO MM (OS OB (20 22 A 
TIME, SEC 
Fic. 22. Typical record of take-off in 2-ft. waves. Length- 
beam ratio, 15. 
@ 
| 
| 
} ‘ ——{—> a \s a 
} | 11+ | 
ot i Se ee ~ 64> - mace 


Fic. 23. Length-beam ratio of 15 with warped forebody. 


ratio, 21 sec. and 1,530 ft., respectively. The overall 
take-off performarce of the low-drag hull with high 
length-beam ratio was therefore approximately the same 
as that of the conventional hull. 

The take-off behavior of the two hulls in rough water 
was investigated by means of take-offs at a constant 
rate of acceleration of approximately 0.lg. The results 
are qualitative, but several points are of interest 
Although the trim cycles were large in 4-ft. waves, there 
was no tendency for the bow to digin. Most of the long 
forebody of the high length-beam ratio was required, 
however, to keep spray from coming over the bow. 
Tracings of typical records made during these take-offs 
in 4-ft. waves are shown in Figs. 20 and 21 for the length 
beam ratios of 6 and 15, respectively. Both hulls 
demonstrated a tendency to follow the waves in trim 
and rise at the lower speeds. The phase relationships 
of trim and rise are also of interest. The rise reached a 
maximum shortly before the trim reached a mini- 
mum. 

The hull of low length-beam ratio exceeded the stall 
angle near hump speed. <A few seconds later the hull 
came clear of the water and reached a stalled attitude 
Since flying speed had not been obtained, the hull fel! 
back into the water with an impact acceleration of 2.5 
Upon contact with a wave, the hull again bounced clear 
of the water and trimmed to a stalled attitude. Flying 
speed was obtained before the hull again entered the 
water. 

The hull of high length-beam ratio (Fig. 21) had take- 
off characteristics in waves superior to those of the hull 
of low length-beam ratio. The oscillations in trim and 
rise at low speeds were large but did not appear to be 
dangerous. At higher speeds the oscillations became 
small as the hull planed over the wave crests and a 
relatively stable take-off was made. 

Fig. 22 shows a take-off of the hull of high length- 
beam ratio in a 2-ft. wave. The oscjllations in rise were 
extremely small. The oscillations in trim were not great 
and the trim did not exceed the stall angle during the 
take-off run. A comparison of the take-offs (Figs. 21 
and 22) of the high length-beam ratio show the marked 
difference in the motions in 2- and 4-ft. waves. 

The take-off investigation in rough water indicated 
that the hull of high length-beam ratio was less likely to 
reach a dangerous attitude than was the hull of low 
length-beam ratio; the behavior of the hull of high 
length-beam ratio was generally less violent. 


WARPED FOREBODY 


Although the range of stable position of the center of 
gravity of length-beam ratio 15 was only slightly less 
than that of length-beam ratio 6, the range of stable 
trims was reduced appreciably. In an effort to increase 
the range of stable trims and to determine the effect of 
increasing this range on the positions of the center of 
gravity for stable take-off, extreme warping of the fore- 
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body planing bottom (progressive increase in dead rise 
from step to bow) was ineorporated in the hull of high 
length-beam ratio. Results of investigations by David- 
son and Locke’ have indicated that warping of the fore- 
body bottom lowered the lower limit without causing an 
appreciable change in the upper limit, thereby increas- 
ing the range of stabie trims. Other possible advantages 
of such an increase in forebody dead rise would be a 
reduction in spray height and a decrease in vertical 
accelerations during rough-water operations. 

The model used for the investigation of forebody 
warp was the same as that having a hull length-beam 
ratio of 15 previously described, with the exception of 
the forebody bottom. Lines of the hull are shown in 
Fig. 23. The dead rise at the step (20°) was the same 
as that of the basic model. From the step forward to 
the forward perpendicular, the dead rise was warped at 
the rate of approximately 7.5° per beam. In order to 
provide straight buttock lines and chines on the planing 
bottom from station 7 to the step (3.6 beams), the 
angle of dead rise was varied between those stations as 
the tangent of the angle. The variation of the angle of 
dead rise, exclusive of chine flare, with length of fore- 
body of the basic and warped hulls is shown in Fig. 24. 
The keel heights, chine half-breadths, and chine flare 
were the same for both forebodies. 


Longitudinal Stability 


The trim limits of stability are presented in Fig. 25. 
The upper limit, increasing trim, for the warped fore- 
body was practically the same as that for the basic 
forebody. At speeds near take-off, the differences in the 
upper limit, decreasing trim, were negligible. The low- 
speed peak of the lower limit occurred at approximately 
the same trim for both forebodies, but, for a given trim, 
the lower limit for the warped forebody occurred at 
much lower speeds. This shift to lower speeds increased 
the range of stable trims between the lower limit and the 
upper limit, increasing trim. 

As noted for the basic forebody, porpoising at con- 
stant speed could be allowed to build up to such a large 
amplitude that the model porpoised across both the 
upper and lower limits. This porpoising was less violent 
than that encountered with the basic forebody and oc- 
curred over a smaller speed range. As in the case of the 
basic forebody, during accelerated take-offs this large- 
amplitude porpoising was encountered only at center- 
of-gravity positions that were definitely ahead of the 
forward center-of-gravity limit. 

At forward positions of the center of gravity where 
lower-limit porpoising occurred during take-off, the 
behavior of the hull with the warped forebody was simi- 
lar to that of the basicforebody. The forward limit for 
the center of gravity is seen in Fig. 26 to have been 
shifted aft approximately 2 per cent mean aerodynamic 
chord by the forebody warp. A more forward position 
of the step would therefore be required for the warped 
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Fic. 26. Center-of-gravity limits of stability. 


forebody in order to make the forward limit coincide 
with that of the basic forebody. 

At after positions of the center of gravity where upper- 
limit porpoising occurred during take-off, the maximum 
porpoising amplitude with the warped forebody did not 
exceed 3.5°, whereas this amplitude with the basic fore- 
body did not exceed 2.5°. A practical definition of the 
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after limit for either forebody therefore becomes difi- 
cult, and, inasmuch as this porpoising did not diverge 
to large amplitudes, an after limit might be considered 
nonexistent. The aft limits of the center of gravity, 
accordingly, have not been shown in Fig. 26. 

The maximum and minimum values of trim and rise 
of the hull with the warped forebody at the greatest 
cycle of oscillation during landings in smooth water were 
determined. The hull did not skip or porpoise on land- 
ing at any contact trim investigated (3° to 14°). At 
contact trims up to 10°, the amplitudes of oscillation in 
trim and rise were approximately the same as those of 
the basic forebody. At contact trims above 10°, the 
amplitudes of oscillation in trim and rise of the warped 
forebody were much less than those of the basic fore- 
body, and the amplitude of oscillation in trim was ap- 
proximately constant at landing trims between 10° and 
14”, 


Spray 


The range of speed over which spray entered the 
propellers and struck the flaps is plotted against gross 
load in Fig. 27. At the design gross load no spray 
entered the propellers or struck the flaps. The gross 
load was increased approximately 25 per cent before the 
heavy blister spray entering the propellers was equiva- 
lent to the spray of the basic forebody at the design 
gross load. The range of light spray in the propellers 
was considerably reduced by warping the forebody 
bottom. 

The range of speed over which spray entered the 
propellers in oncoming waves, 2 ft. high and 110 ft 
long, was determined. Spray entering the propellers 
was generally less and occurred over a smaller speed 
range for the warped forebody than for the basic fore- 
body. 


Accelerations and Motions in Rough Water 


Warping the forebody produced a negligible effect 
upon the vertical accelerations at initial contact. The 
maximum vertical accelerations, which are plotted 
against wave length in Fig. 28, were greatly reduced. 
The peak accelerations occurred at approximately the 
same wave length for both forebodies, but warping the 
forebody resulted in a reduction in maximum vertical 
accelerations at the peak of approximately 35 per cent 
At the long wave .lengths, the maximum accelera- 
tions for the two forebodies were approximately the 
same. 

Warping the forebody reduced the angular accelera- 
tions at initial contact approximately 27 per cent and 
resulted in a reduction in the peak maximum angular 
accelerations (Fig. 29) of approximately 50 per cent. 
The angular accelerations were, therefore, considerably 
less than those obtained for the hull with a length-beam 
ratio of 6. The negative angular accelerations were 
increased slightly by warping the forebody. 
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Fic. 31. Length-beam ratio of 15 with extended afterbody. 


The variation of trim and rise at the greatest cycle of 
oscillation with wave length was determined. The 
oscillation in trim was not appreciably affected by 
warping the forebody, and the maximum change in trim 
was approximately the same for the two forebodies. 
Likewise, the change in rise for both forebodies was 
approximately the same. The maximum rise at wave 
lengths below 200 ft. was the same for both fore- 


bodies. 


Take-Off Performance 


Comparison of the longitudinal acceleration during 
take-off (Fig. 30) indicates that the water resistance was 
approximately the same for the two forebodies up to the 
hump but was slightly greater at high speeds for the 
warped forebody. The take-off time and distance for 
the hull with the warped forebody were 24 sec. and 
1,780 ft., respectively, compared with 21 sec. and 1,530 
ft., respectively, for the basic hull. The overall take- 
off performance of the hull with the warped forebody 
was therefore slightly inferior to that with the basic 
forebody. This increase in high-speed resistance is in 
agreement with the results presented by Bell and Wil- 
lis® on the effect of dead rise. 


EXTENDED AFTERBODY 


Although an increase in hull length-beam ratio from 
6 to 15 appreciably reduced the maximum vertical 
accelerations during landings in rough water, it was 
believed that a further improvement in rough-water 
characteristics would be attained by an increase in 
afterbody length. Unpublished results of other investi- 
gations have shown that an increase in afterbody length 
made possible significant reductions in vertical and 
angular accelerations and decreased the maximum trims 
and height of bounces above the waves. 

Accordingly, an investigation of the hydrodynamic 
qualities was made of the hull with length-beam ratio of 
15 with an extended afterbody. Lines of the hull are 
shown in Fig. 31. The length of the basic afterbody was 
6.4 beams and the extended afterbody 9.3 beams, an 
increase of 45 per cent. The sections of the extended 
afterbody were the same as for the basic afterbody, but 
the distance between sections was increased in order to 
obtain the increase in afterbody length. 

In order to maintain approximately the same smooth- 
water landing stability for both the modified and the 
basic hulls, the ratio of the depth of step to the after- 
body length was held constant as the afterbody length 
was increased, Olson and Land.’ On this basis the 


depth of step was 24 per cent beam for the extended 
afterbody. This resulted in the same sternpost angle 
(6.9°) for the two afterbodies. 


Longitudinal Stability 


The trim limits of stability are presented in Fig. 32. 
Both branches of the upper limit were lowered by the 
increase in length of afterbody. The low-speed peak of 
the lower limit was approximately 1.5° less than that of 
the basic afterbody; otherwise, the lower limits did not 
differ greatly. Below 65 m.p.h., the high trims corre- 
sponding to the upper limits could not be obtained with 
the available elevator moments because of the large 
bow-down hydrodynamic moments of the long after- 
body. The range of speed, over which recovery from 
porpoising across both lower and upper limits could not 
be obtained by use of the elevators, was greatly reduced 
by the increase in length. In general, the trim limits 
of stability were improved by increase in afterbody 
length. 

A plot of elevator deflection against center-of-gravity 
position at which the maximum amplitude of porpoising 
was 2° is presented in Fig. 33. The forward limit for 





























TRIM, 
DEG 
125 
i UNSTABLE 
104 Ms 
. S47 UPPER LIMIT 
x “4 INCREASING TRIM 
85 ee ie 
\ 
é \ x 
>UPPER LIMIT 
6 . KE DECREASING 
stasLe TRIM 
a4 UNSTABLE 
LOWER LIMIT~ an 
~ 
24 AFTERBODY 
EXTENDED —— 
BASIK------- NACA 
Oo T T T T T 1 
30 40 50 60 70 80 90 
SPEED, MPH 
ELEVATOR 
DEFLECTION, 
DEG AF TERBODY 
-30- EXTENDED ——— 
BASIC ---------- 
-25- 
4 
4 
-20- \ 
4 
“15-4 
“104 ~ STABLE 
UNSTABLE 
au 
NACA 
18) T 7 T T my 7 T T T 1 
20 22 24 2 28 30 32 34 36 38 40 


CENTER OF GRAVITY, PERCENT M.A.C. 
Trim limits of stability. Fic. 33 (bottom). Cen- 


Fic. 32 (top). 
ter-of-gravity limits of stability. 








180 JOURNAL OF THE AERONAUTICAL SCIENCES—MARCH, 1949 

















GROSS LOAD, 
PERCENT DESIGN 
GROSS LOAD 
SPRAY IN PROPELLERS SPRAY ON FLAPS 
120-1 LIGHT HEAVY ‘ LIGHT HEAVY 
1 t yo t 
W ; 
1004 ' { a 
1 
i 
i CLEAR | 
B04 CLEAR 1% Jc.eaR 4 1 
Vv + AFTERBODY Vv 
EXTENDED —— 
= ae es <=>" 
60 T T T T 7 T T T T , 
0 10 20 30 4 650 0 t0 20 30 40 50 


SPEED, MPH SPEED, MPH 


Fic. 34. Spray in smooth water. 


the extended afterbody has been moved slightly forward 
of the limit for the basic afterbody. At after positions 
of the center of gravity, the maximum amplitude of 
upper-limit porpoising never exceeded 1° as compared 
with 2'/,° for the basic model, and, therefore, the after 
limit might be considered nonexistent. Thus, the 
longitudinal stability during take-off was slightly im- 
proved by extending the length of afterbody. 

The model with the long afterbody showed no tend- 
ency to skip on landing in smooth water and the depth 
of step, therefore, was considered adequate. Some por- 
poising occurred during the landing runout, but the 
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Fic. 35 (top). Maximum vertical accelerations in 4-ft. waves. 
Fic. 36 (bottom). Maximum angular accelerations in 4-ft. waves. 


afterbody length provided sufficient damping to keep 
the amplitudes small. The smooth-water landing sta- 
bility was approximately the same as that of the basic 
model. 


Spray 


The range of speeds over which spray entered the 
propellers and struck the flaps is plotted against gross 
load in Fig. 34. The speed ranges over which spray 
entered the propellers were approximately the same for 
the two afterbodies. It was apparent from observa- 
tions, however, that the volume of heavy, blister spray 
striking the propellers was greater for the extended 
afterbody. This increase in volume was attributed to 
the fact that the hull with the extended afterbody 
trimmed approximately 2° lower than the hull with the 
basic afterbody over the range of speeds in which 
spray was encountered. With the extended afterbody, 
only light spray struck the flaps at the design gross 
load. 


Accelerations and Motions in Rough Water 


The variation of maximum vertical and maximum 
angular accelerations with wave length is shown in 
Figs. 35 and 36, respectively. The accelerations were 
reduced at all wave lengths by the increase in afterbody 
length. Extending the afterbody resulted in a reduction 
in the maximum vertical accelerations at the peak of 
approximately 35 per cent when compared with those 
for the basic afterbody. This reduction in vertical 
acceleration at the peak was the same as that obtained 
by warping the forebody (Fig. 28). Extending the 
afterbody was less effective in reducing the maximum 
angular acceleration at the peak than was warping the 
forebody (Fig. 29). The negative angular accelerations 
were increased slightly by increasing the afterbody 
length. 

The variation of trim and rise, at the greatest cycle 
of oscillation, with wave length is shown in Figs. 37 and 
38, respectively. The oscillation in trim was reduced 
approximately 20 per cent by the increase in afterbody 
length. The maximum trim did not exceed 13.5° com- 
pared with a maximum of 20.3° for the basic afterbody. 
The danger of reaching a stalled attitude is, therefore, 
considerably lessened for the extended afterbody than 
for the basic afterbody. The oscillation in rise was 
reduced approximately 40 per cent. The maximum rise 
of 19.8 ft. for the basic afterbody was reduced to 10.5 ft. 
by extending the afterbody. The minimum rise for the 
two afterbodies was approximately the same. 


The desirable effects of the long afterbody are associ- 
ated with the restraint imposed by the increased moment 
arm of the planing areanearthesternpost. Observations 
and records showed the influence of this restraint in 
reducing the maximum trim and thereby the height of 
bouncing and the associated accelerations that occurred 
during the landing run. 
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(bottom). Rise at greatest cycle in 4-ft. waves. 


Take-Off Performance 


The excess thrust for the hulls with the extended and 
with the basic afterbodies was approximately the same. 
The overall take-off performance of the two hulls, there- 
fore, did not differ greatly. 


WARPED FOREBODY WITH EXTENDED AFTERBODY 


Warping the forebody has been shown to be effective 
in reducing the spray entering the propellers and strik- 
ing the flaps. Inasmuch as the volume of spray entering 
the propellers was increased considerably when the 
extended afterbody was applied to the length-beam 
ratio of 15, it appeared desirable to combine the ex- 
tended afterbody with the warped forebody in an effort 
to realize the advantages of the extended afterbody 
with respect to the overall longitudinal take-off sta- 
bility without the detrimental effect of large quantities 
of spray in the propellers. It was believed that the 
accelerations and motions in rough water could be fur- 
ther reduced by the combination, inasmuch as the 
warped forebody and the extended afterbody had each 
been effective in reducing the accelerations and the 
maximum trim of the basic hull. 


Accordingly, the hydrodynamic qualities were deter- 
mined for the hull with a length-beam ratio of 15, modi- 
fied to include both the warped forebody and the ex- 
tended afterbody. Lines of the hull are shown in Fig. 
39. 


Longitudinal Stability 


The trim limits of stability are shown in Fig. 40. As 
would be expected, the upper limits were approximately 
the same as those for the hull having the basic forebody 
and extended afterbody. Likewise, the lower limit was 
approximately the same as that of the hull with the 
warped forebody and basic afterbody except at the low- 
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speed peak, which was approximately 1.5° lower for 
the combination of warped forebody and extended 
afterbody. The range of stable trim was considerably 
greater than that for the hull with basic length-beam 
ratio of 15. 

A plot of elevator deflection against center-of-gravity 
position at which the maximum amplitude of porpoising 
was 2° is presented in Fig. 41. The forward limit for the 
combination of warped forebody and extended after- 
body was approximately the same as that of the basic 
hull. There was no after limit up to 40 per cent mean 
aerodynamic chord; the maximum amplitude of upper- 
limit porpoising did not exceed 2°. The combination 
thus had a large stable range of position of the center of 
gravity for satisfactory take-off. 

The landing behavior of the hull with warped fore- 
body and extended afterbody was similar to that of the 
hull with the extended afterbody, and the smooth- 
water landing stability was approximately the same 
as that of the hull with the basic length-beam ratio 
of 15. 


Spray 


The range of speeds over which spray entered the 
propellers and struck the flaps is plotted against gross 
load in Fig. 42. At the design gross load, 75,000 Ibs., no 
spray entered the propellers or struck the flaps for the 
hull combining warped forebody and extended after- 
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body. The spray entering the propellers at overloads 
was similar to that for the hull having the warped fore- 
body and basic afterbody. Only light spray struck the 
flaps at loads up to 95,000 Ibs. The bow spray charac- 
teristics in both smooth and rough water were similar 
to those of the hull having the warped forebody and 
basic afterbody. 


Accelerations and Motions in Rough Water 


The variation of maximum vertical and maximum 
angular accelerations with wave length is shown in 
Figs. 43 and 44, respectively. The accelerations were 
reduced at all wave lengths by warping the forebody 
and extending the afterbody. At most wave lengths, 
the combination gave the minimum accelerations when 
compared with the other configurations having high 
length-beam ratio. The hull combining warped fore- 
body and extended afterbody resulted in a reduction in 
the maximum vertical accelerations at the peak of 
approximately 55 per cent when compared with the 
basic hull and approximately 67 per cent when com- 
pared with the conventional hull with length-beam 
ratio of 6. The combination resulted in a reduction in 
the maximum angular accelerations at the peak of 
approximately 60 per cent when compared with the 
basic hull and approximately 50 per cent when com- 
pared with the hull of length-beam ratio of 6. 


The variation of trim and rise, at the greatest cycle 
of oscillation, with wave length is shown in Figs. 45 and 


























MAXIMUM 
ANGULAR 
ACCELERATION, 
RAD/SEC? FOREBODY AF TERBODY 
165 WARPED EXTENDED——— 
BASIC BASIC ------- 
124 i 
bi 
~ 
84 he 
~ 
iy 
~ 
Oo 
°4 T T T T T T T 1 
80 120 160 200 240 280 320 360 400 
WAVE LENGTH, FT 
TRIM AT 
GREATEST 
CYCLE, DEG FOREBODY AF TERBODY 
30- WARPED EXTENDED —— 
BASIC BASIC ———— 
244 
164 see” ame aia 
84 y-MINIMUM 
— =—=——— == NACA 
0 , T T Tt T t Ls 1 
80 120 160 200 240 280 320 360 400 


WAVE LENGTH, FT 


Fic. 44 (top). Maximum angular accelerations in 4-ft. waves. 
Fic. 45 (bottom). Trim at greatest cycle in 4-ft. waves. 





erloads 
d fore- 
ick the 
‘harac- 
similar 
ly and 


imum 
wn in 
$ were 
ebody 
ngths, 
when 
high 
fore- 
ion in 
ak of 
h the 
com- 
beam 
on in 
ik of 
1 the 
com- 


cycle 
» and 


ves. 





HIGH LENGTH-BEAM RATIO HULLS 183 


46, respectively. The oscillation in trim was some- 
what less than that of the basic hull and approximately 
the same as that of the hull with basic forebody and 
extended afterbody. ‘The maximum trims were slightly 
grater than those of the hull with the basic forebody 
and extended afterbody but were several degrees below 
those of the basic hulJ. The maximum trimattained was 
15.5°. The oscillation in rise was less than that of the 
basic hull but was greater than that of the hull with 
the basic forebody and extended afterbody. 


Take-Off Performance 


The excess thrust for the hull with warped forebody 
and extended afterbody was approximately the same 
as for the hull with the warped forebody and basic after- 
body but was slightly less than that of the basic hull at 
high speeds. The take-off performance of the hull with 
warped forebody and extended afterbody was therefore 
slightly inferior to that of the hull with basic length- 
beam ratio of 15. The take-off behavior in waves was 
similar to that of the basic hull. 


WATERLOOPING 


The waterlooping characteristics of the series of 
length-beam ratio hulls were investigated by catapult- 
ing '/s9-size models having length-beam ratios of 6 and 
12 into Langley tank No. 2. This investigation indi- 
cated that high length-beam ratio had no adverse effect 
on the high-speed directional stability. 


CONCLUDING REMARKS 


The investigation of the effect of length-beam ratio 
on the hydrodynamic qualities of seaplanes leads to the 
general conclusion that high length-beam ratio hulls 
having low aerodynamic drag and reduced structural 
weight have acceptable hydrodynamic qualities and are 
markedly superior to conventional length-beam ratios 
for rough-water operation. 


The longitudinal stability during take-off and land- 
ing, the spray characteristics, and the take-off per- 
formance in smooth water of the basic hull having a 
length-beam ratio of 15. were satisfactory and did not 
differ greatly from the same qualities of the related hull 
having a length-beam ratio of 6. An increase in length- 
beam ratio from 6 to 15 appreciably reduced the maxi- 
mum vertical accelerations and motions during landings 
in rough water but resulted in an increase in angular 
accelerations. 


Extreme warping of the forebody planing bottom of 
the hull having a length-beam ratio of 15 increased the 
longitudinal stability during take-off and landing in 
smooth water and considerably reduced the spray 
entering the propellers and striking the flaps, although 
the take-off time and distance were slightly greater than 
those of the basic hull. Warping the forebody resulted 
in large reductions in the maximum vertical accelera- 
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Fic. 46. Rise at greatest cycle in 4-ft. waves. 


tions during landings in rough water but did not ap- 
preciably affect the motions in trim and rise. 

A 45 per cent increase in length of afterbody of the 
hull of high length-beam ratio resulted in a slight im- 
provement in the longitudinal stability during take-off 
but increased the volume of spray entering the propel- 
lers. Extending the afterbody appreciably reduced the 
maximum accelerations and motions during landings in 
rough water. 

Combining the warped forebody and extended after- 
body on the hull of high length-beam ratio made pos- 
sible the realization of the improved characteristics of 
the two previous modifications. The longitudinal sta- 
bility during take-off was increased when compared 
with the basic length-beam ratio of 15, and the bow 
spray characteristics in both smooth and rough water 
were similar to those of the hull having the warped fore- 
body. Warping the forebody and extending the after- 
body resulted in the minimum accelerations obtained 
during landings in rough water of any of the configura- 
tions with high length-beam ratio. The combination of 
warped forebody and extended afterbody on the length- 
beam ratio of 15 resulted in a 67 per cent reduction in 
the maximum vertical accelerations and a 50 per cent 
reduction in the maximum angular accelerations when 
compared with the conventional hull of length-beam 
ratio of 6. 
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An Empirical Study of Low Aspect Ratio 
Lifting 'Sarkaces with Pastiveine Regard 
to Planing Craft 
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SUMMARY 


This report presents the results of a brief empirical study of the 
lift characteristics in incompressible flow of low aspect ratio 
lifting surfaces. It is indicated that for aspect ratios less than 
ibout 2.0, the lift coefficient can be estimated with reasonable 
.ecuracy by 

Cr = nKa" 
K and n depend only on the aspect ratio, and 7 is primarily a 
function of the operating conditions. Having assumed 7 to be 
unity for rectangular plan-form airfoils, it is then found that 7 
is 2/8 for hulls making leeway and 1/2 for flat-bottom planing 
surfaces. 


INTRODUCTION 


M““ YEARS OF EXTENSIVE THEORETICAL WORK pro- 
duced a simple theory that satisfactorily explains 
the lifting characteristics of wings of conventional 
aspect ratio acting in an incompressible fluid. On the 
other hand, only a rather complex theory! is available 
at present for low aspect ratio wings. The solution of 
such a theory for routine engineering problems is con- 
sidered far too time-consuming for ordinary needs. 

In addition to wings in the ordinary sense, there are 
at least two other types of lifting surfaces of general 
interest. They are the planing surface and the yawed 
hull of a displacement vessel. Both of these lifting 
surfaces are almost invariably of low aspect ratio, and 
both operate in an incompressible fluid. From time 
to time the need arises for estimating the lift character- 
istics of such bodies without entering into elaborate 
calculations. Since no simple theory is available, a 
correlation of available data on the lifting character- 
istics of low aspect ratio surfaces operating under vari- 
ous conditions was attempted. 


NOTATION 


Throughout this report the following notation and nondimen- 
sional coefficients have been used: 


Normal force coefficient Cy = N/S(p/2)V? 
Lift coefficient Cr = L/S(p/2)V? 
Froude Number of planing craft Cy = V/V gb 
Aspect ratio 

Airfoils and planing surfaces A = 68/S 

Hulls A = 262/S 


Received June 23, 1948. 
* Aeronautical Research Engineer (Hydrodynamics), Re- 
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Angle of attack from angle of zero lift, 


deg. 
Airfoils and planing surfaces a 
Hulls a = Xcosé 
where 


N = normal force, Ibs. 
L = lift force, Ibs. 
S = wing area, planing surface wetted area, or projected 
lateral area of hulls, sq.ft. 
p = mass density of fluid exerting lift force, lbs. sec.* per ft.‘ 
V = speed, ft. per sec. 
b = span, beam, or draft, ft. 
g = acceleration of gravity, 32.2 ft. per sec.? 
X = leeway of hulls, measured from angle of zero lift, deg. 
= heel angle of hulls, deg. 


and where 


K = the lift coefficient when a is unity in Eq. (2) 
n = the dimensionless exponent of a in Eq. (2) 
n = ratio between the lift coefficient of a surface and the lift 


coefficient of rectangular plan-form airfoil having the 
same aspect ratio and angle of attack 
¢ = a planing efficiency factor used in Eq. (3) 


ANALYSIS AND DISCUSSION 


The available theory of Bollay? and others says that 
the normal torce coefficient of an airfoil of zero aspect 
ratio is defined by 
Cy = 2.00 sin? a (1) 
This means that Cy is determined by an infinite series, 
since the sine is merely a shorthand method of repre- 
senting one infinite series. It seems possible that air- 
foils of small but finite aspect ratio can also be repre- 
sented by somewhat similar infinite series. Infinite 
series of this type can usually be closely approximated 
by simple power functions for small ranges. 

This study was undertaken primarily to throw a 
little light on the lift characteristics at normal angles 
of attack, so that C, could be substituted for Cy as an 
extremely close approximation. It will be assumed, 
therefore, that the lift characteristics of low aspect ratio 
surfaces may be represented by a simple power func- 
tion of the form: 


Cc, = nK af (2) 


It remains merely to find the constants of Eq. (2) by 
empirical analysis of available data. 
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Airfoils 

Winter's data® on the lift characteristics of rectangu- 
lar flat plates have been replotted in Fig. 1. The 
values of the angles of attack shown in this figure are 
measured from an angle of zero lift determined by ex- 
trapolation. Values of A and m were graphically 
measured on Fig. 1 and then plotted as functions of 
the aspect ratio in Figs. 4 and 5. 


Zimmerman‘ measured the lift characteristics of 
Clark Y airfoils of three different plan forms. His re- 
sults were replotted in charts similar to Fig. 1, which 
have not been presented to conserve space. The values 
of K and m were again measured graphically from these 
auxiliary charts and then plotted in Figs.4and5. The 
exponents m determined for the airfoils of the three 
different plan forms were in good agreement with the 
values determined with the aid of Winter’s data. 
There appeared to be no occasion, therefore, for dis- 
tinguishing the two sets of data on Fig. 5. 


The values of K determined from Zimmerman’s data 
were also in good agreement with Winter’s data, with 
three exceptions. All of the rectangular plan-form air- 
foils tested by Zimmerman agreed well with Winter's 
results. The lowest aspect ratio wing with faired tips 
(A = 0.65) had a value of K, which was about three- 
quarters of the value that might have been expected 
from the mean curve of the other tests in Fig. 4. The 
two elliptical plan form airfoils had values of K, about 
five-eighths of what might have been estimated by 
Fig. 4. These two elliptical plan-form airfoils had 


aspect ratios of 3/4and 1. Withexception of these three 
cases, all the rest of Zimmerman’s results are plotted in 
Fig. 4, and, since they agree with Winter, no attempt 
was made to distinguish the two sets of data. 

The fact that these three airfoils of Zimmerman’s did 
not agree with Winter's results is not regarded too 
seriously. Eq. (2) can be used to estimate the aero- 
dynamic side forces developed by the flying boat hulls 
(A = 0.2) shown by Yates and Riebe,°® provided 7 is 
taken as about 0.6. The reduction of lift can probably 
be attributed to the fact that none of these lifting sur- 
faces is rectangular. Thus, if 7 is assumed to be unity 
for rectangular plan-form airfoils, there is reason to ex 
pect some loss in lifting efficiency of airfoils of other 
plan forms and small*aspect ratio. Hence, caution 
should be exercised in estimating the lift character- 
istics by Eq. (2) of low aspect ratio airfoils having other 
than rectangular plan forms. 


Planing Surfaces 


There is available a great wealth of data on flat 
bottom planing surfaces. By cross plotting some of 
Sottorf’s data,® lift at constant aspect ratio could be 
determined. The results are shown in Fig. 2, and the 
points were determined from cross plots and do not 
represent measurements. 

Values of K and m were measured in Fig. 2 and then 
plotted on Figs. 4 and 5, respectively. It will be noted 
in Fig. 5 that, at given aspect ratio, the values of m are 
approximately the same for airfoils and planing surfaces. 
When first tried, K for the planing surfaces fell uni- 
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formly below the values determined for the airfoils. 
[t was found that, if 7 of Eq. (2) was assumed to be 
0.56 for the planing surface data shown in Fig. 2, the 
values of K for the airfoils and planing surfaces were in 
good agreement when plotted on Fig. 4. 

Analysis of some additional data of Sottorf’s and 
some unpublished tests conducted at the Experimental 
Towing Tank indicated that 7 could vary from 0.56 to 
as low as 0.45. Actually, it is difficult to determine 7 
with any great precision. However, a good average 
value of n for flat-bottom planing surfaces appears to be 
1/2. Bollay? suggests that the lifting efficiency of flat- 
bottom planing surfaces should be about 44 per cent 
from theoretical considerations. It is thought that 
assuming 7 to 1/2 for flat-bottom planing surfaces is in 
good agreement with the facts and not far from Bollay’s 
suggestion. 

There is a rather large range of speeds in which Eq. 
2) will not work unless 7 is taken to be much larger 
than 1/2. The unpublished Experimental Towing 
Tank tests just referred to were conducted in the sum- 
mer of 1940 to extend Sottorf’s work and thus deter- 
mine the boundary of planing. It was found that the 
boundary, defined by » = '/2, was remarkably definite. 
To the left and below the cross-hatched curve in Fig. 6, 
n increases rapidly. In the region to the right of the 
curve, n may be assumed to be 1/2 with good accuracy. 
Why the boundary should be so sharp is not all clear 
but it is presumed that careful study would show that a 
pronounced change in the wave patterns takes place, 
perhaps to allow buoyant forces to become important. 

For the purpose of going a little further into some of 
the details of the lift characteristics of planing bodies, 
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ANGLE OF ATTACK, DEG 


Eq. (2) may be written: 
C. = 0.5 c K a” (3) 


¢ is assumed to be unity for flat-bottom planing sur- 
faces in the region to the right and above the cross- 
hatched curve in Fig. 6. Eq. (3) will allow the study 
of various form parameters on the lift characteristics of 
planing craft. Bollay* presented a curve to show the 
effect of deadrise on the lifting efficiency of planing 
surfaces. His curve has been transposed to Fig. 7, 
The only data readily available to check this curve are 
given by King and Hill’ and Parkinson® for two differ- 
ent hulls. Similar cross plotting was done for these 
two hulls as was performed for the flat-bottom planing 
surfaces. It was found that their values of m in Eq, 
(3) conformed with the others shown in Fig. 5. How- 
ever, K was determined from Fig. 4 and thus allowed 
determination of fin Eq. (3). §¢ did not vary appreci- 
ably over the range of aspect ratios from 0.5 to 2.0, as 
deduced from the data in reference 7. On the other 
hand, the values of ¢ determined from the data given 
in reference 8 varied rather erratically. However, the 
average values of ¢ found for each model were plotted 
in Fig. 7. The accuracy of the values of the deduced 
values of £ actually depend on the accuracy of the as- 
sumed values of K. It is, therefore, perhaps presump- 
tuous to recommend the straight line drawn in Fig. 7 
in preference to Bollay’s curve, but it is believed that 
this line is safer to use until more data establishes 
something better. 

During the handling of the data taken from refer- 
ence 7, it was found that the planing boundary of Fig. 
6 was shifted slightly to the right. The shift amounted 
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TABLE 1 
DATA ON HULLS FROM FILES OF EXPERIMENTAL TOWING TANK 





Model 
No. Class 2b2/S 6 n K n 
205 Barge 0.085 0 0.59 0.00361 1.30 
257 Barge 0.089 0 0.44 0.00279 1.52 
234 Barge 0.107 0 0.95 0.00683 1.38 
329 Ketch 0.298 20 0.685 0.00975 1.265 
131 | 0.570 10,20,30 0.562 0.0123 1.23 
264 12M 0.615 20 0.56 0.0129 1.22 
109 12M 0.628 0 0.725 0.0170 1.22 
195 6M 0.748 10 0.83 0.0214 1.21 
99 6M 0.826 30 0.61 0.0173 1.18 


Average, 0.662 


toabout C, = 0.5. Thisis not considered large enough 
to warrant an additional curve in Fig. 6. However, it 
does seem likely that the shift to the right of the plan- 
ing boundary would continue with further increases of 
the deadrise angle. 


Hulls 


The hull of a sailing yacht sailing close-hauled pro- 
duces lift to counteract the wind forces on the sails. 
Barges and other displacement vessels produce lift to 
counteract the centrifugal forces in steady turns. 
These lift forces are perpendicular to the instantaneous 
course and parallel to the water surface. 


Such data as were readily available in the files of the 
Experimental Towing Tank were analyzed in 1940 
(when the majority of the work reported herein was 
In order to avoid some confusion, only part of 
However, 


done). 
the available data has been plotted in Fig. 3. 
all of the information that was available in 1940 is 
shown in Table 1. It was found that, in order to have 
the values of m in Eq. (2) for hulls fall in line with the 
values of m for airfoils and planing surfaces, it was 
necessary to resort to an effective hydrodynamic aspect 
ratio, which is twice the geometric aspect ratio. It 
will be seen in Fig. 5 that, if the effective hydrodynamic 
aspect ratio is used for the hulls, one line will represent 
the values of » for the three types of lifting surfaces. 
The concept of hydrodynamic aspect ratio was origi- 
nally proposed by Davidson.® The present study may 
be taken as a partial confirmation of the necessity of a 
hydrodynamic aspect ratio for yawed hulls, which 
differs from the usual geometric definition of aspect 
ratio. 


The values of K in Eq. (2) for the hulls were found by 
a trial-and-error process somewhat similar to that used 
for the flat-bottom planing surfaces. However, since 
so much variation in 7 was found, the process was some- 
what less straightforward. As may be seen in Table 1, 
the average value of 7 was found to be about 2/3. This 
average value was then used to determine K. The 
scatter of the values of K for the hulls shown in Fig. 4 
is probably attributable to some inherent variability 
of the lifting efficiency of the hulls of displacement 
vessels in yaw. 
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There are large variations in shape among the hulls 
considered, even in a given class. Since, as already 
noted, nonrectangular plan-form airfoils of aspect ratios 
less than 1 had different lift characteristics than rec- 
tangular plan-form airfoils, it may be inferred that the 
profile shape ofa hullhas an effect on its lifting efficiency 
Further, the thickness of a hull, in an airfoil sense, un- 
doubtedly affects the lift. Finally, differences in the 
wave making characteristics of different hulls could 
affect either K or 7, and it would be difficult to decide 


which. 


Accuracy 


A reasonably good estimate of the accuracy, that may 
be expected from the use of Eq. (2), may be had by 
examination of Figs. 4 and 5. Considering only the 
airfoil data in Fig. 4, the deviation of A will be found 
to be about +7'/. per cent. It should be remembered 
that an error in K will introduce a constant percentage 
error regardless of the angle of attack. In Fig. 5, the 
deviation of the airfoil data on n is about +4 per cent. 
A small percentage error in » will introduce an increas- 
ing percentage error in the lift coefficient as the angle 
of attack is increased. 


Experience with Eq. (2) indicates that, for aspect 
ratios greater than about 1.0, errors greater than about 
= 10 per cent need not ordinarily be expected for angles 
of attack less than about 15°. For the smaller aspect 
ratios, the same sort of error can be expected up to an- 
gles of attack of about 25°. 
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In the case of flat-bottom planing surfaces and hulls, 
errors can be kept to a minimum by a judicious choice 
of the factors 7 and ¢. 

It is believed that there are numerous engineering 
problems connected with low aspect ratio lifting sur- 
faces in which the errors described above would be 
tolerated in exchange for the ease of calculation. 


CONCLUDING REMARKS 


On the basis of this empirical study, it is concluded 
that the lift of low aspect ratio airfoils, planing surfaces, 
and hulls may be represented by a simple power func- 
tion of the form: 


C; = 7 K a” 


where A and n are functions only of the aspect ratio 
and are shown graphically in Figs. 4 and 5. 7 is a 
function only of the type of lifting surface and, if 
issumed to be unity for airfoils, becomes 2/3 for hulls 
and 1/2 for flat-bottom planing surfaces. 

This simple relation should be a useful engineering 
tool for making quick estimates. Much work remains 
to be done to obtain reasonably accurate means for 
determining the lift of various kinds of planing surfaces 
ind different type of yawed hulls. It is hoped that this 
study may provide a foothold on which to expand. 
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Axially Symmetric Laval Nozzle 


(Continued from page 166) 


Eq. (37), substituted into Eq. (36) and by use of Eqs, 
(14), (15), and (19), gives 


y = [D/4 sin (w/2)]rp/reF (6p) (20) 


The exact analytical form of Eq. (36) is not essential 
for computing the transition curve AB. This appears 
to be strange, but it is explained by the fact that the 
mass flow equation was integrated along PK in a lin- 
earized form. 
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On the Hodograph Method 


¥, H. Kuo 

Assistant Professor of Aeronautical Engineering, Cornell University, 
Ithaca, N.Y. 

December 20, 1948 


5 gpa HAVE RECENTLY APPEARED two papers concerning the 
transonic flow of a compressible fluid about a closed cylin- 
jer, using the hodograph method, by M. J. Lighthill'! and T. M. 
Cherry.? As is well known, such a flow is generally two-sheeted, 
vith two branch points, in the hodograph plane. It is clear that 
1 Taylor series is never sufficient to represent a complete solu- 
tion in this situation. For this reason, the question of analytic 
continuation of a Taylor series always arises; it is in this connec- 
tion that Lighthill and Cherry have made important contribu- 
tions. 

The Lighthill-Cherry method of continuation is essentially 
based on the property that the hypergeometric function F,(7r) has 
simple poles at negative integers, so that it is accessible to a rep- 
resentation by a series of partial fractions. By means of this 
form of representation the authors are enabled to transform the 
original Taylor series into one with a larger radius of convergence. 
However, this requires the coefficients of the Taylor series to be 
A,T~"(r;) rather than A,F,~'(r:), where Aq is the coefficient of 
the Taylor series of the corresponding incompressible flow, T(r:) 


a known function of 
: ell. “) 
T = as ae 
2 Co? 


m the circle of convergence, which is always positive and 
nonvanishing in the subsonic range, and m an integer. Thus, 
the simplification in analytical continuation is brought about 
by the sacrifice of a closer resemblance to the incompressible 
dow on the circle of convergence. 

In the appendixes of both papers, a previous paper by Tsien 
ind Kuo* on the same subject is criticized and is, in fact, sus- 
pected of being incorrect, because the solution thereof does not 
agree with the new results. Both authors attempt to demon- 
strate the incorrectness by deriving a solution, for comparison, 
by their method. During the process, apparently, several basic 
facts have been overlooked, and, as a result, the criticism is not 
justified and should be clarified. 

a) To show that a solution, not satisfying his condition C 
reference 2, page 76), continued by the Tsien-Kuo method from 

' Lighthill, M. J., The Hodograph Transformation in Transonic Flow. 
lll, Proc, Roy. Soc. A 191, pp. 352-369, 1947. 

2 Cherry, T. M., Flow of a Compressible Fluid About a Cylinder, Proc. 
Roy. Soc, A 192, pp. 45-79, 1947. 

4 Tsien, H. S., and Kuo, Y. H., Two-Dimensional Irrotational Mixed 
Subsonic and Supersonic Flow of a Compressible Fluid and the Upper Criti- 
cal Mach Number, N.A.C.A. T.N. No. 995, 1946 


a Taylor series, diverges, Cherry considers as an example the 
function 


Wilw) = Yi-w+w 


where w = ge~** is the complex velocity. This function consists 
of two component parts: one is an entire function and the other 
is singular at w = 1. Therefore, Wo(w) has a branch point at 
w= 1. The proper expansion of this function is therefore 


Wow) = p A,w" +w 
0 
for |w| < land 
Wo(w) = i >> Anw-l*- C/O] 4 wy 
0 


for |ew| >1. Both series are absolutely and uniformly convergent 
in their respective domains. Instead of taking the correct ex- 
pansion shown above, Cherry, todemonstrate the supposed short- 
comings of the Tsien-Kuo procedure, forces a second expansion of 


the form: 
Wow) = i D> (Baw) -" + Cyw*-C/] 
0 


and then computes the coefficients B, and C, by the rule of Tsien 
and Kuo. Clearly, by the lack of the proper singularity, this series 
has no choice but to diverge, and it must be noted that this di- 
vergence has nothing to do with the condition C. A difficulty 
like this can always be avoided by sound judgment. 

(b) To show that the continued solution by Tsien and Kuo’s 
method is incorrect, both authors deduce a solution by taking the 
same Taylor series as Tsien and Kuo—namely, for g << 1— 

— F,(r) 
¥ - Pe es id 
1* (g, 3) »e Ang Fan) sin nd 


but use their method of continuation. Indeed, the continued 
solution is different from Tsien and Kuo’s result. The dis- 
crepancies, however, can be traced to an erroneous assumption, 
which seems to have been taken for granted by the authors. It 
must be remembered that the Lighthill-Cherry method of con- 
tinuation depends on the expansibility of F,(r)7~1"(71) as a series 
of partial fractions. In the case of F,(r)/F,n(7ri), the method 
simply fails. Asremarked by Cherry, the function F,(r)/F,(71) is 
now regular at negative integers. Hence, the poles of the func- 
tion are the zeros of F,(7:). But for r< (y — 1)/(y + 1) or in 
the subsonic range in which 7 always lies, the hypergeometric 
function F,(r) is nonoscillatory and is positive for all real m. 
Therefore, F,(7:) has no root in this interval.‘ Consequently, 
F,(r)/F (71) cannot be expanded as a series of partial fractions. 
In other words, Eqs. (53) and (54)! and (7.12) and y*? are fic- 
titious. 

It is undeniable that the Tsien-Kuo method might run into a 
divergent series if not handled properly. This is not a charac- 
teristic feature, and it may be remarked, by way of analogy, 
that the lack of physical consideration will often lead to improper 
integrals and divergence in many problems, for example, in wing 
theory. The same, of course, can be true in the hodograph 
method. Fundamentally, Tsien and Kuo’s method assumes 
nothing beyond the existence of a solution for a linear partial dif- 
ferential equation of two variables. If every requirement for a 
solution to exist is met, there should be no question as to whether 
the series or integral hasa meaning. At this stage of the develop- 
ment of the hodograph method for transonic flows, an important 
criterion to judge any procedure is the amount of labor involved 
in carrying out the numerical calculation. It is this aspect of 
the Tsien-Kuo procedure which was thought unsatisfactory and 
was improved in a later report in 1946.° 


4 Chaplygin, S., Gas Jet, N.A.C.A. T.N. No, 1063, 1944. 
&’ Kuo, Y¥. H., Two-Dimensional Irrotational Transonic Flows of a Com- 


pressible Fluid, N.A.C.A. T.N. No. 1445, 1948. 
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Impact Tube at Low Pressures—Errata 


Paul L. Chambré 
University of California, Berkeley, Calif. 
December 22, 1948 


A“* unfortunate error occurred in the article ‘‘The Theory of 
the Impact Tube at Low Pressures”’ (by P. L. Chambré and 
S. A. Schaaf, Journal of the Aeronautical Sciences, Vol. 15, No. 
12, pp. 735-737, December, 1948). 

In Eq. (1) on page 736, the exponent 2 is missing from the left- 


hand expression: 


Pepe ae ( oe 6446 6 Oe 6 0-6 } { La. cre © oe a (1) 
Pipi W(L/r) 


Aerodynamic Relations with Variable Specific 
Heat 


J. D. Schetzer 
Associate Professor, University of Michigan 
December 30, 1948 


F gee the last 2 years, several publications have appeared 
which give tables and charts for use in compressible flow 
calculations. These tables are quite generally based on a con- 
stant value of 1.4 for the ratio of the specific heats of air. While 
this is a good assumption for a moderate range of temperatures, 
it ceases to be accurate when the temperature variations become 
large. For high-speed flow, the difference between the static 
and stagnation temperatures is extremely large, and the constant 
y formulas on which the charts and tables are based are not 
valid. 

The authors of these publications are, of course, aware of this. 
However, because the charts and tables include values for the 
flow ratios in the very high Mach Number range, it is felt that a 
word of warning to the unwary user is in order. Two years ago, 
C. H. Lauritsen and the writer assessed the errors arising from 
the assumption of a constant specific heat ratio in the calculation 
of static to stagnation ratios and in the ratios of the various flow 
parameters across normal shock waves. Errors were checked for 
a stream static temperature of 520°R. The results are reported 
here. 

The usual formulas for the calculation of the ratio of the static 
to stagnation values of the pressure density and temperature 
based on y = 1.4 are: 


p/po = (1 + 0.2 M2)-88 (1) 
p/po = (1 + 0.2 M2)-25 (2) 
T/T) = (1+ 0.2 M2)-1 (3) 


Values of these ratios, which take into account the variation in 
the specific heats of air, were computed from the enthalpy tables 
given in Keenan & Kaye, Thermodynamic Properties of Air, John 
Wiley & Sons, Inc., New York, 1945. The correction factors 
that must be applied to Eqs. (1), (2), and (3) above are plotted 
versus Mach Number on Fig. 1. For Mach Numbers less than 2, 
the errors are less than 1 per cent. For Mach Numbers greater 
than 2, the errors rise rapidly. At a Mach Number of 6, the er- 
rors in pressure, density, and temperature ratio are 26, 34, and 
13 per cent, respectively. 

The usual formulas for the calculation of the flow ratios across 
a normal shock wave with y = 1.4 are given by: 


p:/pi = 1.167 M2? — 0.167 (4) 
T:/T; = (2+ 0.4 M;?) [(2.8 M,? a 0.4) /5.76 M;?} (5) 
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ratios across a normal shock. 


p/p’ = [2.4/(2.8 M,? — 0.4)]*5 [2.4 1,2/(2 + 0.4 M,2)]*! 


M2? = (5 + M,*)/(7 Mi? — 1) (7 


Subscripts “1”? and “2” refer, respectively, to upstream an( 
downstream conditions. Eq. (6) is the ratio of stagnation pres 
sures. From a graphical solution of the equations of continuity 
equilibrium, state, and energy as computed from, the enthalpy 
tables, true values of the above ratios may be found. The cor 
rection factors that must be applied to Eqs. (4) through (7) are 
plotted versus Mach Number on Fig. 2. For upstream Mach 
Numbers less than 2, the use of a constant value of 1.4 for the 
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specific heat ratio leads to an error of less than 1 per cent in the 
quantities considered. At an upstream Mach Number of 6, the 
errors in the pressure, density, temperature, and stagnation pres- 
sure ratios are 3, 16, 11, and 24 per cent, respectively. The error 
in the downstream value of the Mach Number is 5 per cent. 
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The Simulation of Combustion Models in 
Wind Tunnels 


Bruno W. Augenstein 

Assistant Professor, School of Aeronautics, Purdue University, 
Lafayette, Ind. 

January 7, 1949 


[= THE COURSE of some theoretical investigations for The 

Rand Corporation on problems of unsteady flow in ducts, 
the question arose as to what extent duct flow with burning can 
be simulated without heat addition in a wind tunnel (i.e., a cold 
test). The natural tool for such an investigation is the theory of 
‘haracteristics: if the differential equations are hyperbolic and 
the number of nonperiodic independent variables is two, general 
methods of integration exist, while if the number of nonperiodic 
independent variables is more than two, specialized methods are 
igain required.! In some cases, we can reduce elliptic equations 
to hyperbolic form by use of complex independent variables. 

In the problem at hand, the independent variables are x and £, 
space and time variables, respectively. If we denote the particle 
velocity by u, the speed of sound by a, the entropy by S, and the 
channel area by A, then the basic equations of one-dimensional 


flow are: 


2a a? | 
Momentum UUs + use + — az — ——— S:=0 
y=] (y- 1)Cp 


2u « 
Continuity auz + — az + —_ a; = 
7-1 y¥=- 1 
.: 2 y a — 
—— (A; + uAz) + — — (S; + uSz) = I 
A y¥-1C, 
Energy S:+ uSz = F(x, t) = Il j 


Here y = ratio of specific heats, Cp = specific heat at constant 
pressure, and F(x, ¢) is some known function expressing the rate 
of heat addition. Further, we suppose that we know some curve 
vin the (x, ¢) plane in the parametric form x = x(t) and that be- 


cause of initial or boundary conditions we know u, a, and S and 
their derivatives along y. Then we have a set of six equations 
in the partial derivatives; we write this set as: 











+m +—-a+0 £ | 
UUs ue az a — ——_~ X 
y¥-1 (y — 1)C, 
Sz + 0S;= 0 
2u 2 
Que +Ou; + @e + —az+058, + 
y¥-1 y¥-1 
OS; = I 
Our +0u; + 0a: + Oa +u Sz + 
— (2) 
Ss =I 
dx tz + dtu; + O0az + Oa +0OS; + | 
0 Ss = du | 
Ouz +O0u: + dxaz + dt as + 0S; + | 
OS: = da | 
Ours +0u, + 0az + Oa + dx Ss +f. 
dt Ss = dS } 


Then the theory of characteristics tells us that :! 

(A) Setting the determinant of the six by six array of coeffi- 
cients on the left-hand side of Eq. (2) equal to zero gives us the 
equations of the characteristic lines in the (x, f) plane. We have 
then: 














| 2a a? | | 
*s 1—0 - ——— 
2u 2 } 
a O -— - 0 0 
— y-ly-1 
000 0 u 1 
dx dt 0 0 0 0 | 
0 O dx dt 0 0 | 
000 0 dx dt | 
: : (3) 
<a. | 
— 0 | 
y¥-1 | 
2u 2 
= (udt— dx) la 0 — — 
| y-1 7-1 
0 dt O 0 
dx 0 dx dt | 
2 ' 
= — —— (udt — dx) [(udt — dx)? — (adt)?] | 
73-4 j 


Thus we have the three sets of characteristics: dx/dt = u, the 
particle path lines; and dx/dt = u + a, the wave front path lines, 
If u is'positive to the right, uw + a refers to a right moving wave 
and u — arefers toa left moving wave. 





1 Courant and Hilbert, Mathematische Physik, Chapt. 5, Vol. 2. 


B) If we make use of the fact that wz, 1, az, at, etc., must not be infinite on the characteristics, we obtain certain dynamic condition 


that must hold on the characteristics, the compatibility conditions. 





For example: 





2a a? | } 
0 1 0 — . 
y-1 (y—1C, | 
it eaten ae 0 | | 
y-ly-! ; 
1110 Oo 0 u 1 | 
|du dt O 0 0 o| | 
|\da 0 dx dt 0 | t (4) 
ldSO O 0 dx dt| 
—— ne : said ee = i 
— [2/(y — 1)] (udt — dx) [(udt — dx)? — (adt)?] 
2 J 2a a? 
——— (udt — dx) 4 du(udt — dx) + I a(dt)? — ——— dt da — ———_~ [II (dd)? -— dS dt] 
= 1 l ¥-1 (y— DC» | 


—[2/(y — 1)] (udt — dx) [(udt — dx)? — (adt)?] 


hus, vz is automatically finite on the particle path lines, and upon substituting udt — dx = +adt, 
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d/dtenar. = (d/dt) + (u * a) (d/dx) 


the fact that u, must be finite on the wave fronts gives 


d 2 [ 
oe ip as oe} =| St ———— * 
dtchar. y¥-1 (y -—1)C, 


yau ya 
a ae 
(y¥-0C,/ (y- DC, 


1949 


(d/dtenar. = total derivative along the wave front characteristic) 


Sy ¥ 7 (A; + uAz) ] (5) 


This is the general characteristic equation for one-dimensional flow if injection of matter is neglected; special cases of Eq. (5) have been 


given by various authors.?~¢ 


Note that the method of derivation given is perfectly general and can, for example, be used to derive 


the equations of steady, two-dimensional or axially symmetric nonisentropic flow. 


We are interested in some of the general conclusions that can 
be derived from Eq. (5) without calculation, and these will give 
us the conditions under which combusting models may be simu- 
lated in wind tunnels without actual heat addition. For ex 
ample [these results are obtained by considering the appropriat« 
terms in Eq. (5)]: 

(1) In steady flow, heat addition gives the same upstream ef- 
fects that throttling does; consequently, adding heat readjusts 
the upstream flow conditions in exactly the same way that de- 
creasing the throttle area does. 

(2) If we consider a constant area flow with a combustion zone 
we can simulate the dynamical effects of a pressure pulse travel- 
ing through the duct by the same pressure pulse traveling through 
an area change in the cold duct; thus, if the pulse travels from a 
cold zone to a hot zone, the hot zone can be simulated by an area 
increase in the cold duct. The pressure change and the velocity 
change at the temperature discontinuity cannot, however, be 
simultaneously satisfied by replacing the temperature discon- 
tinuity by an area change; it is not possible to simulate both 
the dynamic and kinetic effects of two regions of different entropy 
simultaneously by area changes. 

(3) The initial upstream effects of unsteady heat addition can 
be simulated by unsteady throttling, and the complicated 
“throttling program” can be derived from Eq. (5). 

(4) However, Eq. (5) shows that the complete dynamics of the 
transient state in going from one steady heat addition rate to 
another steady heat addition rate is not simulated by the throt- 
tling. In any actual problem the upstream and downstream re- 
flections of the initial pressure waves due to the unsteady heat 


? Kahane, A., and Lees, Lester, Unsteady One-Dimensional Flows with Heat 
Addition or Entropy Gradients, Journal of the Aeronautical Sciences, Vol. 15, 
No. 11, p. 665, November, 1948. 

* Guderley, G., ZWB 1744, October, 1942. 

* Sauer, R., F-TS-770-RE, Wright Field. 

§ Schultz-Grunow, F., Ing. Archiv., 14, 1943. 

* Love, Pidduck, Phil. Trans. Roy. Soc., 222, 1922. 








addition will eventually return to the flame front, and the transi- 
tion through the temperature discontinuity cannot then be simu- 
lated in the cold model by a throttle. Thus, in the actual model, 
increasing the rate of heat addition sends out upstream and 
downstream compression waves; if, now, there is a contraction 
at the downstream end of the combustion zone, the downstream 
compression waves will be reflected from the contraction and the 
reflected upstream waves will be amplified both by the reflection 
and by the travel through the upstream flame front. Thus, a 
second compression pulse will travel upstream in the combustion 
model, resulting in larger upstream pressure fluctuations than 
those due to the initial upstream compression waves, and these 
secondary upstream compression pulses will become stronger 
when either the heat addition rate is increased or the combustion 
zone exit area is decreased; these secondary waves cannot be 
simulated by the unsteady throttling. These facts, together 
with the possible instability of the flame front itself, indicate 
that the steady state associated with the new heat addition rate 
may not be attained at all, and in any event the times to reach 
equilibrium will! be very different in the hot and cold tests; it is 
easier to initiate slow instability with actual heat addition. Op- 
eration of several valves would give somewhat better simu: 
lation, but only in a limited way. 

In general, therefore, we must settle with actual combustion 
such questions as coupling or resonance oscillations between the 
combustion chamber and the rest of the duct, times to reach new 
steady flows, complete transient upstream effects of unsteady 
heat addition, etc. The use of a single variable valve to simulate 
unsteady heat addition (or rough burning) gives results of limited 
value. 

Eq. (5) can be solved by setting up difference equations be 
tween the physical boundaries and the characteristics. Some 
applications are: unsteady flow in ducts caused by time varia- 
tion of area (wind-tunnel diffusors), processes involving both 
area variations and heat addition rates varying with time (ex- 
plosion turbines), etc. 
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